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1.  INTRODUCTION 


The  fundamental  goals  of  this  Predoctoral  Traineeship  were  three-fold:  1)  to 
provide  me  with  an  opportunity  to  continue  to  learn  and  apply  the  state  of  the  art  NMR 
technology  and  radiotherapy  techniques  to  cancer  diagnosis  and  treatment;  2)  to  provide 
me  with  a  solid  and  extensive  training,  and  valuable  experience  in  a  modem  NMR 
laboratory  for  a  career  as  a  clinical  medical  physicist  and  a  breast  cancer  research 
scientist;  and  3)  to  develop  and  investigate  a  non-invasive  technique  of  measuring  oxygen 
tension  (p02 )  in  breast  cancers  in  an  animal  model  based  on  19F  MRI  of 
hexafluorobenzene  (HFB),  now  called  the  FREDOM  (Fluorocarbon  Relaxometry  using 
Echo  planar  imaging  for  Dynamic  Oxygen  Mapping).  In  this  report,  I  will  summarize  the 
highlights  of  my  past  three-year's  training  and  research  as  originally  proposed  in  my 
Predoctoral  Traineeship  application. 


2.  TRAINING  ACCOMPLISHMENTS 

For  the  past  three  years,  under  the  guidance  of  my  mentor,  Dr.  Ralph  P.  Mason,  I 
have  gone  through  a  rigorous  training  in  tumor  biology,  tumor  histology,  tumor 
modeling,  tumor  implantation  and  transplantation,  radiation  biology,  MRI  physics,  MRI 
RF  coil  design,  construction,  and  testing,  MRI  data  acquisitions,  computer  programming, 
digital  signal  and  image  processing,  near-infrared  spectroscopy  (NIRS),  experimental 
design,  and  data  analysis.  I  have  learned  how  to  use  many  advanced  medical  systems  and 
instruments.  As  an  important  part  of  my  doctoral  curriculum  training  and  research,  I 
needed  to  use  an  Omega  CSI  4.7  T  MR  system  with  actively  shielded  gradients 
(Acustar7^,  Bruker  Instruments,  Inc.,  Fremont,  CA,  USA).  This  system,  based  on  a  40- 
cm  diameter  bore  horizontal  magnet,  is  located  in  the  Rogers  Magnetic  Resonance  Center 
of  the  University  of  Texas  Southwestern  Medical  Center  at  Dallas,  which  is  an  NIH 
Biotechnology  Resource  Facility.  I  had  immediate  access  to  the  system  and  other 
resources  in  the  center,  and,  on  average,  I  used  the  system  three  days  per  week  in  the  past 
for  either  experiments  or  data  processing  and  programming.  Now  I  can  operate  the 
magnet  independently  for  both  imaging  and  spectroscopy  experiments.  To  meet  our 
particular  experimental  needs,  I  modified  and  wrote  several  data  acquisition  and  post¬ 
processing  programs.  These  include  an  NMR-shell  script  program  for  computing  the 
biological  half-life  of  hexafluorobenzene  (HFB)  in  rat  breast  tumors,  an  NMR-shell 
script  program  for  displaying  the  T1  relaxation  curve  of  individual  voxels  of  MR  EPI 
images,  a  C  and  NMR-shell  script  program  for  acquiring  BOLD  and  Gd-DTPA  data,  a  C 
program  for  converting  MRI  image  files  among  different  platforms,  a  Visual  Basic 
program  for  processing  MRI  EPI  data,  and  a  LabView  program  for  processing  near- 
infrared  spectroscopy  (NIRS)  data.  In  addition,  I  have  also  learned  how  to  use  electronic 
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instruments  commonly  seen  in  a  modem  RF  lab,  including  a  sophisticated  HP  Frequency 
Analyzer,  and  how  to  operate  center  lathes  and  milling  machines  frequently  encountered 
in  the  machine  shop  of  a  radiology  department.  These  are  the  necessary  skills  required  for 
a  professional  medical  physicist  and  biomedical  engineer.  With  these  skills,  I  designed, 
constructed,  and  tested  a  double-tuned  (19F-'H)  birdcage  resonator  and  a  switchable 
slotted  tube  resonator.  Both  phantom  testing  and  in  vivo  imaging  experiments  indicated 
that  the  two  resonators  worked  well.  In  addition,  I  have  been  working  on  a  new  type  of 
RF  coil,  which  combines  a  traditional  RF  coil  with  NIRS.  These  coils  will  allow 
simultaneous  measurement  of  tumor  tissue  pC>2  by  19F  MRI,  tumor  vascular 
oxyhemoglobin  concentration  [Hb02],  and  total  hemoglobin  concentration  [Hb]T0tai  by 
NIRS.  This  is  a  challenging  task  and  successful  completion  of  these  RF  coils  would  be 
significant  in  tumor  oximetry.  Another  goal  of  this  Predoctoral  Traineeship  was  to 
investigate  breast  tumor  physiology  in  response  to  therapeutic  interventions.  This 
required  the  use  of  rat  mammary  adenocarcinoma  13762NF.  Implanting  this  type  of 
tumor  in  female  Fisher  344  rats  in  a  pedicle  model  was  a  complicated  and  lengthy 
surgical  procedure.  Now  I  have  learned  the  techniques  and  can  do  the  surgery  without 
any  difficulties.  I  have  also  learned  the  techniques  of  blood  gas  analysis  using  fiber  optic 
pulse  oximeter  and  the  operation  of  automated  microelectrode  systems.  Besides 
comprehensive  hands-on  training,  I  developed  a  mathematical  model  for  describing 
tumor  hemodynamics  following  physiological  perturbations,  and  a  second  mathematical 
model  for  computing  tumor  oxygen  consumption  following  KCl-induced  cardiac  arrest 
using  NIRS. 


3.  RESEARCH  ACCOMPLISHMENTS 

1)  Background  Information 

Solid  tumors  develop  regions  of  hypoxia  during  their  growth  due  to  an  imbalance 
between  the  rate  of  tumor  cell  proliferation  and  the  proliferation  and  branching  of  the 
blood  vessels  [1-3],  leading  to  diffusion-limited  or  chronic  hypoxia.  However,  tumor 
hypoxia  can  also  occur  due  to  another  important  reason,  the  deficient  oxygen-carrying 
capacity  of  the  blood,  i.e.,  low  blood  hemoglobin  concentration,  as  observed  in  the  case 
of  anemic  cancer  patients  [4-6],  Substantial  clinical  evidence  has  indicated  that  it  is  the 
hypoxia  that  is  responsible  for  the  failure  of  radiotherapy  [7-10],  some  forms  of 
chemotherapy  [11,  12],  and  photodynamic  therapy  [13],  Moreover,  considerable  clinical 
studies  have  also  shown  that  tumor  hypoxia  plays  a  key  role  for  increased  expression  of 
many  genes  that  relate  to  tumor  angiogenesis  and  growth,  including  vascular  endothelial 
growth  factor  (VEGF),  platelet-derived  growth  factor  (PDGF),  and  p53  [14,  15],  and  low 
oxygenation  levels  in  tumors  correlate  with  a  higher  metastatic  rate  [16-18],  Although  the 
mechanism  for  the  latter  effect  has  not  been  fully  understood,  some  studies  have 
suggested  that  it  could  be  caused  by  increased  DNA  mutation  rate  and  overreplication 
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under  the  hypoxic  condition  [19-21],  In  addition,  a  number  of  clinical  trials  have  found 
that  patient  survival,  measured  either  as  tumor  regression  or  as  local  control,  depends 
largely  on  tumor  oxygenation  [10,  16,  22,  23], 

Recently,  the  use  of  erythropoietin  treatment  and  blood  transfusion  has  gained 
increasing  attention  and  popularity  in  tumor  therapy  community  [24,  25]  because  many 
clinical  studies  have  shown  that  there  is  a  relationship  between  hemoglobin 
concentrations  and  therapeutic  outcome  [26,  27],  Results  of  several  other  clinical  trials 
have  indicated  that  higher  hemoglobin  concentrations  correlate  with  improved  local 
tumor  control  and  a  higher  overall  survival  rate  [28,  29],  suggesting  that  hemoglobin 
concentration  could  be  an  independent  prognostic  factor  in  tumor  therapy.  However,  the 
exact  mechanism  how  hemoglobin  concentrations  affect  tumor  therapy  and  radiotherapy 
in  particular  is  still  being  debated  [30],  Available  experimental  and  clinical  data  seem  to 
support  the  hypothesis  that  low  hemoglobin  concentrations  impair  oxygen-transporting 
capacity  of  blood  and  eventually  lead  to  tumor  hypoxia.  In  view  of  the  important  role  of 
hemoglobin  in  tumor  therapy,  it  is  necessary  to  have  a  means  to  measure  and  manipulate 
hemoglobin  levels  in  tumor. 

The  critical  effects  of  oxygenation  on  tumor  therapy,  angiogenesis,  metastasis,  and 
prognosis  have  stimulated  the  development  of  novel  tumor  oximetry  techniques.  Over  the 
past  several  decades,  substantial  progress  has  been  made  in  developing  techniques  for 
measuring  tumor  tissue  oxygen  tension  (p02 )  and  hypoxia.  These  include 
microelectrodes  [31,  32],  ESR/EPR  [33,  34],  the  comet  assay  [35],  phosphorescence 
quenching  imaging  [36],  nitroimidazole  binding  assays  [37],  paired  survival  assay  [37] 
and  NMR  [38-40],  However,  none  of  these  has  been  recognized  generally  as  a  perfect 
noninvasive  method.  Electrodes  are  highly  invasive  and  only  sample  a  limited  region. 
ESR  has  two  major  drawbacks  compared  with  NMR:  lack  of  tissue  penetration  at  very 
high  frequencies  (GHz)  and  the  inability  to  obtain  functional  information  at  multiple 
locations  (mapping).  31 P  NMR  is  perhaps  the  most  attractive  indicator  of  tumor 
oxygenation  since  it  is  entirely  non-invasive  with  observation  of  endogenous  metabolites. 
However,  it  is  relatively  insensitive,  sampling  a  large  volume,  and  metabolic  hypoxia 
may  occur  at  an  oxygen  tension  that  considerably  exceeds  radiobiological  hypoxia  [41]. 
Recent  studies  have  indicated  that  tissue  contrast  changes  in  *H  MRI  on  the  basis  of  the 
blood  oxygen  level  dependent  (BOLD)  method  provide  a  qualitative  approach  to  tissue 
oxygenation.  However,  changes  in  blood  flow  affect  the  signal  intensity  and  therefore, 
complicate  the  interpretation  of  results  [42,  43],  19F  NMR  techniques  can  provide  a  direct 
measurement  of  p02  based  on  the  principle  that  the  19F  NMR  spin-lattice  relaxation  rates 
R1  (=1/T1)  of  perfluorocarbon  (PFC)  emulsions  are  linearly  proportional  to  oxygen 
tension  [44-46],  19F  has  100%  natural  isotopic  abundance  and  an  83%  sensitivity  relative 
to  *H.  i9F  occurs  in  exceedingly  low  concentrations  in  biological  systems  (as  the  fluoride 
ion),  thus,  there  is  essentially  no  background  noise  to  interfere  with  in  vivo  studies.  I9F 
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NMR  techniques  have  the  potential  advantages  of  being  non-invasive,  repeatable,  and 
able  to  produce  local  tumor  p02  maps  with  high  spatial  resolution. 

Tumor  oxyhemoglobin  concentration  [Hb02]  or  .s02  is  another  important  indicator 
of  tumor  oxygenation.  However,  for  many  years,  very  little  progress  has  been  made  in 
developing  techniques  to  measure  it.  One  promising  technique  is  near-infrared 
spectroscopy  (NIRS).  The  NIRS  tumor  oximetry  is  based  on  the  fact  that  there  exists  a 
substantial  absorption  difference  of  light  in  NIR  region  (700  ~  900  nm)  between 
deoxyhemoglobin  [Hb]  and  oxyhemoglogin  [HbOs],  while  absorption  of  light  by  other 
macromolecules  and  water  is  insignificant.  Based  on  light  modulation  mechanisms,  NIRS 
techniques  can  be  classified  as:  continuous  wave  (CW)  light  spectroscopy  in  DC  format, 
amplitude-modulated  laser  light  spectroscopy  in  frequency  domain,  and  pulsed-laser  light 
spectroscopy  in  time  domain.  Using  NIRS  techniques,  in  vivo  measurements  of  5O2  have 
been  carried  out  in  a  wide  variety  of  biological  systems  such  as  exercised  muscles  [47- 
49]  and  brain  [50-52],  However,  so  far,  very  limited  studies  have  been  done  in  in  vivo 
measurements  of  tumor  ,v02.  Among  the  published  studies,  NIRS  has  been  used  to 
evaluate  the  effects  of  anesthetics  on  vascular  s02  in  RIF-1  tumors  [53]  and  in  9L 
gliosarcoma  [54],  and  to  assess  the  effects  of  hypoxia,  hyperoxia,  and  asphyxia  on 
BA1112  rhabdomyosarcoma  [55],  Recently,  Mariya  et.  al.  [56]  has  published  their  data 
on  monitoring  tumor  oxygenation  status  during  fractionated  irradiation  in  two  murine 
tumor  cell  lines  using  NIR  reflection  spectroscopy.  NIRS  is  completely  non-invasive, 
inexpensive,  portable,  and  amenable  to  real-time  measurements.  It  can  be  used  to 
measure  tumor  vascular  hemoglobin  oxygen  saturation  s02,  hemoglobin  concentration 
[Hb],  and  tumor  blood  flow,  and  to  monitor  tumor  transient  response  to  therapeutic 
interventions.  However,  because  of  inhomogeneous  nature  and  limited  dimension  of 
tumors,  classical  diffusion  theory  does  not  hold.  Absolute  quantification  of  tumor 
oxygenation  based  on  photon  diffusion  approximation  approach  still  remains  a 
challenging  issue.  An  alternative  approach  is  to  modify  Beer-Lambert's  law  and  use  the 
measured  transmitted  light  amplitude  to  compute  the  trends  in  the  changing  absorption 
coefficients  and,  thus,  changes  in  [Hb02]  or  s02  and  [Hb]Totai- 

We  have  surveyed  the  relative  sensitivity  of  several  PFCs  and  found  that 
hexafluorobenzene  (HFB)  offers  exceptional  sensitivity  to  changes  in  p02  with  relatively 
little  response  to  temperature.  HFB  has  a  single  resonance,  providing  optimal  signal-to- 
noise  ratio  (SNR).  I  have  applied  the  19F  MRI  of  HFB  technique  to  investigate  dynamic 
changes  in  p02  in  tumors  in  response  to  respiratory  challenges  and  the  feasibility  of 
mapping  the  clearance  rate  of  HFB.  In  addition,  I  have  also  investigated  hemodynamic 
changes  in  [Hb02]  and  [Hb]Totai  in  tumors  in  response  to  therapeutic  interventions  using  a 
newly  developed  NIRS  system. 

2)  Tumor  Transplantation  and  Handling,  Tumor  Volume  Doubling  Time 
(VDT),  and  Tumor  Histology 
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Tumor  Model 

Murine  mammary  adenocarcinomas  13762NF  [57]  were  used  in  this  study.  The 
13762NF  is  a  subline  of  mammary  adenocarcinoma  13762.  This  tumor  model  was  chosen 
because  it  demonstrates  substantially  different  therapeutic  sensitivity  and  metastatic 
characteristics  from  other  tumor  types.  The  13762NF  is  less  differentiated.  It  is  inhibited 
by  estrogen  in  young  female  rats  and  stimulated  in  adult  female  rats.  It  metastasizes  to 
regional  lymph  nodes,  and  lung,  and  occasionally  to  liver,  but  not  to  brain.  It  is  highly 
responsive  to  alkylating  agents  and  platinum  chemotherapeutic  agents. 

Tumor  Transplantation  and  Handling 

Our  lab  has  successfully  developed  a  tumor  transplantation  model,  the  pedicle 
model  [58],  Murine  mammary  adenocarcinomas  13762NF  were  implanted  in  skin 
pedicles  on  the  forebacks  of  adult  female  Fischer  344  rats  (-250  g).  To  identify  them, 
each  tumor  was  assigned  a  unique  code.  The  surgical  procedure  for  creating  a  pedicle 
tumor  model  is  as  follows.  A  flap  of  depilated  skin  was  raised  from  the  body  of  the  rat 
and  held  in  position  with  a  non-traumatic  curved  bull-dog  clip.  A  3 -cm  incision  was 
made  through  die  skin  using  the  curved  edge  of  the  clip  as  a  guide.  Wound  clips  were 
used  to  joint  the  edge  of  the  skin,  producing  a  tube  resembling  a  suitcase  handle.  Animals 
were  housed  separately  after  the  surgery  and  during  the  entire  course  of  experiments. 
Two  weeks  later,  the  clips  were  removed  and  the  distal  end  of  the  pedicle  severed.  A 
piece  of  fresh  tumor  tissue  (~  2x  2  x  2  mm3)  from  the  first  generation  of  mammary 
adenocarcinomas  13762NF  was  implanted  in  the  lumen  and  cut  closed  with  a  wound  clip. 

One  of  the  major  advantages  of  the  pedicle  model  is  that  it  is  basically  isolated 
from  the  body  proper  and  perfectly  suited  for  in  vivo  NIRS  and  MRI  studies,  therapy,  and 
manipulation.  In  addition,  the  pedicle  model  allows  accurate  measurement  of  tumor  size. 
It  has  been  shown  to  have  no  significant  difference  from  the  traditional  subcutaneous  site 
in  the  thigh  in  terms  of  growth  [58], 

Tumor  Growth  and  Mathematical  Model  for  Computing  VDT 

After  tumors  were  implanted,  they  first  entered  a  two-week  silent  period  of 
undetected  growth  and  then  an  accelerating  period  in  which  most  tumors  grew 
exponentially.  Some,  though  they  may  not  appear  to  grow  exponentially  over  their  entire 
life  spans,  did  show  exponential  growth  over  short  periods  of  time.  The  time  for  the 
tumors  to  grow  to  a  predetermined  size  spanned  a  wide  range.  In  some  cases,  the  silent 
periods  lasted  more  than  three  weeks.  Once  the  implanted  tumor  tissue  grew  into  a 
detectable  tumor,  tumor's  three  orthogonal  dimensions  were  measured  at  least  once  every 
two  days  with  a  caliper.  Depending  on  tumor  growth  rate,  tumors  were  measured  daily  in 
some  cases  to  improve  the  goodness  of  curve  fit  when  computing  volume  doubling  time 
(VDT).  Prior  to  the  measurements,  the  rats  were  anesthetized  with  50 pi  ketamine 
hydrochloride  (100  mg/ml)  and  the  tumors’  hair  was  cut  with  a  pair  of  surgical  scissors  to 
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improve  the  accuracy  of  the  measurements.  Using  an  ellipsoidal  approximation,  tumor 
volume  was  determined  using  formula 


T/  An.  ,a  b  c  n  , 

K=(TH2H2H2)“(6 yC"b'C 


(1) 


where  a,  b,  and  c  were  the  diameters  along  the  three  major  orthogonal  axes  of  the  tumor. 
A  very  important  proliferative  feature  of  the  solid  tumors  is  their  volume  doubling  time 
( VDT)  (day),  a  kinetic  parameter  closely  related  to  the  underlying  mathematical  growth 
model. 


Historically,  a  growing  tumor  is  modeled  as  a  deterministic  dynamic  system 
mathematically  described  by  ordinary  differential  equations.  The  central  concept  of  a 
dynamic  system  is  the  trajectory.  In  the  case  of  tumor  growth,  the  trajectory  is  the  growth 
curve  that  describes  the  change  in  tumor  size  as  a  function  of  time  from  the  start  of 
proliferation  of  initial  tumor  cells.  The  tumor  size  can  be  expressed  as  dimensions  or 
volume  or  mass  or  cellularity,  depending  on  personal  choice.  For  tumor  modeling,  these 
quantities  are  interchangeably  used,  because  they  are  linearly  proportional  to  each  other 
for  solid  tumors.  Unfortunately,  most  of  these  tumor  growth  models  have  been  seldom 
validated  against  experimental  tumor  growth  curves,  either  because  of  the  relative 
scarcity  of  high  quality  tumor  growth  data  [59],  Limited  accurately  measured  tumor  data 
revealed  that  growth  curves  for  some  human  tumors  are  very  close  to  exponential 
function  [60],  yielding  a  constant  volume  doubling  time.  However,  other  studies 
indicated  that  many  human  tumors  also  show  irregular  or  decelerating  growth,  giving  a 
progressively  longer  volume  doubling  time  [60],  Such  curves  can  be  fitted  with  a  variety 
of  mathematical  models,  but  the  most  prominent  one  is  the  Gompertz  model  [61], 
However,  the  empirical  Gompertz  model  lacks  a  truly  fundamental  biological 
explanation. 

In  an  attempt  to  obtain  an  analytical  expression  for  computing  the  tumor  volume 
doubling  time,  I  used  a  first  order,  autonomous  differential  equation  to  model  the  tumor 
growth 


dV 

—  =  XV,  A>0  (2) 

dt 

where  X  is  the  tumor  growth  rate  constant.  The  model  is  based  on  the  fundamental 
biological  argument  that  tumor  growth  results  from  exponential  cell  proliferation.  If,  for 
example,  each  healthy  tumor  cell  divides  to  produce  two  proliferative  daughters,  then  the 
growth  of  the  tumor  cell  population  is:  1,  2,  4,  8,  ....  2n,  an  exponential  function,  i.e.,  the 
population  size  will  increase  exponentially  with  time.  This  model  describes  unrestricted 
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tumor  growth  as  the  time  goes  to  infinity,  a  phenomenon  not  supported  by  experimental 
data.  However,  for  a  reasonable  period  of  growth,  it  gives  excellent  description  of  tumor 
growth  curves.  Most  importantly,  we  can  derive  an  analytical  expression  for  computing 
the  tumor  volume  doubling  time  from  this  model. 

The  solution  of  Equation  (2)  is  an  exponential  function:  V  =  V0  ■  exp(At)  or 
V  =  V0  •  exp(7 /  r) ,  where  V0(cm 3)  is  the  initial  tumor  volume  at  time  t  -  0  (day)  and 

t  (day)  is  the  time  constant  of  the  exponential  growth  curve.  A  very  important  property 
of  exponential  tumor  growth  is  that  the  tumor  volume  doubling  time  ( VDT)  is  a  constant. 
Thus,  when  time  increases  by  one  VDT,  the  tumor  volume  will  be  doubled. 
Mathematically,  this  is  written  as 


2V  =  V0  ■  exp( 


H-FDT 

T 


(3) 


Dividing  Equation  (3)  by  V  =  V()  -exp(t/r)  gives 


2V 

V 


T/  rt  +  VDT N 
K  '  exp( - ) 


V0  -exp(-) 

T 


2  =  exp( 


VDT^ 

T 


(4) 


Taking  natural  log  on  both  sides  of  Equation  (4)  gives  an  equation  for  computing  VDT 

VDT  =  rln2  (5) 


where  r  is  obtained  by  fitting  the  experimental  tumor  data  to  V  =  V0  ■  exp (t/  r) . 


Figure  1  shows  the  growth  curve  of  a  representative  mammary  adenocarcinoma 
13762NF.  This  particular  tumor  had  a  “silent  interval”  of  about  10  days.  When  its  growth 
was  detected,  the  tumor  volume  increased  exponentially  with  time  ( R  =  0.954)  over  a 
fairly  long  period  of  time  as  seen  in  Figure  1,  indicating  that  the  assumption  for  the 
exponential  growth  model  was  correct.  The  volume  doubling  time  (VDT)  as  determined 
by  Equation  (5)  was  approximately  5.2  days.  A  few  points  worth  mentioning  here.  First 
of  all,  the  long  silent  phase  of  undetected  growth  was  not  included  in  the  plot.  If  it  were 
to  be  included,  the  entire  curve  would  shift  10  days  to  the  right  and  the  shape  of  the 
exponential  part  of  the  growth  curve  would  not  changed  at  all.  This  would  not,  in  any 
way,  affect  the  result  of  VDT  computation  as  the  volume  doubling  time  is  a  kinetic 
parameter  used  to  describe  the  accelerating  phase  or,  in  most  cases,  the  exponential  phase 
of  a  tumor  growth.  Secondly,  in  measuring  the  tumor  dimensions,  especially  when 
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tumors  were  small,  the  skin  thickness  was  subtracted  from  the  measured  dimensions. 
Last,  the  VDTs  were  computed  solely  based  on  the  data  obtained  from  intact  tumors. 
During  their  late  growth  stage,  some  tumors  might  bleed  due  to  a  variety  of  reasons. 
When  this  happened,  the  data  gathered  after  this  point  were  not  included  in  the  final 
computation  of  VDTs. 


Figure  1.  Tumor  volume  of  a  representative  13762NF  breast  tumor  as  a 
function  of  time  after  implantation.  The  tumor  volume  was  computed  using 
an  ellipsoidal  model:  V  =  (7i/6)abc.  An  excellent  curve  fit  (R  =  0.954)  was 
obtained  using  an  exponential  growth  model. 


Table  1.  A  summary  of  tumor  volume  doubling 
time  for  a  group  of  six  13762NF  breast  tumors. 


Tumor  No. 

VDT  (Day) 

1 

3.2 

2 

4.0 

3 

3.3 

4 

3.6 

5 

5.2 

6 

4.4 

Mean  VDT  (Day) 

3.95  ±0.76 

Tumor  Histology 

Conventional  tumor  histology  was  performed  to  gather  information  regarding 
tumor  size,  gross  morphology,  presence  and  extension  of  tumor  necrosis,  histological 
type  and  grade.  Randomly  chosen  tumors  were  sacrificed  by  tail  injection  of  a  lethal  dose 
of  KC1.  Immediately  following  the  death  of  the  rats,  the  tumors  were  excised  and  fixed  in 
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10%  neutral  buffered  formalin.  The  tumors  were  subsequently  processed  through  graded 
ethanols  and  xylene,  and  embedded  in  paraffin.  After  embedding,  sections  were  cut  at 
5  pm  thickness  and  routine  hematoxylin  and  eosin  (H&E)  staining  performed  according  to 
established  protocols.  Selected  central  and  peripheral  regions  of  the  stained  histological 
sections  were  photomicrographed  using  brightfield  optics  with  a  magnification  of  ten. 
Figure  2  shows  histological  sections  of  a  representative  13762NF  breast  tumor. 


Figure  2.  Histological  sections  of  a  representative  13762NF  breast 
tumor:  central  (top)  and  peripheral  (bottom)  regions  of  the  tumor. 


3)  Design  and  Construction  of  MRI RF  Resonators 

This  project  involved  19F  and  ‘H  MR  imaging  with  living  rats  in  a  single 
experiment,  so  it  was  highly  desirable  to  be  able  to  change  the  resonant  frequency 
without  retuning  the  resonator  and  disturbing  the  animal.  The  desired  imaging  and 
spectroscopy  resonator  for  our  purposes  would  be  the  one  that  provides  the  highest 
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possible  signal-to-noise  ratio  (SNR),  a  good  filling  factor,  a  high  Q  factor,  minimum 
resistance  losses,  and  a  highly  homogeneous  field.  To  address  these  issues,  I  designed, 
constructed,  and  tested  two  RF  resonators,  one  being  a  double-tuned  birdcage  resonator 
(’H/19F)  and  another  one  a  swithable  slotted  tube  resonator. 

Double-Tuned  ^H-^F)  Birdcage  Resonator 

A  birdcage  resonator  consists  of  a  set  of  N  (2,  4,  8,  16,  etc.)  copper  wires  or  sheets 
arranged  axially  on  the  surface  of  a  plastic  cylinder  and  connected  by  high  quality  RF 
capacitors  at  each  end.  In  this  configuration,  the  effective  current  density  in  the  wires  or 
sheets  varies  in  proportion  to  cos($),  where  <f)  is  the  azimuthal  angle  in  cylindrical 
coordinates  [62],  and  there  exist  several  possible  modes  of  resonance,  two  of  which  have 
the  desired  sinusoidal  dependence  of  current  on  $  [63], 

•  Design  Theory 

Assuming  perfectly  conducting  copper  sheets,  Figure  3  shows  part  of  the 
simplified  circuit  model  of  a  birdcage  resonator  with  16  identical  legs  (copper  sheets), 
where  L,  is  the  self  inductance  of  the  /th  copper  sheet  and  C,  is  the  capacitor  connected 
between  the  /th  and  the  (H-l)th  copper  sheet.  The  self  inductance  of  the  short  copper  sheet 
used  to  connect  the  capacitors  at  the  ends  is  neglected.  This  is  based  on  the  fact  that  the 
electric  current  in  the  resonator  is  longitudinal  (along  the  copper  sheet)  and  maximal  at 
the  ends  of  the  resonator. 


Cj-2  Cm  Cj  Cj+i  Cj+2 


Ci-2  Cm  Cj  Cj+i  Cj+2 


Figure  3.  Part  of  the  equivalent  circuit  for  the  high-pass 
double-tuned  whole  body  birdcage  resonator. 


By  applying  Kirchhoffs  Law  to  the  /th  loop  made  up  of  /th  and  (/+l)th  copper 
sheets  and  /th  capacitor  C,  and  using  the  complex  number  notation  for  AC  current  and 
impedance,  we  can  get  a  linear  system  of  equations: 
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(6) 


-I -—J,  =0 

k  k  i 

where  Ml  k  is  the  mutual  inductance  between  the  /th  and  &th  copper  sheets,  «/,  is  the 
current  in  the  /th  capacitor,  Ik  is  the  current  in  the  kx\\  copper  sheet,  /  is  the  imaginary 

number,  and  the  summation  is  carried  out  over  all  possible  values  of  index  k  from  0  to  TV 
-  1. 


There  are  many  ways  of  solving  Equation  (6),  either  analytically  or  numerically. 
The  solution  will  represent  the  various  resonant  modes  of  operation  of  the  resonator.  This 
circuit  model  is  a  simplified  one  and  does  not  take  into  account  the  complications  that 
exist  in  any  real  implementation,  such  as  the  self-inductance  of  the  capacitors  and  the 
inductance  of  the  end  rings.  Nevertheless,  the  general  form  of  the  model  is  not  altered  by 
these  factors,  and  the  qualitative  description  is  still  valid. 

•  Computer  Simulation 

Instead  of  solving  the  above-described  complex  equation,  I  simulated  this  circuit 
model  on  a  PC  computer  using  P-Spice  [64],  The  desired  resonant  frequencies  are  188.22 
MHz  for  19F  and  200. 16  MHz  for  ’H  in  a  4.7  T  magnet.  The  inductance  produced  by  the 
end  rings  was  neglected  in  this  simulation  circuit,  since  it  only  accounts  for  ~  5%  of  the 
total  inductance.  However,  the  circuit  model  is  adequate  for  analyzing  the  essential 
frequency  characteristics  of  the  resonator.  The  first  step  in  simulating  this  problem  was  to 
estimate  the  inductance  of  the  copper  sheets,  which  is  given  by  [63], 

L  =  2h[ln  (: 2h/b )  +  1/2  ]  (7) 

where  L  is  the  inductance  in  nanoHenries,  h  is  the  height  in  centimeters,  and  b  is  the 
width  in  centimeters.  The  copper  sheets  have  a  height  of  16.5  cm  and  a  width  of  1.0  cm, 
with  a  thickness  of  ~  0.2  mm.  The  calculated  value  of  the  inductance  based  on  Equation 
(7)  was  113  nH. 

The  second  step  was  to  substitute  this  value  into  the  simulation  circuit,  apply  an 
AC  current  source  to  the  circuit,  vary  the  value  of  each  of  the  32  capacitors,  and  observe 
the  amplitudes  of  output  voltages  and  resonant  frequencies  of  the  circuit  until  two  sharp 
resonant  peaks  (188.22  MHz  for  19F  and  200.16  MHz  for  ]H  in  a  4.7  T  magnet)  with  the 
maximum  output  voltages  were  achieved.  The  capacitance  values  for  the  end  ring 
capacitors  obtained  from  the  P-Spice  simulation  are  tabulated  in  Table  1.  Theoretically 
speaking,  all  the  capacitance  values  should  be  equal  for  a  single-tuned  resonator  and 
close  to  each  other  for  a  double-tuned  resonator  in  order  to  generate  a  homogeneous  Bx 
field.  The  circuit  model  presented  here  is  a  simplified  one.  Therefore,  the  capacitance 
values  obtained  are  not  quite  close  to  each  other.  The  simulation  study  indicates  that  CIO, 
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C26,  Cll,  C27,  Cl 2,  and  C28  are  relatively  insensitive  to  l9F  resonant  frequency,  thus, 
they  could  be  connected  to  the  'H  tuning  circuit.  I  chose  to  place  the  tuning  circuit  in 
parallel  to  CIO  since  this  placement  provides  both  the  best  tuning  for  'H  mode  and  also 
stabilizes  the  19F  mode. 


Table  1  Capacitance  values  determined  by  P-Spice  simulation  (unit:  pF) 


Cl 

C2 

C3 

C4 

C5 

C6 

Cl 

C8 

pi Igl 

Jjifc 

130 

120- 

pTl5  . 

108 

115 

C9 

CIO 

Cll 

C12 

C13 

C14 

C15 

C16 

13 pg! 

ipS! 

185 

120 

Ip§t02 

90 

90 

C17 

C18 

C19 

C20 

C21 

Cll 

C23 

C24 

1  40 

Si’..  !•. 

130 

130 

Wm  4 

§415 

108: 

106 

C25 

C26 

Cll 

C28 

C29 

C30 

C31 

C32 

1  42'  : 

IMSl 

HUB 

120 

120 

S§02  . 

90-il 

'  ^0  •  : 

Fine  tuning  of  I9F  mode  is  achieved  by  connecting  another  capacitor  CF19T  in 
parallel  with  C6.  When  19F  tuning  capacitor  CF19T  varies  from  0  pF  to  2000  pF,  19F 
resonant  peak  varies  from  (188.224  MHZ,  14.727  V)  to  (184.871  MHZ,  16.606  V)  while 
'H  resonant  peak  changes  only  from  (200.435  MHZ,  13.905  V)  to  (200.535  MHz,  13.602 
V),  which  is  a  clear  indication  that  each  mode  is  independently  tuned. 

•  Construction 

A  16-leg  high-pass  double-tuned  birdcage  resonator  was  constructed  according  to 
the  design  principles  and  the  simulation  results  described  above.  The  resonator  is  tunable 
to  both  188.2  MHZ  and  200.1  MHZ  for  operation  at  19F  and  ’H  frequencies  on  a  4.7  T 
magnet.  The  cylindrical  resonator  body  is  24.60  cm  in  length  and  has  an  outer  diameter 
of  24.00  cm.  The  outer  structure  of  the  resonator  is  made  of  a  5  mm  thick  plexiglass 
cylinder,  which  provides  mechanical  support  and  structural  stability  for  the  resonator.  A 
0.2  mm  thick  copper  sheet  shield,  connected  to  the  ground,  is  mounted  on  the  inner  wall 
of  the  outer  plexiglass  cylinder. 

The  critical  part  of  the  resonator  is  its  inner  structure.  The  inner  structure,  made 
from  a  5  mm  thick  plexiglass  cylinder,  has  an  outer  diameter  of  10.50  cm.  On  the  inner 
wall  of  the  plexiglass  cylinder  are  mounted  directly  16  copper  sheets  (legs)  with  1.5  mm 
gaps  and  32  capacitors.  The  size  of  the  uniform  region  of  B ,  magnetic  field  generated  by 
the  resonator  is  directly  related  to  the  number  of  legs.  The  more  numerous  the  legs,  the 
more  uniform  the  Bx  magnetic  field  in  the  radial  direction.  I  chose  to  build  a  birdcage 
resonator  with  sixteen  legs  because  the  Bx  magnetic  field  generated  should,  in  principle, 
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be  uniform  enough  for  my  applications  and  should  produce  strong  NMR  signal.  The  legs 
are  made  from  0. 1  mm  thick  copper  foil.  Each  of  the  sixteen  legs  is  16.5  cm  in  length,  1.0 
cm  in  width  in  the  central  part  and  1.5  cm  in  width  at  the  two  ends  (Figure  4). 


I  i 

^  2.0  cm  f*" 

I  I 


1.5  cih 


Figure  4.  Geometry  of  the  legs  (copper  sheets). 


The  structure  of  the  two  end  rings,  i.e.,  the  two  circular  bands  at  the  ends  of  the 
resonator,  consisted  of  capacitors  and  was  made  identical  to  achieve  electrical  balance  so 
that  the  loading  effects  of  the  sample  could  be  reduced.  The  resonator  is  double-tuned, 
and  is  driven  between  equivalent  points  on  the  two  rings  halfway  between  two  adjacent 
parallel  legs.  The  resonator  was  tuned  by  individually  adjusting  the  CF19T  and  CHIT  at 
each  of  the  resonant  frequencies.  The  capacitively  coupled  tuning  loops  were  placed 
between  legs  11  and  13  and  legs  19  and  21.  These  loops  were  connected  by  standard 
coaxial  cables  to  the  external  tuning  capacitors  CF19T  and  CHIT  that  allowed  each 
resonant  mode  to  be  independently  tuned.  Special  care  was  taken  to  ensure  that  there  was 
little  tuning  interaction  between  the  two  resonant  frequencies  as  possible.  The  resonator 
was  matched  capacitively  through  a  variable  capacitor  Cmatch 

•  Bench  Testing 

The  laboratory  bench  testing  was  carried  out  using  a  frequency  analyzer 
(HP8752A).  To  test  the  resonator’s  performance,  it  was  essential  to  have  a  set  of  standard 
samples  that  simulated  the  type  of  in  vivo  samples  used  in  the  magnet.  Two  NMR 
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phantoms  were  prepared  for  this  purpose,  which  roughly  had  the  same  loading  effects 
(lowering  Q  and  shifting  resonant  frequencies)  as  a  male  adult  rat  (~  250  g).  A  444  ml 
phantom  was  made  of  an  NMR  compatible  plastic  bottle  filled  with  a  solution  of  saline 
and  trifluoroacetic  acid  (TFA,  5%  v/v).  A  smaller  250  ml  phantom,  also  made  of  an  NMR 
compatible  plastic  bottle,  was  filled  with  a  solution  of  TFA  (5%  V/V)  and  CuS04  (1 
g/liter)  in  distilled  water.  The  major  chemical  components  and  their  concentrations  of  the 
phantoms  are  listed  in  Table  2. 


Table  2.  Major  chemical  components  and  their 
concentrations  of  the  saline  phantoms. 


niumsun? 

5SSBB 

NaCl 

6.90 

KCl 

0.375 

MgS04 

0.145 

CaCh 

0.185 

Na2HC03 

2.10 

Glucose 

1.80 

The  resonator’s  resonant  frequencies,  unloaded  and  loaded  Q  values  (quality 
factors)  and  Q  damping  for  each  of  the  resonant  modes  as  measured  by  the  frequency 
analyzer  are  presented  in  Table  3.  go  is  the  unloaded  Q  value  measured  with  the  resonator 
in  air  (without  the  phantom),  while  Qh  is  the  Q  value  measured  when  the  saline  phantom 
was  placed  centrally  inside  the  resonator.  The  Q  damping  is  defined  as  the  ratio  of  Qo/Ql  ■ 
The  results  are  the  means  of  three  measurements.  The  loaded  Q  values  are  affected  by  the 
dimensions  of  the  phantom  and  its  position  in  the  resonator.  In  general,  the  loaded  Q 
values  are  dominated  by  the  effect  of  the  RF  lossy  sample.  The  Q  values  decrease  as  the 
dimensions  of  the  sample  increases. 


Table  3.  Results  of  the  laboratory  bench  testing  of  the  birdcage 
resonator.  Qn  and  Q\J  were  measured  with  the  444  ml  phantom. 


Resonant 

Mode 

Resonant 

Frequency 

(MHz) 

Unloaded 

Qo 

Loaded 

Ql 

Qo/Ql 

19f 

188.22 

14 

9 

1.56 

’H 

200.16 

13 

8 

1.63 
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It  can  be  seen  from  Table  3  that  the  unloaded  Q0,  loaded  QL,  and  the  Q  damping  Q0 
/<2Lfor  both  resonant  modes  are  similar.  The  relatively  low  Q  values  may  be  attributed  to 
the  relatively  long  coaxial  cables  used  to  connect  the  tuning  and  match  circuits.  Another 
possible  reason  could  be  the  large  size  of  the  resonator  because  the  longer  the  legs,  the 
larger  their  resistance,  causing  increased  resistance  losses  and  significantly  lowering  Q 
values.  One  of  the  consequences  of  the  low  Q  values  on  the  NMR  experiments  would  be 
to  increase  the  90°  pulse  width  for  a  given  RF  power  level.  This  can  be  verified  by  the 
following  equation  [64]: 


( 

PW=K 


Q) 


(8) 


where  PW  is  the  90°  pulse  width,  Vr  is  the  resonator  volume,  Q  is  the  quality  factor  of 
the  resonator,  and  the  factor  K  is  roughly  a  constant  for  all  resonator  designs.  This 
equation  indicates  that  the  90°  pulse  width  is  linearly  proportional  to  the  square  root  of 
resonator’  volume  and  inversely  proportional  to  the  square  root  of  resonator’s  Q  value. 

•  Phantom  Imaging 

The  double-tuned  birdcage  resonator’s  NMR  performance  was  tested  on  the 
Omega  CSI  4.7  T  magnet,  which  operates  at  a  proton  frequency  of  200.106  MHZ  and  a 
fluorine  frequency  of  188.273  MHZ.  The  phantom  was  centered  in  the  resonator  and  the 
resonator  was  positioned  in  the  isocenter  of  the  magnet,  with  its  axis  being  aligned  with 
the  direction  of  the  magnet’s  static  magnetic  field  B0.  The  90°  pulse  widths  were 
determined  from  a  180°  null  of  the  whole  sample.  Typical  90°  pulse  widths  were  700  ps 
for  proton  and  350  ps  for  fluorine.  Shimming  was  accomplished  on  the  water  proton  FID 
of  the  saline  solution  in  the  phantom,  with  the  resonator  tuned  to  the  double  resonant 
modes.  The  goal  was  to  obtain  a  B0  of  homogeneous  strength  over  the  entire  imaging 
sample  volume.  Fourteen  shim  currents  were  adjusted  on  the  proton  FID  throught  the 
system  automated  shimming  utility  until  a  symmetrical  and  long  FID  was  obtained.  The 
proton  spectral  line  width  at  half  height  after  shimming  was  34  Hz,  which  was  sufficient 
for  imaging  experiments. 

To  minimize  acquisition  time,  a  driven  equilibrium  spin-echo  pulse  sequence  with 
a  fairly  short  echo  time  (TE  =  8ms)  was  employed  so  that  the  T2  dephasing  effects  were 
not  significant.  Proton  transaxial  (transverse)  images  were  acquired  as  a  3-D  data  set  with 
a  repetition  time  TR  =  150  ms,  echo  time  TE  =  8  ms,  and  a  128  x  64  x  8  matrix  size  over 
100  mm  x  100  mm  field  of  view  (FOV),  providing  0.78  mm  x  1.56  mm  x  25  mm  digital 
resolution.  To  improve  signal-to-noise  ratio  (SNR),  two  acquisitions  were  averaged  for 
each  image,  thus,  giving  a  total  acquisition  time  of  2  min  33  sec.  Proton  coronal  and 
sagittal  images  were  also  acquired  with  the  same  imaging  parameters. 
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SNRs  for  the  images  with  2  averages  were  measured  using  a  standard  routine.  The 
values  for  signal  were  measured  from  a  randomly  chosen  portion  of  images  by  a  square 
crop  and  the  values  for  noise  were  taken  from  the  four  comers  of  the  image  background. 
For  each  measurement,  the  square  crop  covered  a  reasonably  large  area  of  the  image  or 
background,  so  that  the  data  obtained  were  unbiased.  The  results  of  SNR  measurements 
for  the  proton  images  are  listed  in  Table  4.  The  resonator’s  Bi  field  homogeneity  was 
estimated  by  visual  inspection  of  the  images. 


Table  4.  Results  of  SNR  measurements  for  the  proton  images. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower- 

Left 

SNR 

Lower 

-Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

88.0 

97.4 

93.0 

93.3 

92.9 

3.9 

Sagittal  Section 

279.8 

272.7 

278.7 

294.0 

281.3 

9.0 

Coronal  Section 

287.5 

114.3 

267.7 

125.8 

198.8 

91.4 

The  magnet  frequency  was  then  set  on  resonance  for  the  CF3  group  of  TFA, 
without  retuning  the  resonator.  The  90°  pulse  was  360  ms.  The  ,9F  spin-echo  images  were 
acquired  as  transaxial,  saggital  and  coronal  projections.  The  images  were  fluorine 
density- weighted  with  TR  =  150  ms  and  TE  =  8  ms  and  were  averaged  16  times  to 
achieve  an  acceptable  SNR.  Thus,  each  image  was  acquired  in  2  min  33  sec.  The  image 
matrix  size  was  128  x  64  with  a  field  of  view  of  100  x  100  mm  for  the  tranaxial 
projection  and  200  x  200  mm  for  the  sagittal  and  coronal  projections,  respectively. 


Table  5.  Results  of  SNR  measurements  for  the  19F  images. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower 

-Left 

SNR 

Lower 

-Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

3.2 

3.3 

3.4 

3.5 

3.3 

0.1 

Sagittal  Section 

5.9 

6.7 

7.3 

6.9 

6.7 

0.6 

Coronal  Section 

4.3 

4.5 

3.8 

3.7 

4.1 

0.4 

The  image  SNRs  were  measured  in  the  same  way  as  the  proton  images  and  the 
results  are  presented  in  Table  5.  The  spin-lattice  relaxation  time  Tl  of  the  CF3  group  was 
estimated  by  applying  a  non-spatially  selective  inversion  recovery  (IR)  RF  pulse.  Seven 
different  IR  delays  (x)  increasing  in  the  range  between  800  ms  to  20  sec  were  used  (Table 
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6).  The  Tl  value  was  calculated  based  on  the  Levenberg-Marquardt  three-parameter 
fitting  algorithm  [65]  on  peak  intensity  values  (Table  6). 


Table  6.  Results  of  Tl  measurement  for  TFA  solution  and 
delay  list  ( r )  of  the  inversion  recovery  pulse  sequence. 

Function:  y  =  A*(l-(W+l)*exp(-x/Tl) 

A  =  1.92e+06  ±26484 
W  =  0.026  ±0.042 
Tl  =2.3  ±0.19 
Standard  Error  =  41129.6 
Y  Standard  Deviation  =  31091.1 


T 

Intensity 

Calculated 

Difference 

SD 

20.000 

1923050.524 

1923610.434 

559.911 

0.018 

16.000 

1932423.278 

1921888.260 

-10535.018 

-0.339 

8.000 

1857983.273 

1859578.430 

1595.157 

0.051 

4.000 

1534324.010 

1567419.792 

33095.782 

1.064 

2.000 

1135584.863 

1085041.321 

-50543.542 

-1.626 

1.000 

587685.972 

637132.244 

49446.272 

1.590 

0.800 

545796.816 

522178.645 

-23618.170 

-0.760 

Slotted  Tube  Resonator 

The  slotted-tube  resonator  was  originally  developed  to  efficiently  provide  'H 
decoupling  to  samples  in  high-field  magnets  [66],  The  common  slotted  tube  resonator 
consists  of  a  conducting  tube  with  two  slots  cut  symmetrically  along  both  sides  of  the 
tube.  The  surface  current  distribution  produced  by  the  slotted  tube  resonator  peaks  near 
the  edges  of  the  conducting  tube  adjacent  to  the  slots.  The  conducting  tube  of  the  slotted 
tube  resonator  is  opaque  to  RF  field  and  diverts  all  the  Bx  field  flux  through  the  two  slots. 
Thus,  the  aperture  angle  of  the  slots  must  be  optimized  in  order  for  the  resonator  to 
generate  a  homogeneous  field.  The  fundamental  feature  of  the  slotted  tube  resonator  is 
that  the  conducting  component  has  a  low  inductance,  and  hence,  large  volume  resonators 
can  be  constructed.  At  the  RF  frequencies  encountered  in  clinical  MRI  systems,  this 
means  that  the  slotted  tube  resonators  can  be  designed  to  contain  samples  of  size 
comparable  to  that  of  the  human  head.  The  major  applications  of  the  slotted  tube 
resonator  are  for  imaging  and  for  localized  spectroscopy  techniques  demanding  a  uniform 
flip  angle  over  the  entire  sample. 

•  Design  Theory 


21 


The  circuit  model  is  illustrated  in  Figure  5.  Ls  and  Rs  are  sample  inductance  and 
resistance,  respectively.  This  circuit  model  can  be  further  simplified  in  terms  of 
equivalent  circuit  (Figure  6)  with 


z,  i(L,CW  -1) 

2  co(C3  -C2(L3C3co2  -1)) 


(9) 


where  co  is  the  resonant  frequency  of  the  circuit  and  i  is  the  imaginary  number.  Thus,  the 
tuning  impedance  ZT  is  given  by  the  following  expression: 


Zj  +  Z2  -  +  Z2 

CX(D 

.c,q,,c>2  -l)-o(C,  -c2(z,cy  -1)) 

fflC,(C,-C,(L,C,C2-l» 


(10) 


Rs 


Ls 

rrm 


c2 

fYYY\ 


Figure  5.  The  circuit  model  of  the  switchable  slotted  tube  resonator. 

Rs  is  the  sample  resistance  and  Ls  the  sample  inductance.  R3  is  the 
equivalent  resistance  of  the  tuning  loop. 

When  the  resonator  operates  at  the  higher  frequency  co  =  coh,  the  tuning  impedance 
ZT  is  a  function  of  the  total  capacitance  of  the  tuning  circuit  and  is  given  by  the 
expression 


where 


ZT  =-/ 


;(0+O 

cohCxC\ 


(11) 
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(12) 


^-.1  _  C2(L36>,  1)  C3 

L3(>,2-1 

Thus,  the  tuning  is  achieved  by  adjusting  all  the  capacitors.  When  the  resonator 
operates  at  the  lower  frequency  co=  coi,,  the  series  L3 C3  circuit  has  a  zero  impedance  and 
the  total  tuning  impedance  is  reduced  to 


ZT  = 


i 

CjtUi 


(13) 


Thus,  tuning  is  obtained  by  adjusting  C,  only. 

Zs 


Figure  6.  Equivalent  circuit  model  of  the  slotted  tube  resonator. 


•  Construction 

A  switchable  slotted  tube  resonator  was  constructed  based  on  the  above  design 
theory  using  two  copper  sheets,  which  cover  an  80°  arc  [67]  for  optimal  Bx  field 
homogeneity.  The  cylindrical  resonator  body  is  23.5  cm  long  with  outer  diameter  25.4 
cm.  The  outer  structure  of  the  resonator  is  made  of  a  5mm  thick  plexiglass  cylinder.  The 
inner  structure  is  made  of  a  3  mm  thick  plexiglass  cylinder,  having  a  small  outer  diameter 
of  7.6  cm  to  increase  the  filling  factor  and  reduce  the  pulse  width.  The  two  ends  of  the 
resonator  are  covered  by  two  6  mm  thick  plexiglass  boards  stabilized  by  6  nonmagnetic 
screws  on  each  end,  respectively.  The  two  boards,  along  with  the  outer  structure,  provide 
very  stable  mechanical  support  for  the  resonator.  Three  capacitors  of  20.0  pF  were 
positioned  symmetrically  between  the  guard  rings  and  the  copper  sheets  to  obtain 
uniform  surface  current  distribution. 

The  resonator  has  a  usable  length  of  15.5  cm,  sufficient  for  in  vivo  adult  rat  NMR 
experiments.  The  two  guard  rings  and  two  copper  sheets  are  made  from  0.3  mm  thick 
copper  foil.  The  width  of  each  of  the  guard  rings  is  1.0  cm  and  there  is  a  spacing  of 
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approximately  0.5  cm  between  the  guard  ring  and  the  vertical  copper  sheets,  which  are 
joined  with  capacitors.  The  copper  sheets  have  a  length  of  12.5  cm  and  their  width  covers 
about  80°.  The  resonator  is  tunable  to  l9F  (188.22  MHz)  and  'H  (200.16  MHz),  and  is 
driven  between  one  copper  sheet  and  one  guard  ring.  The  tuning  of  the  resonator  was 
achieved  by  three  variable  capacitors  Ci,  C2,  and  C3,  which  are  V2145H  family  high  Q 
RF  capacitors  (Voltronics  Corporation,  Denville,  New  Jersey).  Ci  and  C2  (NMTM38GE) 
have  a  capacitance  range  of  1.0  ~  38.0  pF,  C3  (NMTM120CE)  has  a  capacitance  range  of 
2.0  ~  120.0  pF,  and  CM  (NMQM22GE)  has  a  capacitance  range  of  1.0  ~  22.0  pF.  19F 
resonance  was  obtained  by  adjusting  C,  and  'H  resonance  was  achieved  by  adjusting  Q, 
C2,  and  C3,  respectively.  The  resonator  was  matched  capacitively  through  Cm- 

•  Bench  Testing 

The  resonator  was  tuned  and  matched  on  the  laboratory  bench  with  the  help  of  a 
frequency  analyzer  (HP8752A).  The  resonant  frequencies,  unloaded  Q0  and  loaded  QL 
values  and  Q  damping  for  each  of  the  resonant  modes  are  presented  in  Table  7.  Q0  is  the 
unloaded  Q  value  measured  with  the  resonator  in  air,  while  QL  is  the  Q  value  measured 
when  the  phantom  was  placed  centrally  inside  the  resonator.  The  Q  damping  is  defined  as 
the  ratio  of  Q0  /Qt .  The  results  are  the  means  of  three  measurements. 


Table  7.  Results  of  the  laboratory  bench  testing  of  the  slotted  tube 
resonator.  Q0  and  QL  were  measured  with  the  250  ml  phantom. 


Resonant 

Mode 

Resonant 

Frequency 

(MHz) 

Unloaded 

Qo 

Loaded 

Ql 

Qo/Ql 

19p 

188.22 

97.5 

54.2 

200.16 

63.8 

53.6 

1.19 

As  can  be  seen  from  Table  7,  the  unloaded  Q0  for  l9F  is  much  higher  than  that  for 
'H,  and  the  loaded  QL  for  both  resonant  frequencies  is  essentially  same.  The  Q  damping, 
Qo  /QL,  for  19F  is  1.5  times  higher  than  that  for  'H.  This  means  that  the  RF  lossy  phantom 
has  a  much  stronger  effect  on  the  I9F  resonance  than  on  the  'H  resonance.  The  high  Q 
values  could  be  the  result  of  several  factors,  including  a  relatively  small  usable  volume, 
minimized  length  of  connecting  wires,  and  a  structure  without  a  copper  shield.  The  high 
Q  values  guarantee  that  the  resonator  has  a  higher  excitation  efficiency,  i.  e.,  the  90°  pulse 
width  for  a  given  RF  power  level  is  short  or  the  magnitude  of  the  Bi  field  generated  with 
a  given  RF  power  level  is  high.  The  results  of  the  laboratory  bench  testing  also  indicated 
that  the  resonator  was  easy  to  tune  and  match  and  both  resonant  frequencies  were  stable 
and  immune  to  external  electromagnetic  interference. 
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•  Phantom  Imaging 

Once  again,  the  NMR  performance  of  the  slotted  tube  resonator  was  tested  on  the 
Omega  CSI  4.7  T  magnet.  The  in  vitro  resonator  sensitivity  for  proton  and  fluorine  was 
determined  using  the  250  ml  phantom,  which  was  filled  with  a  solution  of  TFA  (5%v/v) 
and  CuS04  (lg/liter)  in  distilled  water.  Measurements  of  90°  pulse  width  were  performed 
from  a  180°  null  of  the  whole  sample.  Typical  90°  pulse  widths  were  around  170  ps  for 
proton  and  1 10  ps  for  fluorine.  Shimming  was  performed  on  the  water  proton  FID  of  the 
saline  solution  with  resonator  tuned  to  the  ’H  resonance,  to  a  typical  spectral  linewidth  of 
32  Hz.  The  driven  equilibrium  spin-echo  pulse  sequence  was  used  for  the  imaging 
experiment.  Imaging  parameters  were:  TR=  150ms,  TE  =  8ms,  and  matrix  size  =128  x  64 
x  8.  The  field  of  view  (FOV)  was  80  mm  x  80  mm  for  transaxial  images,  100  mm  x  100 
mm  for  coronal,  and  200  mm  x  200  mm  for  sagittal  images.  The  images  were  acquired 
with  one  excitation,  thus  giving  a  total  acquisition  time  of  1. 16  min.  The  image  signal-to- 
noise  ratios  (SNRs)  were  measured  as  described  before.  The  results  of  SNR 
measurements  for  the  proton  images  are  listed  in  Table  8. 


Table  8.  Results  of  SNR  measurements  for  the  proton  images. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower- 

Left 

SNR 

Lower- 

Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

389.3 

199.5 

897.8 

666.3 

538.2 

177.1 

Sagittal  Section 

494.2 

517.5 

478.8 

450.5 

485.3 

28.1 

Coronal  Section 

190.3 

30.2 

164.6 

22.8 

102.0 

87.9 

Following  the  'H  imaging  experiment,  the  resonator  was  retuned  in  place  to 
188.27  MHz  and  corresponding  19F  imaging  experiment  was  carried  out.  The  90°  pulse 
was  108  ps  and  the  spectral  linewidth  after  shimming  was  46  Hz.  The  spin-echo  (SE) 
images  were  acquired  with  the  imaging  parameters:  TR  =  150  ms,  TE  =  8  ms,  FOV  =  100 
x  100  mm,  NA  =  16,  and  matrix  size  =  64  x  32.  The  image  SNR  measurements  were 
performed  in  the  same  way  as  the  proton  images  and  the  results  are  given  in  Table  9. 


Table  9.  Results  of  SNR  measurements  for  the  fluorine  images. 


Upper- 

Left 

SNR 

Upper- 

Right 

SNR 

Lower 

-Left 

SNR 

Lower 

-Right 

SNR 

Mean 

SNR 

SD 

Transaxial  Section 

3.1 

2.9 

2.6 

3.0 

2.9 

0.2 

Sagittal  Section 

4.2 

5.1 

4.0 

4.2 

4.4 

0.5 

Coronal  Section 

4.2 

4.3 

4.7 

4.6 

4.5 

0.2 
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4)  Assessment  of  Hexafluorobenzene  (HFB)  Distribution 
Animal  Preparation 

Once  the  tumors  reached  ~1  cm  diameter  (~  0.5  cm3),  corresponding  to  a  typical 
lower  limit  of  tumor  detected  in  patients,  the  rats  were  anesthetized  with  200  pi  ketamine 
hydrochloride  i.p.  (100  mg/ml;  Aveco,  Fort  Dodge,  1A)  and  were  maintained  under 
general  gaseous  anesthesia  using  a  small  animal  anesthesia  unit  with  air  (1.0  dm3/min) 
and  1.0%  isoflurane  (Ohmeda  PPD  Inc.,  Fort  Dodge,  IA).  Tumor  hair  was  cut  with  a  pair 
of  surgical  scissors  for  reduction  of  the  NIR  light  scattering  and  ease  of  HFB  injection. 
The  rats  were  placed  on  their  sides  in  a  specially  designed  NMR  bed.  The  body 
temperature  was  maintained  at  about  37°C  by  a  warm  water  blanket  with  a  feedback 
system  (K-MOD  100,  Baxter  Healthcare  Co.,  Deerfield,  IL).  A  fiber  optic  pulse  oximeter 
(Nonin  Medical,  Inc.,  Plymouth,  MN)  was  placed  on  the  hind  foot  to  monitor  arterial 
hemoglobin  saturation  (sa02)  and  heart  rate  (HR),  and  a  thermocouple  (Cole-Parmer 
Instrument  Co.,  Vernon  Hills,  IL)  was  inserted  rectally  to  monitor  core  temperature.  40  pi 
HFB  (99.9%,  Aldrich  Chemical  Co.,  St.  Louis,  MO)  was  injected  directly  into  selected 
areas  of  tumor  central  and  peripheral  regions  at  the  same  plane  using  a  Hamilton  syringe 
with  a  32  G  needle.  The  needle  was  inserted  manually  to  penetrate  across  the  whole 
tumor  and  withdrawn  ~  1  mm  to  reduce  the  tissue  pressure  and  2~5  pi  HFB  injected.  The 
needle  was  repeatedly  withdrawn  at  a  step  size  1~2  mm  and  further  HFB  injected.  A  total 
of  2  ~  3  tracks  of  HFB  was  injected  in  the  form  of  a  fan. 


90°  180°  180°  90° 


Figure  7.  Pulse  diagram  of  the  driven  equilibrium  spin-echo  pulse  sequence 
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JH  Magnetic  Resonance  Imaging 

For  each  tumor,  high  resolution  3-D  *H  MR  imaging  was  first  performed  for 
anatomical  reference.  A  frequency  switchable  (*H/19F)  RF  coil  was  placed  around  the 
tumor  with  the  tumor  centered  in  the  coil.  The  rat  was  then  be  placed  on  its  side  in  the 
bed  and  positioned  in  the  isocenter  of  the  magnet,  with  the  coil  axis  being  aligned 
perpendicular  to  the  direction  of  the  magnet's  static  magnetic  field  B0.  The  90°  pulse 
width  was  determined  from  a  180°  null  of  the  whole  tumor.  Shimming  was  performed  on 
the  tumor  tissue  water  proton  FID  to  a  typical  line  width  of  70  Hz.  This  was  accomplished 
by  adjusting  fourteen  shim  currents  on  the  FID  through  the  system  automated  shimming 
utility  until  a  symmetrical  and  long  FID  was  obtained.  The  goal  was  to  obtain  a  B0  of 
homogeneous  strength  over  the  entire  tumor. 


Table  10.  ’H  imaging  parameters 


Repetition  Time  (TR) 

1500  ms 

Echo  Time  (TE) 

80  ms 

Matrix  Size 

128x64x8 

Field  of  View  (FOV) 

20x20x20  mm 

Digital  Resolution 

156  pm  x  312  pm  x  2.5  mm 

Number  of  Acquisitions  (NA) 

2 

Total  Acquisition  Time 

25  min  36  sec 

To  minimize  acquisition  time,  the  driven-equilibrium  spin-echo  pulse  sequence 
with  a  short  echo  time  (TE  =  80  ms)  (Figure  7)  was  used  so  that  the  T2  dephasing  effects 
was  not  significant.  *H  images  were  acquired  as  a  3-D  date  set.  Imaging  parameters  used 
are  shown  in  Table  10.  These  T2- weighted  images  showed  the  tumor  anatomy  and  its 
position  relative  to  the  back  of  the  rat.  Following  the  *H  imaging,  the  corresponding  19F 
imaging  was  then  performed  to  show  the  distribution  of  HFB  in  the  tumor. 

19F  MRI  and  Assessment  of  HFB  Distribution 

The  magnet  frequency  was  then  set  on  resonance  for  the  CF  group  of  HFB  and  the 
coil  was  retuned  in  place  to  188.273  MHz.  Corresponding  fluorine  density-weighted 
images  was  acquired  as  a  3-D  data  set  using  the  driven-equilibrium  spin-echo  pulse 
sequence.  The  imaging  parameters  are  shown  in  Table  11.  Gradients  were  compensated 
to  account  for  the  difference  in  gyromagnetic  ratios.  Data  were  processed  using  sine-bell 
apodization  to  improve  SNR  and  zero-filling  in  the  first  phase  encoding  direction  for  the 
execution  of  fast  Fourier  transform  (FFT).  Images  were  transferred  to  a  PC  and  further 
processed  off  line  using  a  SCION  imaging  software.  Figure  8  shows  conventional  SE  *H 
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images  (top  four)  and  corresponding  19F  SE  images  of  a  representative  breast  tumor. 
Comparison  of  the  *H  and  19F  images  reveals  that  the  tumor  was  centrally  labeled  in  this 
case. 


Table  11.  19F  imaging  parameters 


Repetition  Time  (TR) 

150  ms 

Echo  Time  (TE) 

8  ms 

Matrix  Size 

128x64x8 

Field  of  View  (FOV) 

20x20x20  mm 

Digital  Resolution 

156  pm  x  312  pm  x  2.5  mm 

Number  of  Acquisitions  (NA) 

8 

Total  Acquisition  Time 

10  min  15  sec 

Figure  8.  ‘H  and  19F  coronal  images  of  a  representative 
breast  tumor.  Bar  represents  1  cm. 


Determination  of  HFB  Redistribution  and  Clearance 

In  contrast  to  traditional  perfluorocarbons  (PFCs),  which  exhibit  excessive  tissue 
retention,  HFB  clears  from  tissue  relatively  rapidly.  To  verify  that  HFB  global  clearance 
was  not  so  rapid  as  to  interfere  with  relaxometry,  19F  NMR  spectroscopy  was  performed 
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interleaved  between  19F  EPI  relaxometry.  HFB  spectroscopic  signal  intensity  was 
compared  with  a  benzene-doped  standard.  In  addition,  19F  3-D  MRI  was  also  repeated  to 
determine  the  global  and  regional  clearance  of  HFB  and  to  assess  any  redistribution  of 
signal  intensity  within  the  tumor.  To  compute  HFB  clearance  time  (biological  halflife),  a 
single  parameter  exponential  decay  model  was  used: 

e(0  =  g(0)exp(-i)  (14) 

where  Q(t)  was  the  signal  intensity  at  any  given  time  t,  (7(0)  was  the  signal  intensity  at 
time  zero,  and  T  was  the  clearance  time  in  minute.  Since  HFB  is  a  non-ionic  freely 
diffusable  tracer,  clearance  provides  an  indication  of  relative  tumor  blood  flow  (TBF).  . 


Figure  9,  HFB  clearance  curve  of  a  typical  voxel 
from  a  13762NF  breast  tumor,  t  =  160  min. 


Figure  9  shows  the  HFB  clearance  curve  of  a  representative  voxel  from  a 
13762NF  breast  tumor  during  the  3  hours  period  of  a  typical  experiment.  For  many 
regions,  the  HFB  signal  intensity  was  found  to  decay  exponentially  with  a  typical 
biological  half-life  ranging  from  Ti/2  =  700  to  1200  min.  Since  our  19F  MR  EPI  oximetry 
experiments  lasted  about  three  hours,  this  means  the  global  and  regional  clearance  and 
redistribution  of  HFB  within  the  tumors  did  not  interfere  with  19F  MR  EPI  oximetry.  By 
analyzing  the  clearance  of  individual  voxels,  it  was  found  that  as  a  whole  19F  signal 
intensity  decreased  with  time,  some  voxels,  however,  showed  increases  in  signal  intensity 
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as  shown  in  Figure  10.  This  could  be  due  to  the  inflow  of  HFB  from  the  surrounding 
voxels.  Since  HFB  concentration  only  affects  the  SNR  of  19F  signal  intensity  and  does 
not  have  a  direct  impact  on  relaxation  rate  Rl,  this  local  redistribution  of  HFB  will  not 
compromise  the  quality  of  our  tumor  oxygenation  studies. 


Figure  10.  HFB  clearance  curve  of  a  voxel  from  a  13762NF  breast 
tumor  that  showed  an  increase  in  signal  intensity  with  time 


5)  Software  Development 

A  Computer  Program  to  Assess  the  Goodness  of  T1  Relaxation  Curve  Fit 

19F  EPI  was  used  as  a  basic  building  block  for  19F  MR  EPI  oximetry.  To  measure 
spin-lattice  relaxation  time  (Tl)  and  hence  p02  maps,  a  pulse  burst  saturating  (PBSR) 
pulse  train  was  applied  immediately  prior  to  the  EPI  sequence.  This  led  to  some  loss  of 
dynamic  range  in  the  results,  but  substantially  shortened  the  time  required  for  a  Tl 
experiment  and  thus,  improved  the  temporal  resolution.  Fourteen  data  sets  were  acquired 
with  delay  time  in  a  geometric  progression  ranging  from  200  ms  to  90  sec.  Since  the 
longest  and  shortest  delays  were  alternated  during  data  acquisitions  to  reduce  systematic 
bias,  it  was  necessary  to  restore  the  images  to  the  correct  order,  i.e.,  from  the  shortest  to 
the  longest  delay  (200  ms — »  90  sec)  prior  to  the  curve  fit 

In  order  to  obtain  the  spatial  distribution  (map)  of p02in  breast  tumors,  we  needed 
to  solve  for  the  Rl  map  of  injected  HFB  first.  This  was  done  by  fitting  the  signal  intensity 
in  each  of  the  voxels  of  the  fourteen  images  to  a  three  parameter  exponential  relaxation 
model: 
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(15) 


y„ O'.  J)  -  A(i,  j)  •  [1  -  (1  +  W)  •  exp(-i?l(7, y)  •  r„ )] 

(«  =  1, 2, . ,14) 

(/,  7  =  1,2, . 32) 

by  the  Levenberg-Marquardt  least-squares  algorithm,  where  yn(i,j)  was  the  measured 
signal  intensity  corresponding  to  delay  time  r„(the  «th  image)  of  voxel  (ij),  A(i,j) 
was  the  fully  relaxed  signal  intensity  amplitude  of  voxel  (ij),  W  was  a  dimensionless 
scaling  factor  allowing  for  imperfect  signal  conversion,  and  Rl  (i,j)  was  the  relaxation 
rate  of  voxel  (/,  j)  in  unit  of  sec'1 . 


,5  T1  Relaxation  Curve  of  a  Well-Oxygenated  Voxel 


T1  Relaxation  Curve  of  a  Hypoxic  Voxel 


Figure  11.  T1  relaxation  curves  of  a  well-oxygenated 
voxel  (left)  and  a  hypoxic  voxel  (right). 


Good  T1  relaxation  curve  fit  was  crucial  in  19F  PBSR-EPI  oximetry.  It  was 
relatively  easy  to  obtain  excellent  curve  fits,  provided  SNR  >  10  for  most  intense  signals. 
To  assess  the  goodness  of  curve  fit,  I  wrote  a  Unix  based  NMR-shell  script  program  to 
display  the  T1  relaxation  curve  of  individual  voxels  of  PBSR-EPI  images.  The  graphical 
representation  of  the  data,  relaxation  model,  and  goodness  of  fit  provided  a  quick  and 
convenient  way  of  assessing  the  quality  of  the  data.  Based  on  the  rms  error  of  fit,  the 
goodness  of  curve  fit  could  be  classified  into  three  categories:  good,  intermediate,  and 
poor. 


For  typical  good  relaxation  data,  the  global  least  squares  minimization  was 
achieved  with  a  very  small  rms  error.  These  data,  related  to  strong  19F  signal,  represented 
those  voxels  within  the  region  of  interest  (ROI).  For  typical  intermediate  data,  the  global 
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least  squares  minimization  was  achieved  with  a  relatively  large  rms  error  due  to  a  single 
bad  data  point,  which  led  to  a  big  T1  error.  To  improve  the  relaxation  curve  fit,  a 
threshold  was  applied  to  the  raw  data,  which  eliminated  those  data  points  >  3cr  from  the 
curve.  The  remaining  data,  thus,  provided  a  better  curve  fit.  In  applying  the  threshold, 
only  one  data  point  was  eliminated  from  the  raw  data  per  curve.  These  relaxation  data 
represented  those  voxels  within  the  ROI,  but  had  relatively  strong  19F  signals.  In  general, 
they  gave  good  T1  values.  For  poor  relaxation  data,  the  global  least  squares  minimization 
either  failed  or  converged  slowly  with  a  very  large  rms  error.  These  relaxation  data, 
associated  with  those  voxels  outside  the  ROI,  represented  either  very  week  19F  signals  or 
background  noise.  They  gave  unreliable  T1  values.  Figure  11  shows  T1  relaxation  curves 
of  a  well-oxygenated  voxel  (left)  and  a  hypoxic  voxel  (right). 

A  Computer  Program  to  Compute  HFB  Clearance  Rate 

To  compute  HFB  clearance  rate  (biological  halflife),  I  used  a  single  parameter 
exponential  decay  model: 


Qi,i  (0  =  Qt,i  (0)exp(-*  ITUI)  (16) 

where  Qt  J  (t)  was  the  signal  intensity  of  voxel  (/,  j)  at  time  t,  Qt  j  (0)  was  the  signal 
intensity  of  voxel  (/',  j)  at  time  zero,  and  Ti}  was  the  clearance  rate  of  voxel  (/,  j). 
Equation  (16)  was  transformed  into  a  logarithmic  form: 

ln(Q,  j  (0)  =  ln(<2,j  (0))  ~  t  /  TUj  (17) 

Based  on  above  theory,  1  wrote  a  program  to  compute  the  clearance  rate  TUj  of 

individual  voxels  by  fitting  the  signal  intensity  of  PBSR-EPI  images  to  Equation  (17).  A 
clearance  rate  map  was  then  produced  by  displaying  TUf  as  a  color-coded  image.  This 

program  is  capable  of  determining  the  global  and  regional  clearance  rate  of  HFB. 

A  Computer  Program  to  Process  p02  Data 

A  windows-based  Visual  Basic  program  was  written  to  process  p02  data.  The 
program  can  display  and  process  images  in  the  following  formats: 


BMP 

Bitmap 

GIP 

Graphics  Interchange  Format 

JPG 

Joint  Photographic  Experts  Group 

DIB 

Device  Independent  Bitmap 

WMF 

Windows  MetaFile 

EMF 

Enhanced  MetaFile 

ICO 

Icons 
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More  functions  are  being  added  to  the  program. 

A  Computer  Program  to  Process  NIRS  Data 

NIRS  signal  and  biomedical  signals  in  general  are  commonly  contaminated  by  a 
number  of  sources.  To  reduce  the  contamination  (noise),  I  employed  digital  signal 
processing  techniques:  averaging  and  filtering.  Averaging  was  to  improve  the  SNR  of  the 
signal  and  filtering  was  to  remove  the  unwanted  frequency  components  contained  in  the 
signal.  Digital  filters  are  mathematical  algorithms  implemented  in  hardware  and/or 
software  that  operate  on  a  digital  signal  to  produce  a  desired  digital  output  signal.  I  wrote 
a  LabView™  program  to  implement  these  techniques.  Lab  View™  is  a  graphical 
programming  development  environment  based  on  the  G-programming  language  for  data 
acquisition  and  control,  data  analysis,  and  data  presentation.  It  is  more  powerful  and,  at 
the  same  time,  easier  to  learn  and  implement  than  the  traditional  text-based  programming 
languages,  such  as  C. 


•  Averaging 

The  first  step  in  processing  the  NIRS  data  was  averaging.  Four  adjacent  data 
points  were  averaged  to  produce  a  new  data  point  according  to: 


y{n)  = 


x(n)  +  x  (n  + 1)  +  x(n  +  2)  +  x(n  +  3) 

”T~ 


(18) 


where  x{n)  represented  the  nth  input  data  point  and  y(n)  the  nth  output  data  point.  Since 
the  NIRS  data  were  acquired  at  relatively  high  sampling  frequency,  this  averaging 
operation  improved  the  SNR  of  the  signal  but  not  at  the  expense  of  losing  detailed 
information  in  the  data. 


•  Digital  Filtering 

Following  averaging,  the  data  were  fed  into  a  high  order  Butterworth  low-pass 
filter  with  a  cut-off  frequency  of  1.0  Hz.  Butterworth  filters  are  one  of  the  infinite 
impulse  response  (HR)  filters.  The  main  characteristics  of  the  IIR  filters  is  that  the  output 
depends  not  only  on  the  current  and  past  input  data,  but  also  on  the  past  output  data  as 
shown  in  the  following  recursive  equation  [68]: 

N  M 

y(n)=Yjakx(n-k)-Yjbky(n-lc )  (19) 

k= 0  k=l 

where  ak  are  called  the  forward  coefficients,  bk  the  reverse  coefficients,  x(n)  and  y(n) 

are  the  input  and  output  to  the  filter,  and  N  and  M  are  the  number  of  forward  and  reverse 
coefficients,  respectively.  Given  the  filter  specifications  (Table  12),  ak,  bk,  N  and  M 

were  computed  using  either  a  standard  procedure. 
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The  rationale  for  me  to  use  a  high  order  low-pass  Butterworth  filter  was  two-fold. 
Firstly,  a  Butterworth  filter  requires  fewer  coefficients.  Thus,  it  executes  faster  and  does 
not  require  extra  memory.  Secondly,  a  Butterworth  filter  has  no  ripple  in  either  the 
passband  or  the  stopband,  and  has  a  smooth,  monotonically  decreasing  frequency 
response  in  the  transition  band. 

Table  12.  Specifications  for  the  low-pass  Butterworth  filter 


Parameters 

Specifications 

Passband 

0  -1.0  Hz 

Stopband 

>2.0  Hz 

Stoppband  attenuation 

>20  dB 

Filter  Order 

15 

6)  Investigation  of  Tumor  Oyxgenation  during  Untreated  Growth 
Investigation  of  Tumor  Tissue  pC>2 
•  19F  MR  EPI  Oximetry 

EPI  is  a  high-speed  imaging  technique  first  developed  by  Mansfield  [69],  EPI  can 
be  divided  into  single-shot  EPI  and  multiple-shot  EPI.  For  this  project,  I  used  a  variation 
of  the  original  single-shot  EPI,  pulse  burst  saturation  recovery  echo-planar  imaging 
(PBSR-EPI).  It  has  been  shown  to  be  very  effective  and  efficient  in  oxygen  dynamic 
studies  [70-72].  A  typical  PBSR-EPI  pulse  sequence  diagram  is  shown  in  Figure  12. 
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Figure  12.  A  typical  PBSR-EPI  pulse  sequence  diagram. 


34 


The  PBSR-EPI  consists  of  three  basic  building  blocks:  a  preparation  pulse 
sequence  (PBSR),  a  variable  delay  time  x,  and  a  single  shot  spin-echo  EPI  with  “blipped” 
phase  encoding  gradient.  The  PBSR  is  a  pulse  train  of  20  non-spatially  selective  90° 
pulses  with  50  ms  spacing  to  saturate  19F  nuclei.  Fourteen  different  x  values,  ranging 
from  200  ms  to  90  s,  are  used  to  modulate  image  intensities  so  that  good  T1  relaxation 
curves  can  be  obtained.  The  function  of  the  brief  “blipped”  phase  encoding  gradient 
between  each  echo  is  to  increment  phase  in  the  ky  direction  to  form  a  blipped  echo  image. 
Fourteen  EPI  images  are  acquired  with  32x32  in-plane  resolution.  PBSR-EPI  may  result 
in  some  loss  of  dynamic  range  in  signal  intensity,  but  substantially  shortens  the  data 
acquisition  time  for  a  T1  experiment.  By  incorporating  a  PBSR  preparation  pulse 
sequence  into  EPI,  each  R1  (1/T1)  map  can  be  produced  with  a  temporal  resolution  of  ~  6 
min,  facilitating  measurements  of  dynamic  changes  in  p02  accompanying  therapeutic 
interventions  and  allowing  the  fate  of  individual  voxels  to  be  traced. 


Figure  13.  A  schematic  representation  of  the  ARDVARC  protocol. 


To  enhance  SNR  of  the  19F  EPI  images,  and  thus,  improve  the  precision  of  R1 
measurements,  I  used  a  novel  data  acquisition  protocol  developed  in  our  lab:  ARDVARC 
(Alternated  R1  Delays  with  Variable  Acquisitions  to  Reduce  Clearance  effects)  (Figure 
13).  The  ARDVARC  has  two  new  features  in  comparison  with  the  traditional  approach. 
Firstly,  the  number  of  acquisitions  (NA)  varies  with  x  values.  This  innovation 
significantly  improves  the  SNR  for  short  delays  and  provides  a  better  curve  fit  by 
additional  acquisitions.  Data  are  amplitude-corrected,  /.<?.,  divided  by  the  number  of 
acquisitions  to  maintain  correct  signal  amplitude.  Secondly,  longest  and  shortest  delays 
are  alternated  to  reduce  bias  resulting  from  HFB  clearance  or  fluctuations  in  R1  (Table 
13). 
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In  order  to  determine  Rl„  and  hence  p02,  signal  intensity  in  each  voxel  was  fitted 
to  a  three-parameter  exponential  relaxation  model: 

y„  0‘,  j )  =  A(U  j)-[\-{\  +  W  (/,  j))  •  exp(-/?l(/,  j)  ■  t„  )] 

(n  ~  1,  2, . ,14)  (20) 

0,  7  =  1,2, . 32) 

by  the  Levenberg-Marquardt  least-squares  algorithm,  where  y„(i,  j )  was  the  measured 
signal  intensity  of  voxel  (/,  j)  at  delay  time  x„  (the  nth  image),  A(i,j)  was  the  fully  relaxed 
signal  intensity  of  voxel  (/,  j),  W{i,  j )  was  a  dimensionless  scaling  factor  used  for 
imperfect  signal  conversion,  and  Rl(i,  j)  is  the  relaxation  rate  of  voxel  (/,  j)  in  unit  of 
sec'1. 


Table  13.  ARDVARC  parameters 


Order  of 
Acquisition 

Order  in 
Curve  Fit  (n) 

Delay  Time 
x  (sec) 

Number  of 
Acquisitions 
(NA) 

i 

14 

90 

1 

2 

1 

0.2 

12 

3 

13 

60 

1 

4 

2 

0.4 

12 

5 

12 

40 

1 

6 

3 

0.6 

12 

7 

11 

20 

1 

8 

4 

0.8 

12 

9 

10 

16 

2 

10 

5 

1 

8 

11 

9 

8 

2 

12 

6 

1.5 

4 

13 

8 

4 

4 

14 

7 

2 

4 

Regional  tumor  p02  maps  were  produced  by  applying  the  calibration  curve  to  Rl 
maps.  It  was  found  that  at  37C°  and  4.7  T: 


P02(}J ) 


Rl(i,j)-  0.0836 
0.00188 


(torr) 


(21) 


using  PBSR-EPI  for  HFB. 

•  Respiratory  Challenge  Paradigms 
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One  of  the  important  issues  in  radiotherapy  is  how  to  manipulate  oxygenation  in 
tumors  and  monitor  its  dynamic  response  to  various  therapeutic  interventions 
noninvasively.  Many  protocols  and  techniques  have  been  proposed.  A  simple 
intervention  is  respiratory  challenge,  i.e.,  attempting  to  elevate  tumor  oxygenation  with 
inhaled  gas.  For  this  project,  two  respiratory  challenge  paradigms  were  used  to 
manipulate  oxygenation  in  tumors.  19F  MR  EPI  and  NIRS  were  used  to  monitor  oxygen 
dynamic  response  and  assess  its  temporal  characteristics  and  changes  in  the  extent  of 
hypoxia. 


(1) 

^4/r 

(21%  02) 

-»  Carbogen  —> 

(95%  02  +  5%  02) 

Air  — > 

Carbogen 

— »  Air 

(2) 

Air 

->  Carbogen  -» 

Air  -» 

100%  02 

— >  ^4ir 

Most  experiments  were  performed  using  paradigm  (1).  However,  in  some  cases, 
repeated  carbogen  interventions  were  performed  sequentially  to  evaluate  the 
reproducibility  of  the  time  course  profiles  of  the  tumors. 

•  Investigation  of  Tumor  Tissue  Oxygen  Dynamics  by  19F  MR  EPI 

Following  *H  and  19F  MRI  to  determine  the  tumor  anatomy  and  the  distribution  of 
HFB,  a  series  of  19F  MR  EPI  oximetry  experiments  were  performed  according  to  the 
respiratory  challenge  paradigms  described  above.  During  the  imaging  sessions,  initially, 
rats  inhaled  medical  grade  compressed  air  (21%  02)  and  three  measurements  were  taken 
for  the  first  phase.  Subsequently,  five  measurements  were  taken  for  each  gas  switch, 
giving  a  total  of  twenty-three  measurements.  The  complete  five-phase  imaging  session 
took  about  2.5  hours. 


Figure  14.  *H  and  19F  coronal  images  of  a  breast  tumor. 
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Figure  14  shows  conventional  spin  echo  (SE)  'H  images  (top  four)  and 
corresponding  19F  SE  images  of  a  representative  breast  tumor.  ]H  images  were  acquired 
to  show  the  tumor  anatomy  and  19F  images  to  show  the  distribution  of  HFB  within  the 
tumor.  In  this  case,  the  tumor  was  centrally  labeled  Figure  15  shows  representative  p02 
maps  obtained  from  a  rat  mammary  adenocarcinoma  13762NF  using  19F  MR  EPI 
oximetry.  Each  map  was  acquired  in  6.5  min,  with  (A)  being  rat  breathing  33%  02  and 
(B)  carbogen  (95%  02  +  5%  C02). 


Figure  15.  Representative  p02  maps  of  a  breast  tumor. 


pO 2  (torr) 


Figure  16.  A  representative  histogram  for  a  breast  tumor. 

Figure  16  shows  histograms  of  p02  distributions  in  response  to  respiratory 
challenge  in  a  rat  mammary  adenocarcinoma  13762NF.  Elevating  inspired  02  (carbogen 
in  this  case)  caused  a  substantial  change  in  the  p02  distribution,  resulting  in  a  shift  toward 
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the  higher  p02  values.  Arrows  indicate  mean  (x)  p02.  Blue  histogram  represents  baseline 
33%  02,  and  red  histogram  carbogen.  Figure  17  shows  dynamic  changes  in  breast  tumor 
mean  p02  corresponding  to  Figure  15.  Each  data  point  represents  mean  p02  ±  SD.  The 
mean  baseline  p02  value  was  18.1  ±2.1  torr  when  the  rat  was  breathing  33%  02  and  was 
108  ±  42  torr  after  the  gas  was  switched  to  100%  02  (p  <  0.009).  The  p02  dropped  to  75 
±19  torr  after  the  gas  was  switched  back  to  33%  02.  During  the  second  phase  of 
respiratory  challenge,  the  p02  increased  to  110.9  ±  54  torr  (p< 0.02)  when  the  gas  was 
switched  to  carbogen  and  dropped  to  32.4  ±  33  torr  after  the  gas  was  switched  back  to 
33%  02.  As  can  be  seen  from  Figure  17,  19F  MR  EPI  oximetry  showed  good 
reproducibility. 


Figure  17.  Dynamic  response  in  tumor  tissue p02. 


Figure  18.  Maximum  p02  vs.  mean  baseline  p02  for  a  breast  tumor 
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A  strong  linear  correlation  has  been  observed  in  some  of  13762NF  breast  tumors 
between  initial  mean  baseline  p02  and  the  maximum  p02  for  a  given  group  of  voxels 
with  respect  to  carbogen  or  oxygen  breathing.  For  the  breast  tumor  shown  in  Figure  18, 
an  ROI  of  83  voxels  was  selected.  In  this  case,  the  linear  coefficient  r  was  found  to  be 
0.95.  This  observation  was  of  clinical  significance  and  could  help  us  predict  tumor 
response  to  elevated  oxygen  breathing  based  on  its  initial  baseline  value.  Figure  19  shows 
a  group  of  1 1  voxels  of  a  second  breast  tumor.  Once  again,  a  strong  linear  correlation  was 
observed  and,  in  this  case,  r  was  found  to  be  0.86. 


Figure  19.  Maximum  p02  vs.  mean  baseline  p02  for 
a  group  of  1 1  voxels  of  a  second  breast  tumor. 


Figure  20.  Dynamic  changes  in  p02  for  a  group 
of  1 1  voxels  of  a  second  breast  tumor. 
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The  19F  MR-EPI  oximetry  of  tumor  has  the  distinct  advantage  over  other 
techniques  that  subsequent  measurements  are  completely  non-invasive.  The  greatest 
strength  of  this  method  is  the  ability  to  trace  the  fate  of  individual  voxels  (regions)  with 
respect  to  therapeutic  interventions.  Figure  20  shows  dynamic  changes  in  p02  of  11 
specific  voxels  of  a  second  breast  tumor  with  respect  to  different  inhaled  gases.  It  is 
noteworthy  that  voxels  with  high  baseline  p02  had  significantly  different  response 
characteristics  from  those  with  initially  low  p02i  which  showed  small  changes.  To  further 
investigate  the  temporal  response  characteristics  of  individual  voxels,  we  modeled  the 
temporal  response  in  p02  using  exponential  equations: 

1 )  y  =  a  +  b-(l-e~'/T),  for  increasing  trend 

2)  y  =  a+b-  e~“T  ,  for  decreasing  trend 

where  y  is  p02 ,  a  and  b  are  two  constants,  t  is  time,  and  t  is  the  time  constant. 


Mean  BaseHne  pQ2  (torr) 


Figure  21.  Time  constant  vs.  mean  baseline  p02  of  44 
specific  voxels  of  a  third  breast  tumor. 


Figure  21  shows  the  relationship  between  time  constant  and  mean  baseline  p02  of 
44  specific  voxels  of  a  third  breast  tumor  (r  =  0.68).  It  was  found  that,  in  general,  the  time 
constants  of  well-oxygenated  voxels  (10  ~  20  min)  were  much  shorter  than  those  of 
hypoxic  voxels  (>  50  min)  and  the  global  p02  time  constant  (60  ~  80  min)  was  much 
longer  than  the  vascular  hemoglobin  saturation  (s02)  time  constant  (10  ~  20  min).  In 
general,  changes  in  tumor  vascular  s02  preceded  tumor  tissue  pQz?  particularly  for 
smaller  tumors. 


41 


Investigation  of  Tumor  Vascular  [Hb02]  or  SO2  and  [Hbjxotai 
•  Oximetry  Based  on  Near-Infrared  Spectroscopy  (NIRS) 

Near-Infrared  Spectroscopy  (NIRS)  can  be  used  to  measure  total  hemoglobin 
concentration  [Hb]Totai,  oxyhemoglobin  concentration  [Hb02],  and  oxygen  saturation 
(s02)  because  of  two  things:  firstly,  the  absorption  of  light  by  deoxyhemoglobin  (Hb)  and 
oxyhemoglobin  (Hb02)  predominates  over  water  and  other  macromolecules  at  the 
selected  wavelengths;  secondly,  the  absorption  coefficient  of  deoxyhemoglobin  (Hb) 
differs  substantially  from  that  of  oxyhemoglobin  (Hb02)  in  the  NIR  region  (700  -  900 
nm).  The  traditional  continuous  wave  spectroscopy  (CWS)  uses  the  classical  Beer- 
Lambert’s  law  to  determine  the  absorption  coefficient  of  the  sample  over  a  known  optical 
path  length  and  thus  the  concentration  of  the  sample.  Current  available  CWS  techniques 
include  dual-wavelength  and  multi-wavelength  light  spectroscopy,  pulse  oximetry,  and 
CW  brain  functional  imaging  systems. 

Theoretically  speaking,  however,  Beer-Lambert’s  law  is  valid  only  for 
homogenous  translucent  solutions  that  do  not  scatter  light.  In  a  continuous  medium,  light 
photons  travel  at  the  speed  of  light  that  is  dependent  only  on  the  refractive  index  of  the 
medium.  If  this  is  a  scattering  medium,  photons  travel  with  the  constant  speed  and 
direction  until  they  collide  elastically  with  a  scatterer,  resulting  in  a  change  in  traveling 
direction.  Depending  on  the  scattering  properties  of  the  medium,  the  scattered  photons 
may  travel  in  a  random  direction,  or  in  a  preferential  forward  or  backward  direction.  The 
total  distance  or  optical  path  length  traveled  by  photons  between  the  source  and  detector 
is  much  longer  than  the  geomebic  distance  between  the  source  and  detector.  In  this  case, 
Beer-Lamberf  s  law  is  not  valid  since  exact  L  is  not  known.  In  the  absence  of  scatterers, 
the  total  optical  path  length  traveled  by  photons  before  detection  is  simply  the  geometric 
distance  between  the  source  and  detector.  Thus,  L  is  known  and  Beer-Lamberf  s  law  is 
valid.  Since  biological  tissues  are  highly  inhomogeneous  and  optically  turbid,  it  is 
obvious  that  we  can  not  directly  apply  Beer-Lambert’s  law  to  biological  tissues.  Light 
attenuation  results  not  only  from  tissue  absorption,  but  also  from  tissue  scattering.  Light 
absorption  results  mainly  from  absorption  of  oxyhemoglobin  and  deoxyhemoglobin, 
myoglobin,  and  cytochrome  oxidase.  The  amount  of  light  absorption  depends  on  the 
concentrations  of  these  molecules.  Light  scattering  results  from  mitochondrion,  protein, 
and  various  ions.  Scattering  gives  rise  to  distributed  path  lengths  that  are  significantly 
longer  than  the  geometric  distance  between  the  source  and  detector  used  in  Beer- 
Lambert’s  law,  leading  to  a  distorted  absorption  spectrum  of  the  sample.  To  take 
scattering  into  account,  we  need  to  use  the  photon  diffusion  approximation  theory  [73] 
based  on  a  general  mass  transfer  diffusion  equation.  This  equation  allows  us  to  describe 
photon  migration  in  scattering  media  and  compute  both  absorption  and  scattering 
coefficients.  For  this  study,  I  used  a  new  dual  wavelength,  homodyne  frequency-domain 
near-infrared  spectroscopy  (NIRS)  system  to  measure  changes  in  tumor  vascular 
oxyhemoglobin  concentration  [Hb02]  and  total  hemoglobin  concentration  [Hb]T0tai-  The 
theoretical  formulation  of  the  system  is  based  on  a  modified  Beer-Lambert's  law. 
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•  In-Phase  and  Quadrature  Phase  Detection  System 

The  dual  wavelength,  homodyne  NIRS  system  (wavelengths  758  nm  and  782  nm) 
is  based  on  an  In-phase  and  Quadrature-phase  chip  (I&Q  chip).  These  wavelengths  were 
chosen  because  they  not  only  allow  the  calculation  of  [Hb02]  and  [Hb]  iotai,  but  also  fall 
into  the  range  of  wavelengths  compatible  with  the  low  cost  photo  multiplier  tube  (PMT). 
The  system  uses  only  one  RF  source  to  modulate  the  light  emitted  by  two  laser  diodes 
and  lets  the  I&Q  chip  to  determine  amplitude  and  phase  changes  of  light  attenuated  by 
the  sample  [74], 
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Figure  22.  A  schematic  diagram  of  the  NIRS  I&Q  system. 


Figure  22  shows  a  schematic  diagram  of  the  NIRS  system.  An  RF  source 
modulates  the  light  from  two  laser  diodes  (LD1  and  LD2)  at  140  MHz  through  a  time¬ 
sharing  system.  The  light  passes  through  fiber  optic  cables,  is  transmitted  through  the 
tumor  tissue,  and  collected  by  a  second  fiber  bundle.  The  light  is  then  amplified  by  a 
photo  multiplier  tube  (PMT)  and  filtered  by  a  bandpass  (BP)  filter  to  pick  up  only  the 
signal  being  modulated  at  140  MHz.  The  signal  is  demodulated  into  I  and  Q  components 
by  the  I&Q  chip  and  filtered  again  by  two  lowpass  filters  (LPFs)  to  select  the  DC 
components.  These  final  DC  signals  are  digitized  by  a  12-bit  analog  to  digital  converter 
(ADC)  and  stored  in  a  computer.  Light  amplitude  and  phase  changes  caused  by  the  tumor 
attenuation  are  used  to  compute  changes  in  tumor  vascular  oxyhemoglobin  concentration 
[Hb02]  and  total  hemoglobin  concentration  [Hbjiotai-  The  NIRS  system  is  interfaced  to  a 
laptop  via  a  National  Instruments  data  acquisition  card  (DAQCard-1200)  that  is  fully 
software-configurable.  A  National  Instruments  DAQ  driver,  NI-DAQ,  is  used  to  control 
the  operations  of  the  card.  The  graphical  programming  language  LabViewrv/  is  used  to 
create  a  graphical  user  interface;  (GUI)  for  controlling  data  acquisition,  processing,  and 
display  [75],  The  raw  data  are  saved  in  a  text  format  and  can  be  further  processed  using  a 
variety  of  digital  signal  processing  (DSP)  techniques. 
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•  Algorithms  for  Computing  Tumor  Vascular  [HbOi]  and  [HbjTotai 

To  obtain  the  absorption  coefficients  of  deoxyhemoglobin  and  oxyhemoglobin,  it 
is  assumed  that  background  absorbance  is  negligible  and  deoxyhemoglobin  and 
oxyhemoglobin  are  the  predominant  light  absorbing  molecules  in  the  tumor  tissue. 
Therefore,  the  absorption  coefficient  of  the  tumor  can  be  approximated  as  the  product  of 
the  extinction  coefficients  for  deoxyhemoglobin  and  oxyhemoglobin  and  their  respective 
concentrations: 

A  =4,[Hb]  +  4b0l[Hb02]  (22) 

Equation  (22)  has  two  unknowns:  [Hb]  and  [Hb02]  and  we  can  not  solve  for  them  with 
only  one  equation.  If  we  use  N1R  light  at  two  different  wavelengths,  then  we  have  two 
equations  and  [Hb]  and  [Hb02]  can  be  obtained  by  solving  following  linear  systems  of 
equations 


^758  =4?[Hb]  +  ^o2[Hb02] 

A'I82=4f[Hb]+4bo2[Hb02] 


(23) 


where  juJ52,  and  //f,782  are  the  absorption  coefficients,  and  s1^  the  extinction 

coefficients  for  deoxyhemoglobin,  and  £^bo2  316  the  extinction  coefficients  for 

oxyhemoglobin  at  the  wavelengths  758  nm  and  782  nm,  respectively,  and  [Hb]  and 
[Hb02]  are  the  deoxyhemoglobin  and  oxyhemoglobin  concentrations,  respectively.  Since 

^Hb8  >  £ub  >  £Hbo2  >  £mo2  are  physical  constants,  changes  in  [Hb02]  and  [Hb]  in  tumor 

tissue  vasculature  cause  changes  in  jua  and  /Ja  according  to  Equations  (23).  Thus,  by 
measuring  changes  in  jua  ancl  jua  ,  we  can  determine  changes  in  [Hb02]  and  [Hb], 
Combining  Equation  (23)  with  modified  Beer-Lambert’s  law  yields  equations  for 
computing  changes  in  [Hb]  and  [Hb02][76], 
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(25) 


where  letters  P  and  B  stand  for  perturbation  and  baseline,  respectively,  AB  is  the  baseline 
light  amplitude,  Ap  is  the  light  amplitude  under  the  physiological  perturbation,  L  is  the 
distance  between  the  source  and  detector,  and  A[]  represents  a  change  in  concentration. 
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The  quantity  (AB  /  Ap)  can  be  obtained  by  direct  measurement.  Given  s™  =  0.359 

cm'mM'1,  ^Hbo,  =0.1496  cm  'mM'1,  £-^2  =  0.265  cm^mM'1,  and  £^q2  =0.178  cm" 

1mM"1  [77],  the  final  expressions  for  computing  changes  in  [Hb]  (mM),  [Hb02]  (mM), 
and  [Hbjiotai  (mM)  are: 
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(28) 


where  Equation  (28)  represents  a  change  in  total  hemoglobin  concentration  in  the  tumor, 
a  physiological  parameter  that  is  linearly  proportional  to  a  change  in  blood  volume  in  the 
tumor. 

•  Assessment  of  Temporal  Characteristics  of  Oxygen  Dynamic  Response 

As  a  first  step  in  trying  to  quantitatively  characterize  the  temporal  response  in 
[HbC>2],  a  simple  mono-exponential  model  was  employed  to  compute  the  time  constant 
and  the  amplitude  of  the  response: 


y  =  A  -|l-exp[-A_Jol 

where  y  was  measured  [Hb02]  (mM),  A  was  a  constant,  representing  the  amplitude  of  the 
response,  t  (min)  was  time,  t0  (min)  was  the  beginning  time  of  the  exponential  response, 
and  r  (min)  was  the  time  constant  of  the  exponential  response. 

This  mono-exponential  model  was  based  on  an  implicit  assumption  that  a  tumor  is 
composed  of  a  homogenous  tissue  with  a  single  absorption  coefficient  ju  of  light  or  a 
single  uniform  compartment.  Therefore,  the  whole  tumor  should  respond  to  carbogen 
intervention  in  a  single  mode  and  this  dynamic  process  can  be  completely  characterized 
by  only  two  parameters:  amplitude  A  and  time  constant  r.  This  simple  mono-exponential 
model,  however,  did  not  take  into  account  the  important  fact  that  most  tumors  exhibit 
substantially  spatial  heterogeneity  in  both  tissue  density  and  blood  vessel  distribution. 
Therefore,  a  more  sophisticated  and  realistic  mathematical  model,  the  bi-exponential 
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model,  was  also  used.  The  model  assumed  that  a  tumor  consists  of  two  distinct 
homogenous  tissue  regions,  a  well-perfused  region  and  a  poorly-perfused  region.  The  two 
tissue  regions  can  be  characterized  optically  by  two  different  absorption  coefficients, 
/j\  and  /J-2  and  kinetically  by  two  different  time  constants,  T\  and  r2.  The  bi-exponential 
model  is 


y  =  A 


•  \  1  -  exp[- 


+  A2 


•si-  exp[- 


(^o)]l 
*2  j 


(30) 


where  t0  (min)  was  the  beginning  time  of  the  fast  and  slow  exponential  response,  and  xi 
and  x2  (min)  were  the  fast  and  slow  time  constants,  respectively.  The  parameters  in  the 
two  models  were  obtained  by  fitting  the  data  to  Equations  (29)  and  (30),  respectively, 
using  KaleidaGraph. 

•  Investigation  of  Tumor  Vascular  [HbOj]  and  [Hbjxotai  by  NIRS 


Figure  23.  A  schematic  diagram  of  transmittance  optical  mode 

Animal  preparation  was  similar  to  19F  MR  EP1  studies.  The  experiments  were 
performed  in  transmittance  optical  mode.  In  this  configuration,  the  light  source  and 
detector  were  positioned  on  opposite  sides  of  the  tumor  and  stabilized  by  a  mechanical 
mechanism.  The  light  passed  through  and  was  attenuated  by  the  tumor.  The  transmitted 
light  was  then  picked  up  by  the  detector.  In  this  way,  the  NIR  spectroscopy  measured 
A[Hb02]  and  A[Hb]Totai  in  both  peripheral  and  central  regions  of  the  tumor,  yielding  an 
average  A[Hb02]  and  A[Hb]Totai  (Figure  23).  It  was  important  to  position  the  light  source 
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and  detector  in  such  a  way  that  they  made  very  good  optical  contact  with  the  tumor  but 
did  not  exert  excessive  force  on  it.  This  was  because  excessive  pressure  on  the  tumor 
could  hinder  the  tumor  blood  flow  and  thus,  affect  A[Hb02]  and  A[Hb]Totai-  To  avoid 
instrument  drift,  the  NIRS  system  was  first  wanned  up  for  about  30  minutes.  Then,  5  ■*- 
10  minute  test  data  were  acquired  to  check  the  stability  of  the  system  using  a  tissue 
equivalent  phantom  with  stable  optical  properties  [78].  Once  the  system  reached  a  steady 
state,  a  series  of  NIRS  experiments  were  performed  according  to  the  respiratory 
challenge  paradigms  described  earlier. 


Figure  24.  Result  of  a  typical  routine  NIRS  system  drift  test  for 
A[Hb]totai  which  was  obtained  by  using  Equation  (28). 


Table  14.  Summery  of  five  system  drift 
tests  obtained  using  a  tissue  phantom 


NO 

SD  of[Hb02](mM) 

SD  of  [lib]  jotai  (mM) 

1 

0.0053 

0.0037 

2 

0.0043 

0.0025 

3 

0.0056 

0.0034 

4 

0.0066 

0.0032 

5 

0.0046 

0.0023 

Figure  24  shows  the  results  of  a  typical  routine  system  drift  test  for  A[Hb]Totai  over 
a  period  of  33  minutes.  The  test  was  performed  using  a  tissue  equivalent  phantom  with  a 
source-to-detector  separation  of  5.0  cm.  The  result  was  obtained  using  Equation  (28).  For 
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this  particular  drift  test,  the  standard  deviation  (SD)  for  A[Hb]Totai  was  0.0034  (mM), 
indicating  that  the  NIRS  system  had  a  superior  stability.  Table  14  is  a  summery  of  five 
system  drift  tests  obtained  under  identical  experimental  conditions,  but  at  different  dates. 
As  can  be  seen  from  Table  14,  the 


dtauuoiu  ugviauuas  iui  cm  wclc  Oi  Hit'  aaiiiv 


order,  suggesting  the  NIRS  system  was  very  stable  in  terms  of  variability.  In  addition,  T 
also  used  a  chunk  of  fresh  red  port  meat  (3x3x4  cm3)  to  further  evaluate  the  system 
stability.  In  this  case,  the  source-to-detector  separation  was  3.0  cm  and  die  drift  tests  'were 
repeated  for  four  times,  with  a  time  interval  of  one  hour  between  the  tests.  Table  15 
summarizes  the  results. 


Table  15.  Summery  of  four  system  drift  tests  obtained 
using  a  chunk  of  fresh  port  meat  (3x3x4  cm3). 


NO 

SD  of  [Hb02]  (mM) 

SD  of  [Hbjxotai  (mM) 

1 

0.0047 

0.0018 

2 

0.0048 

0.0019 

3 

0.0053 

0.0023 

4 

0.0043 

0.0016 
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Figure  25.  Time  course  profile  of  A[Hb02]  and  A[Hb]Totai  in  response  to  carbogen 
intervention  for  a  representative  13762NF  breast  tumor  (14.8  cm3). 


Figure  25  shows  the  time  course  impacts  of  inhaled  gases  on  changes  in  tumor 
vascular  A[Hb02]  and  A[Hb]x0tai  for  a  representative  13762NF  breast  tumor  (14.8  cm3). 
Repeated  carbogen  interventions  were  performed  sequentially  to  evaluate  the 
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reproducibility  of  the  time  course  profile  of  the  tumor.  The  data  were  acquired  in 
transmittance  mode  with  a  source-to-detector  distance  of  2.0  cm.  To  show  the  quality  of 
the  raw  data,  averaging  and  filtering  were  not  applied  to  the  curve.  The  vertical  dotted 
lines  mark  the  beginning  of  each  gas  switch.  The  measurement  uncertainties  in  A[Hb02] 
and  A[Hb]Totai  were  estimated  using  baseline  data  based  on  error  propagation  theory  and 
were  only  labeled  at  representative  locations.  The  initial  baseline  A[Hb02]  value  was 
about  0  mM  when  the  rat  was  breathing  air  over  a  period  of  10  minutes.  Immediately 
following  a  gas  switch  from  air  to  carbogen,  A[Hb02]  rose  momentarily  and  significantly 
( p  <  0.0001)  from  0  mM  to  about  0.2  mM  within  the  first  70  second  and  then  increased 
further  at  a  slower,  but  still  significant  rate  (p  <  0.0001)  to  about  0.25  mM  for  the  next  20 
minutes  until  an  apparent  saturation  was  reached.  After  the  gas  was  switched  back  to  air, 
A[Hb02]  did  not  decrease  until  about  30-40  seconds  later  when  a  sudden  and  significant 
drop  occurred,  followed  by  a  gradual  return  to  the  baseline.  Note  that  there  was  a  small 
overshoot  in  A[Hb02]  as  it  was  returning  to  the  baseline  and  this  overshoot  also  appeared 
in  the  second  carbogen  intervention.  It  was  found  that,  for  some  tumors,  this  final 
baseline  value  was  little  bit  higher  than  the  initial  baseline  value.  Detailed  analysis  also 
revealed  that  A[Hb]T0tai  showed  a  similar  time  course  profile,  but  at  a  much  less  dramatic 
and  significant  level  in  terms  of  both  time  constant  and  amplitude.  It  is  important  to  point 
out  that  similar  time  course  response  patterns  were  observed  in  all  13762NF  breast 
tumors,  though  the  time  constants  and  amplitudes  of  A[Hb02]  and  A[Hb]iotai  varied  with 
tumor  size. 


Figure  26.  Time  course  profile  of  the  systemic  arterial  sa02  obtained 
with  a  pulse  oximeter  in  response  to  carbogen  intervention.  The  error 
bars  represent  standard  deviation  estimated  using  the  baseline  data 


The  systemic  arterial  hemoglobin  oxygen  saturation  sa02  of  the  rat  was  also 
monitored  using  a  pulse  oximeter  placed  on  the  hind  foot.  This  was  designed  to 
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investigate  sa 02  response  to  carbogen  intervention  and  determine  the  relationship 
between  sa 02  and  A[Hb02].  Figure  26  shows  sa 02  time  course  profile  of  the  same  rat  in 
response  to  carbogen  intervention.  The  data  were  acquired  simultaneously  with  A[Hb02] 
and  A[Hb]Totai  data  and  smoothed  during  postprocessing.  In  comparison  to  A[Hb02] 
response,  sa02  rose  almost  immediately  from  76%  to  97%  in  about  20  seconds  following 
the  onset  of  carbogen  breathing,  indicating  a  much  shorter  time  constant  than  that  of 
A[Hb02].  sa02  remained  at  97%  during  the  course  of  carbogen  administration  over  the 
next  20  minutes  until  the  cessation  of  carbogen  breathing.  Then  sa02  dropped  very 
rapidly  to  its  baseline  value  of  76%.  Twenty  minutes  later,  the  second  cycle  of  carbogen 
intervention  was  initiated  and  sa 02  showed  an  almost  identical  response  profile. 


Figure  27.  The  best  bi-exponential  fit  gave:  0.031  +  0.177{l-exp[-(/-9.6)/0.293]} 
+0.23{l-exp[-(f-  11.43)/3.985]},  with  r  -  0.991,  whereas  the  best  mono¬ 
exponential  fit  gave:  0.033  +  0.179(1  -  exp[-(t  -  9.61  )/0.53]},  with  r  —  0.971. 


To  quantitatively  characterize  the  dynamic  features  of  A[Hb02],  we  computed  the 
time  constants  of  A[Hb02]  response  using  two  different  mathematical  models:  a  mono¬ 
exponential  model  and  a  bi-exponential  model.  The  models  were  fit  to  the  rising  portion 
of  the  raw  data  corresponding  to  the  first  carbogen  intervention.  Figure  27  shows  only  the 
bi-exponential  curve  fit.  The  bi-exponential  model  gave  two  time  constants:  Xi  = 
0.293  ±  0.023  (min)  and  x2  =  3 .985  ±  0.056  (min)  (r  =  0.99 1),  representing  the  two  phases 
of  a  bi-modal  response:  the  fast  phase  and  the  slow  phase.  In  terms  of  tumor  perfusion 
based  on  the  tumor  model  depicted  in  Figure  23,  the  former  could  be  a  physiological 
indicator  of  tumor  blood  perfUsion  in  well-oxygenated  regions  and  the  latter  in  poorly- 
oxygenated  regions.  The  mono-exponential  model  yielded  only  one  time  constant:  x  = 
0.530±0.015  (min)  (r  =  0.971),  a  possible  physiological  indicator  of  tumor  blood 
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perfusion  for  a  homogeneous  tumor  model.  The  bi-exponential  model  gave  a  better  curve 
fit  than  the  mono-exponential  model,  as  is  manifested  by  their  respective  r  values.  This 
was  true  for  all  tumors,  especially  for  larger  tumors. 


Figure  28.  Time  course  profile  of  AfHbCy  and  A[Hb]xotai  in  response 
to  carbogen  intervention  for  a  second  breast  tumor  (9.9  cm3). 


Figure  28  shows  the  time  course  profile  of  A[Hb02]  and  A[Hb]xotai  in  response  to 
carbogen  intervention  for  a  second  breast  tumor  (9.9  cm3).  The  data  were  also  acquired  in 
transmittance  mode  with  a  source-to-detector  distance  of  1.9  cm.  Again,  we  saw  a  similar 
response  pattern  here.  After  the  gas  was  switched  from  air  to  carbogen,  A[Hb02]  rose 
sharply  and  significantly  (p  <  0,0001)  from  0  mM  to  about  0.7  mM  within  the  first  2,5 
minutes.  This  was  a  three-fold  increase  in  amplitude  as  compared  to  the  first  tumor, 
suggesting  that  this  is  a  better-perfused  tumor.  A[Hb02]  continued  to  increase  further  at  a 
slower  rate  for  the  rest  duration  of  carbogen  administration.  After  the  gas  was  switched 
back  to  air,  A[Hb02]  decreased  monotonically  and  exponentially  back  to  the  baseline. 
Here,  we  did  not  observe  the  small  overshoot  as  seen  in  the  first  tumor  during  the  course 
of  air  breathing.  Interestingly,  in  this  case,  A[Hb]T(>tai  responded  differently  to  carbogen 
intervention  from  the  first  tumor.  It  did  not  exhibit  a  regular  response  pattern.  Initially, 
A[Hb]xotai  increased  sharply  in  response  to  carbogen  intervention,  but  showed  very  little 
changes  for  the  subsequent  gas  switches.  Again,  the  dynamic  response  in  A[Hb02]  was 
analyzed  using  the  two  exponential  models.  The  bi-exponential  model  gave:  = 

0.853  ±0.0079  (min)  and  x2  =  7.555  ±0.089  (min)  (r  =  0.991),  whereas  the  mono¬ 
exponential  model  gave:  x  =  1.655  ±0.022  (min)  (r  =  0.938).  As  in  the  first  case,  the  bi¬ 
exponential  model  provided  a  better  curve  fit. 
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We  found  that  there  was  a  strong  linear  relationship  (r  =  0.843)  between  AHb]Totai 
and  maximum  A[Hb02]  for  a  group  of  38  breast  tumors  as  shown  in  Figure  29.  This 
suggests  that  an  increase  in  A[Hb02]  was  caused  in  part  by  an  increase  in  A[Hb]T0tai ,  i-e-, 
an  increase  in  total  blood  volume  in  tumors.  This  means  that  breathing  carbogen  can 
increase  tumor  blood  flow  and,  thus,  improve  tumor  oxygenation. 


Figure  29.  Relationship  between  A[Hb]r0tai  and  AfHbCh] 
for  a  group  of  38  breast  tumors. 


Figure  30.  Relationship  between  fast  time  constant  tland  tumor  volume. 


We  also  systematically  investigated  the  relationships  between  time  constants  and 
tumor  volume.  Figure  30  shows  the  correlation  between  fast  time  constant  t1  and  tumor 
volume  and  Figure  3 1  shows  the  correlation  between  slow  time  constant  x2  and  tumor 


52 


volume  for  a  group  of  31  tumors.  As  can  be  seen  from  Figures  30  and  31,  there  were  no 
clear  relationships  between  time  constants  and  tumor  volume. 
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Figure  31.  Relationship  between  slow  time  constant  t2  and  tumor  volume. 


To  Investigate  Tumor  Physiology  in  Response  to  Irradiation 

One  of  the  reasons  for  the  birth  of  modem  fractionated  radiotherapy  was  to  exploit 
tumor  reoxygenation  that  had  been  found  to  occur  following  each  irradiation.  Since  early 
1980s,  many  new  fractionation  protocols  have  been  designed,  including 
hyperfractionation  and  accelerated  fractionation.  Currently,  however,  the  timing  of 
successive  doses  is  not  optimized  for  individual  patients/tumors,  but  rather  is  based  on 
experience  with  cell  culture,  animal  tumor  model  systems,  and  clinical  trials.  The 
possibility  of  monitoring  changes  in  tumor  oxygenation  following  radiotherapy  would  of 
great  significance  in  predicting  the  optimal  time  for  delivering  sequential  radiation  doses 
during  a  course  of  fractionated  radiotherapy. 

•  Irradiation 

Mammary  adenocarcinomas  13762NF  were  implanted  in  pedicles  as  described 
earlier.  When  tumors  reached  about  1  cm  in  diameter,  they  were  irradiated  using  a  Varian 
CLIN  AC  4-100  (4  MeV  photon  beam)  in  the  Department  of  Radiation  Oncology  outside 
normal  business  hours.  Rats  were  anesthetized  with  200  pi  ketamine  hydrochloride  i.p. 
and  maintained  under  general  gaseous  anesthesia  with  air  (1.0  dm3/min)  and  1.0% 
isoflurane  for  1  hour  prior  to  irradiation  to  stabilize  tumor  oxygenation.  Tissue  equivalent 
bolus  material  was  placed  around  the  tumors  to  ensure  uniformity  of  dose  to  the  tumors. 
A  treatment  plan  was  computed  based  on  a  source-to-axis  distance  (SAD)  of  100  cm  and 
a  field  size  of  5  x  5  cm.  Tumors  were  irradiated  with  a  4  MeV  photon  beam  from  both 
sides  with  10  Gy  delivered  from  each  side  at  a  rate  of  2  Gy /min,  giving  a  total  dose  of  20 
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Gy.  Prior  to  irradiation,  baseline  tumor  oximetry  was  performed  and  then,  the 
measurements  were  repeated  again  following  irradiation  in  order  to  examine  the  extent 
and  time  course  of  tumor  hypoxiation  and  reoxygenation.  Tumor  size  was  measured  once 
every  two  day  following  irradiation. 
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Figure  32.  Tumor  volume  as  a  function  of  time 
following  a  single  dose  irradiation  (20  Gy). 


Figure  32  shows  the  growth  curve  of  a  representative  mammary  adenocarcinoma 
13762NF  before  and  after  a  single  dose  irradiation  (20  Gy).  Irradiation  caused  a 
significant  decrease  in  tumor  volume  (p  <  0.007).  Changes  in  tumor  vascular 
oxyhemoglobin  concentration  A[HbC>2]  were  measured  at  2  hours,  24  hours,  48  hours, 
and  96  hours  following  radiotherapy  and  the  results,  obtained  by  fitting  the  raw  data  to 
the  bi-exponential  model,  are  listed  in  Table  16. 


Table  16.  Results  of  tumor  oxygenation  following  radiotherapy 


Vd(cnf) 

Al(niV) 

Tl(nin) 

A2(itM) 

r2(nin) 

A1/A2 

(A1/A2V(t1/t2) 

r 

Ftejrrad&ion 

3. 85  err? 

0.471 

0.241 

0.161 

5.357 

2932 

65184 

0.947 

2HsFtet-liTada 

285  err? 

1.004 

0.256 

0.304 

61.305 

3.365 

806.494 

Q951 

24HsRost-lrrada 

245  cn? 

Q140 

0.071 

0.348 

4157 

0.401 

23.358 

0.957 

48RsFte-lrrada 

236  err? 

0.002 

0113 

0.080 

4.185 

0.309 

14.817 

0987 

96HsFbst-li7ada 

1.99  cn? 

0.117 

4.871 

0006 

4.886 

1.355 

1.360 
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7)  Dissertation 

The  studies  described  in  this  report  constitute  a  part  of  my  dissertation,  which  will 
be  completed  by  the  end  of  the  year. 
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4.  KEY  RESEARCH  ACCOMPLISHMENTS 


•  The  pedicle  tumor  model  w;is  proven  to  be  ideally  suited  for  in  vivo  NIRS  and  MRI 
studies,  therapy,  and  manipulation. 

•  The  pedicle  tumor  model  allowed  accurate  measurement  of  tumor  size  and  had  no 
significant  difference  from  the  traditional  subcutaneous  site  in  the  thigh  in  terms  of 
growth. 

•  Mammary  adenocarcinomas  13762NF  were  found  to  have  a  viable,  well-oxygenated 
peripheral  region,  and  a  necrotic,  poorly-oxygenated  central  region. 

•  Mammary  adenocarcinomas  13762NF  were  found  to  have  a  silent  phase  ranging  from 
2  to  3  weeks  and  a  mean  volume  doubling  time  (VDT)  of  4  days. 

•  The  first  order,  autonomous  differential  equation  was  proven  to  be  a  good  tumor 
growth  model  at  least  for  a  short  period  of  growth  for  mammary  adenocarcinomas 
13762NF. 

•  For  the  birdcage  resonator,  laboratory  bench  testing  indicated  that  the  unloaded  Qs  for 
both  proton  and  fluorine  resonant  modes  were  similar  and  relatively  low.  The  loaded 
Q  dropped  by  about  36%  for  fluorine  and  39%  for  proton.  Phantom  imaging  found 
that  the  Bx  fields  were  somewhat  heterogeneous  across  usable  volume  of  the  resonator 
for  both  resonant  modes. 

•  The  slotted  tube  resonator  had  the  capability  of  being  continuously  tuned  in  the 
frequency  range  of  150  MHz  ~  220  MHz.  The  resonator  had  high  Q  values  and  short 
90°  pulse  widths  for  both  proton  and  fluorine.  In  addition,  the  field  homogeneity 
was  found  to  be  excellent  for  proton  resonance  and  reasonably  good  for  fluorine 
resonance. 

•  Both  resonant  frequencies  of  the  slotted  tube  resonator  were  stable  and  immune  to 
external  electromagnetic  interference. 

•  The  copper  shield  not  only  lowered  the  Q  value  of  the  resonator,  but  also  increased 
the  RF  pulse  width. 

•  For  the  slotted  tube  resonator,  NMR  samples  had  a  much  stronger  effect  on  the  19F 
resonance  than  on  the  'H  resonance. 


56 


•  The  HFB  signal  intensity  was  found  to  decay  exponentially  with  a  typical  biological 
half-life  ranging  from  Ti/2  =  700  to  1200  min,  which,  we  believe,  would  provide  an 
indication  of  relative  tumor  blood  flow  (TBF). 

•  The  global  and  regional  clearance  and  redistribution  of  HFB  within  the  tumors  did  not 
interfere  with  19F  MR  EPI  oximetry. 

•  Tumor  voxels  with  high  baseline  p02  had  significantly  different  response 
characteristics  from  those  with  initially  low  p02,  with  voxels  of  high  baseline  p02 
showing  significant  changes  in  p02  while  voxles  of  low  baseline  p02  showing  small 
changes. 

•  Time  constants  (x)  of  well-oxygenated  voxels  (10  ~  20  min)  were  much  shorter  than 
those  of  hypoxic  voxels  (>  50  min).  The  global  p02  time  constant  (60  ~  80  min)  was 
much  longer  than  the  blood  hemoglobin  saturation  (s02)  time  constant  (10-20  min). 

•  NIR  spectroscopy  showed  significant  changes  in  tumor  vascular  oxygenation  (s02) 
accompanying  respiratory  interventions.  19F  MR-EPI  showed  significant  changes  in 
tumor  tissue  p02,  with  considerable  regional  heterogeneity  in  both  absolute  values 
and  rate  of  change  accompanying  interventions.  Changes  in  tumor  vascular  s02 
preceded  tumor  tissue  p02,  particularly  for  smaller  tumors. 

•  Strong  correlation  existed  between  the  maximum  p02  value  attained  during  the  course 
of  an  experiment  and  mean  baseline  p02  and  between  mean  p02  and  mean  baseline 
p02. 

•  19F  EPI  oximetry  of  HFB  was  proven  to  be  a  useful  technique  for  measuring  tumor 
oxygenation. 

•  p02  and  the  distribution  of  p02  in  breast  adenocarcinomas  13762  NF  changed  with 
tumor  growth  and  there  existed  heterogeneity  in  p02  distribution. 

•  Tumor  oxygenation  could  be  manipulated  by  inhaling  different  gases. 

•  We  found  that  breathing  elevated  F02  had  a  significant  effect  on  arterial  sa02,  tumor 
vascular  [H02],  and  tumor  tissue  p02.  sa02  had  the  fastest  response,  followed  by 
[H02],  and  p02. 

•  NIR  spectroscopy  is  completely  non-invasive,  inexpensive,  portable,  and  amenable  to 
real-time  measurements. 
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•  We  found  that  there  was  a  slrong  linear  relationship  between  AHbJjotai  and  maximum 
A[HbC>2]  following  carbogen  intervention. 

•  We  found  that  there  was  no  clear  relationship  between  time  constants,  as  determined 
by  NIRS,  and  tumor  volume. 

•  Irradiation  caused  a  significant  decrease  in  tumor  volume  (p  <  0.007). 

•  Tumor  vascular  oxyhemoglobin  concentration  A[HbC>2]  changed  following 
radiotherapy. 
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(Presented  at  the  Era  of  Hope  Meeting  sponsored  by  the  Department  of  Defense,  Atlanta,  GA, 
June  8-1 1,2000). 

9)  "FREDOM  (Fluorine  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen 
Mapping):  a  contextual  review" 

Ralph  P  Mason.  Dawen  Zhao  Anca  Constantinescu,  Yulin  Song,  Lan  Jiang  and  Eric  W.  Hahn 
(Presentated  at  ISMRM  Workshop  on  MR  in  Experimental  and  Clinical  Cancer  Research  in  the 
New  Millennium). 

10)  "Diverse  Approaches  to  Monitoring  Oxygen  Dynamics  in  Rat  Breast  and  Prostate  Tumors" 

Dawen.  Zhao,  Yulin  Song,  Han  Liu,  Amca.  Constantinescu.  Eric  W.  Hahn,  and  Ralph.  P.  Mason, 
(Presented  at  the  Conference  on  Chemical  Modifiers  of  Cancer  Treatment,  Banff,  Canada,  Oct. 
2000). 


— DEGREE  OBTAINED  THAT  IS  SUPPORTED  BY  THIS  AWARD 

Ph.D.  in  Biomedical  Engineering 

—  EMPLOYMENT  OR  RESEARCH  OPPORTUNITIES  RECEIVED  ON 
TRAINING  SUPPORTED  BY  THIS  AWARD 

Post-doctoral  fellow  in  the  Department  of  Radiation  Oncology,  School  of  Medicine, 
Stanford  University  (will  start  on  March  1,  2001) 
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6  CONCLUSIONS 


Regional  tumor  tissue  p02,  vascular  [Hb02]  or  s02,  and  [Hb]Totai  are  critical 
physiological  parameters  in  radiotherapy  and  some  forms  of  chemotherapy.  The 
capability  to  measure  and  manipulate  them  will  provide  insight  into  progressive 
physiological  changes  in  a  tumor  accompanying  interventions  and  enhance  therapeutic 
outcomes.  NIRS  has  the  advantages  of  being  entirely  noninvasive,  inexpensive,  portable, 
and  real-time.  But  the  19F  MR  EPI  approach  clearly  reveals  detailed  oxygenation 
heterogeneity.  The  correlation  of  NIRS  and  19F  MR  EPI  technologies  will  help  us 
understand  issues  of  oxygen  transport,  perfusion,  and  consumption.  I  believe  that 
synergistic  application  of  multiple  approaches  to  tumor  oxygenation  can  lead  to 
optimized  tumor  therapy. 
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Noninvasive  investigation  of  blood  oxygenation 
dynamics  of  tumors  by  near-infrared  spectroscopy 


Hanli  Liu,  Yulin  Song,  Katherine  L.  Worden,  Xin  Jiang,  Anca  Constantinescu,  and 
Ralph  P.  Mason 


The  measurement  of  dynamic  changes  in  the  blood  oxygenation  of  tumor  vasculature  could  be  valuable 
for  tumor  prognosis  and  optimizing  tumor  treatment  plans.  In  this  study  we  employed  near-infrared 
spectroscopy  (NIRS)  to  measure  changes  in  the  total  hemoglobin  concentration  together  with  the  degree 
of  hemoglobin  oxygenation  in  the  vascular  bed  of  breast  and  prostate  tumors  implanted  in  rats.  Mea¬ 
surements  were  made  while  inhaled  gas  was  alternated  between  33%  oxygen  and  carbogen  (95%  02,  5% 
C02).  Significant  dynamic  changes  in  tumor  oxygenation  were  observed  to  accompany  respiratory 
challenge,  and  these  changes  could  be  modeled  with  two  exponential  components,  yielding  two  time 
constants.  Following  the  Fick  principle,  we  derived  a  simplified  model  to  relate  the  time  constants  to 
tumor  blood-perfusion  rates.  This  study  demonstrates  that  the  NIRS  technology  can  provide  an  effi¬ 
cient,  real-time,  noninvasive  means  of  monitoring  the  vascular  oxygenation  dynamics  of  tumors  and 
facilitate  investigations  of  tumor  vascular  perfusion.  This  may  have  prognostic  value  and  promises 
insight  into  tumor  vascular  development.  ©  2000  Optical  Society  of  America 
OCIS  codes:  170.1470,  170.3660,  170.4580,  170.5280,  290.1990,  290.7050. 


1 .  introduction 

The  presence  and  the  significance  of  tumor  hypoxia 
have  been  recognized  since  the  1950's.  There  is  in¬ 
creasing  evidence  that  tumor  oxygenation  is  clini¬ 
cally  important  in  predicting  tumor  response  to 
radiation,  tumor  response  to  chemotherapy,  overall 
prognosis,  or  all  three.  Hypoxic  cells  in  vitro  and  in 
animal  tumors  in  vivo  are  documented  to  be  3  times 
more  resistant  to  radiation-induced  killing  compared 
with  aerobic  cells.1  Recent  studies  show  that  hy¬ 
poxia  may  have  a  profound  impact  on  malignant  pro¬ 
gression  and  on  responsiveness  to  therapy.23 
Numerous  studies  on  tumor  oxygen  tension  (p02) 
measurements  have  been  conducted  in  recent  years 
by  use  of  a  variety  of  methods,  such  as  microelec¬ 
trodes,2  phosphors,4  electron  paramagnetic  reso¬ 
nance,5  or  magnetic  resonance  imaging6  (MRI). 
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Comparing  needle-based,  oxygen-sensitive  electrodes 
or  electron  paramagnetic  resonance  and  MRI  for 
measuring  p02  shows  that  the  latter  two  offer  the 
advantage  of  facilitating  multiple  repeated  measure¬ 
ments  to  map  p02  noninvasively.  However,  mag¬ 
nets  are  large,  and  the  methods  are  not  readily 
portable.  A  versatile  method  for  monitoring  intra¬ 
tumor  oxygenation  rapidly  and  noninvasively  is 
therefore  very  desirable  for  tumor  prognosis  and  tu¬ 
mor  treatment  planning. 

In  the  near-infrared  (NIR)  region  (700-900  nm) 
the  major  chromophores  in  tissue  are  oxygenated  he¬ 
moglobin  and  deoxygenated  hemoglobin,  which  differ 
in  their  light  absorption.  Measurements  of  the  ab¬ 
sorption  of  light  travelling  through  the  tissue  under 
study  allow  us  to  evaluate  or  quantify  blood  oxygen¬ 
ation,  such  as  the  concentrations  of  oxygenated  he¬ 
moglobin  (Hb02),  and  deoxygenated  hemoglobin  (Hb) 
and  the  hemoglobin  saturation  S02.  In  the  past  de¬ 
cade,  three  forms  of  NIR  spectroscopy  (NIRS)  that 
use  pulsed-laser  light  in  the  time  domain,  amplitude- 
modulated  laser  light  in  the  frequency  domain,  and 
cw  light  in  a  dc  form  were  developed  for  blood  oxy¬ 
genation  quantification  in  tissue.7  Significant  in¬ 
vestigations  in  both  laboratory  and  clinical  settings 
by  use  of  NIRS  were  conducted  for  noninvasive,  quan¬ 
titative  measurements  and  imaging  of  cerebral 
oxygenation8-12  and  blood  oxygenation  of  exercised 
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muscle13-17  in  vivo.  Although  NIR  techniques  were 
used  extensively  in  conjunction  with  cryospectropho- 
tometry  to  investigate  tumor  blood-vessel  oxygen¬ 
ation  in  biopsies,18  only  a  few  reports19-22  were 
published  on  using  the  NIR  techniques  for  monitor¬ 
ing  tumor  oxygenation  in  vivo.  In  principle,  the  the¬ 
oretical  model,  i.e.,  the  diffusion  approximation  to  the 
photon  transport  theory,  works  well  for  only  large 
and  homogeneous  media.23’24  Accurate  quantifica¬ 
tion  of  tumor  oxygenation  by  use  of  the  NIR  approach 
is  limited  because  of  the  considerable  heterogeneity 
and  the  finite  sizes  of  tumors. 

It  is  understood  and  documented25  that  the  NIR 
technique  used  for  blood  oxygenation  monitoring  is 
sensitive  to  vascular  absorption  in  the  measured  or¬ 
gan.  The  NIR  method  is  not  limited  to  measure¬ 
ments  of  blood  oxygenation  in  arteries  (c.f.,  pulse 
oximetry)  or  in  veins  but  interrogates  blood  in  the 
entire  vascular  compartment,  including  capillaries, 
arterioles,  and  venules,  i.e.,  the  vascular  bed.  A  va¬ 
riety  of  terms  like  cerebral  oxygenation,  tissue  hemo¬ 
globin  oxygenation,  and  mean  hemoglobin 
oxygenation  are  used  in  the  literature7’24’26  to  indi¬ 
cate  this  concept.  Although  tissue  hemoglobin  oxy¬ 
genation  is  not  rigorous  because  hemoglobin 
molecules  are  located  in  only  blood,  the  term  is  used 
specifically  to  differentiate  between  the  hemoglobin 
saturation  in  the  tissue  vascular  bed,  as  measured  by 
the  NIR  method,  and  the  arterial  hemoglobin  satu¬ 
ration  Sa02,  as  measured  by  a  pulse  oximeter. 

The  goal  of  this  paper  is  to  demonstrate  the  NIR 
technique  as  a  real-time,  noninvasive  means  of  mon¬ 
itoring  hemoglobin  oxygenation  dynamics,  i.e., 
changes  in  the  concentrations  of  total  hemoglobin 
(Hb;)  and  oxygenated  hemoglobin  (Hb02),  in  the  vas¬ 
cular  bed  of  breast  and  prostate  rat  tumors  in  re¬ 
sponse  to  respiratory  challenge.  Compared  with 
previous  NIR  studies  of  tumors  in  vivo,  our  approach 
has  the  following  features:  (1)  The  transmission 
mode,  as  opposed  to  the  reflectance  mode  used  by 
Hull  et  al .,22  interrogates  deeper  regions  (central 
parts)  of  the  tumor.  (2)  Only  two  wavelengths,  as 
opposed  to  the  spectrum  of  300-1100  nm  used  by 
Steen  et  al.,21  are  employed  and  provide  a  fast  and 
low-cost  instrument.  (3)  A  source- detector  separa¬ 
tion  of  1—2  cm  interrogates  a  large  tumor  noninva- 
sively,  as  opposed  to  the  needlelike  probe  used  by 
Steinberg  et  al.20  More  innovatively,  on  the  basis  of 
the  experimental  observation  of  tumor  hemoglobin 
oxygenation  dynamics,  we  developed  a  tumor  he- 
moperfusion  model  that  provides  important  insight 
into  tumor  blood  perfusion. 

This  paper  is  organized  as  follows:  In  Section  2, 
we  describe  our  animal  model,  the  NIR  instrument, 
and  the  algorithm  for  calculations  of  tumor  blood 
oxygenation.  In  Section  3,  we  show  experimental 
results  measured  from  both  breast  and  prostate  tu¬ 
mors  under  respiratory  interventions  and  calculate 
time  constants  for  the  hemoglobin  oxygenation  dy¬ 
namics  of  the  tumors.  In  Section  4,  we  develop  a 
tumor  hemoperfusion  model  to  interpret  the  experi¬ 
mental  data  obtained  in  the  tumor-intervention  stud¬ 


ies  and  to  relate  the  time  constants  to  tumor  blood 
perfusion.  Finally,  in  Section  5,  we  discuss  the  re¬ 
sults,  the  future  extensions,  and  the  potential  uses  of 
the  NIR  technique  as  a  novel  diagnostic-prognostic 
tool  for  tumor  therapy  and  cancer  research. 

2.  Materials  and  Methods 

A.  Animal  Model  and  Measurement  Geometry 

NF13762  breast  tumor  was  implanted  in  adult  fe¬ 
male  Fisher  rats,  and  Dunning  prostate  adenocarci¬ 
noma  R3327-AT1  was  implanted  in  adult  male 
Copenhagen  rats.  The  tumors  were  grown  in 
pedicles27  on  the  forebacks  of  the  rats  until  the  tu¬ 
mors  were  approximately  1-2  cm  in  diameter.  Rats 
were  anesthetized  with  200- pi  ketamine  hydrochlo¬ 
ride  (100  mg/ml)  and  maintained  under  general  gas¬ 
eous  anesthesia  with  33%  inhaled  02  (0.3  dm3 /min 
02,  0.6  dm3/min  N20,  and  0.5%  methoxyflurane) 
through  a  mask  placed  over  the  mouth  and  nose. 
Tumors  were  shaved  to  improve  the  optical  contact 
for  transmitting  light.  Body  temperature  was  main¬ 
tained  with  a  warm-water  blanket.  In  some  cases,  a 
fiber-optic  pulse  oximeter  (Nonin,  Inc.,  Model  8600V) 
that  was  manufacturer  calibrated  was  placed  on  the 
hind  foot  to  monitor  arterial  oxygenation  Sa02,  and  a 
fiber-optic  probe  was  inserted  rectally  to  measure 
temperature.  The  tumor  volume  V  (in  centimeters 
cubed)  was  estimated  as  V  =  (4tt/3)  [(L  4-  W  +  H)/ 
6]3,  where  L,  W,  and  H  are  the  three  respective  or¬ 
thogonal  dimensions. 

Most  measurements  were  performed  with  33%  oxy¬ 
gen  as  inhaled  gas  to  achieve  a  stable  baseline  for  a 
period  of  5  to  15  min.  The  inhaled  gas  was  then 
switched  to  carbogen  (95%  oxygen,  5%  carbon  diox¬ 
ide)  for  at  least  20  min  and  then  switched  back  to  33% 
02  for  approximately  15  min.  The  complete  cycle 
lasted  1  hour.  Sometimes  repeated  carbogen  inter¬ 
ventions  were  performed  sequentially  to  evaluate  the 
reproducibility  of  the  time  profiles  of  the  tumors.  In 
certain  cases  alternative  gases  were  used,  as  defined 
in  the  results  and  figures,  and  some  rats  were  sacri¬ 
ficed  by  KCl-induced  cardiac  arrest. 

Figure  1  shows  the  measurement  geometry:  Hor¬ 
izontally,  the  delivering  and  the  detecting  fiber  bun¬ 
dles  were  face  to  face  in  the  transmittance  mode,  and 
both  were  in  contact  with  the  tumor  surface  without 
hard  compression.  The  separation  of  the  two  bundle 
surfaces  was  between  1.0  and  2.5  cm,  depending  on 
the  tumor  size.  Vertically,  the  two  bundle  tips  (with 
diameters  of  0.5  cm)  were  placed  around  the  middle  of 
the  tumor.  Thus  the  current  setup  of  the  probes 
provides  an  optimal  geometry  for  the  NIR  light  to 
interrogate  deep  tumor  tissue  with  minimal  interfer¬ 
ence  from  the  foreback  of  the  rat. 

B.  Near- Infrared  Instrument  and  Data  Analysis 

As  shown  in  Fig.  1,  we  used  a  homodyne  frequency- 
domain  photon-migration  system28*29  that  was  capa¬ 
ble  of  determining  the  amplitude  and  the  phase 
changes  of  amplitude-modulated  light  passing 
through  tumors.  In  this  setup  a  rf  source  modulates 
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Fig.  1.  Experimental  setup  of  a  one-channel,  NIR,  frequency- domain  IQ  instrument  for  tumor  oxygenation  measurement.  PMT, 
photomultiplier  tube  for  detecting  light;  IQ,  in-phase  and  quadrature  demodulator  for  retrieving  amplitude  and  phase  information;  Low 
Pass,  low-pass  filter.  The  5-mm-diameter  fiber  bundles  deliver  and  detect  the  laser  light  through  the  tumor  in  transmittance  geometry. 


the  light  from  two  laser  diodes  (wavelengths  of  758 
and  782  nm)  at  140  MHz.  The  laser  light  passes 
through  a  combined  fiber-optic  bundle,  is  transmitted 
through  the  tumor  tissue,  and  is  collected  by  a  second 
fiber  bundle.  The  light  is  then  detected  by  a  photo¬ 
multiplier  tube  and  demodulated  with  a  commer¬ 
cially  available  in-phase  and  quadrature  (IQ) 
demodulator  chip  into  its  I  and  Q  components.  After 
these  components  are  put  through  a  low-pass  filter 
they  can  be  used  to  calculate  the  amplitude  and  the 
phase  changes  caused  by  the  tumor.  These  steps  are 
expressed  mathematically  by 

I(t)  -  2 A  sin(otf  -1-  0)sin(otf) 

low 

=  A  eos(0)  -  A  cos (u>t  +  0) - >  7dc 

pass 

=  A  cos(0),  (1) 

Q{t)  =  2A  sin(w£  +  0)cos(utf) 

low 

=  A  sin(0)  +  A  sin(co£  +  0) - »  Qdc 

pass 

=  A  sin(0),  (2) 

0  =  tan_1(Qdc//dc),  (3) 

A  =  (Ij  +  Qj)v\  (4) 

where  A  and  0  are  the  amplitude  and  the  phase  of  the 
detected  light,  respectively,  and  a)  is  the  angular  mod¬ 
ulation  frequency  (=2^  X  140  MHz). 

The  two  laser  lights  were  time  shared,  and  the  con¬ 
trolling  process  and  the  data  acquisition  both  inter¬ 
faced  through  a  12-bit  analog-to -digital  board  (Real 
Time  Devices,  Inc.,  Model  AD2100)  with  a  maximum 
sampling  rate  of  4  Hz.28  However,  slower  sampling 
rates  were  used  in  measurements  to  compensate  for 
experimental  noise.  Simple  time  averaging  among  a 
few  adjacent  data  points  was  performed  during  data 
analysis  to  further  decrease  the  noise.  However,  data 
smoothing  was  not  applied  for  (1)  calculating  the  ex¬ 
perimental  uncertainty  (error  bars)  or  (2)  fitting  the 
time  constants  to  prevent  the  fast-changing  compo¬ 
nent  from  being  oversmoothed  and  overlooked.  The 
pulse-oximeter  data  were  not  averaged  because  they 
were  recorded  manually  and  appear  discrete  compared 


with  the  NIR  data.  The  experimental  uncertainties 
for  arterial  saturation  and  changes  in  hemoglobin  con¬ 
centrations  were  calculated  by  use  of  the  baseline  data 
taken  over  5-10  min  without  respiratory  perturbation 
to  the  rat.  Nonlinear  curve  fitting  based  on  the  Mar- 
quardt  algorithm30 >31  was  performed  by  use  of  Kalei- 
daGraph.32  The  software  also  provided  the  errors  (or 
uncertainties)  for  each  fitted  parameter,  the  optimized 
X2  values,  and  the  fitting  correlation  coefficient  R ,  to¬ 
gether  with  the  goodness  of  the  fit  R2.33  The  signifi¬ 
cance  of  changes  was  assessed  on  the  basis  of  Fisher 
protected-least-significant-difference  analysis  of  vari¬ 
ance  by  use  of  Statview  software. 

C.  Calculation  for  Changes  in  the  Hemoglobin 
Concentration 

It  is  well  known  that  the  NIRS  of  tissue  can  be  used 
to  determine  the  total  hemoglobin  concentration  Hty 
and  the  hemoglobin  oxygen  saturation  S02  of  an  or¬ 
gan  in  vivo.7  When  two  NIR  wavelengths  are  used 
(758  and  782  nm,  in  this  case)  it  is  assumed  that 
tissue  background  absorbance  is  negligible  and  that 
the  major  chromophores  in  organs  are  oxygenated 
and  deoxygenated  hemoglobin  molecules.  In  princi¬ 
ple,  because  the  IQ  system  can  give  both  phase  and 
amplitude  values,  we  should  be  able  to  obtain  abso¬ 
lute  calculations  of  Hb02,  Hb,  and  S02.7>28  How¬ 
ever,  given  the  tumor’s  small  size  and  large  spatial 
heterogeneity,  it  is  very  difficult  to  obtain  such  ab¬ 
solute  quantification  accurately  with  conventional 
algorithms34  that  are  based  on  the  diffusion  approx¬ 
imation.  Instead,  on  the  basis  of  the  modified  Beer- 
Lambert  law,  we  can  use  the  amplitude  of  the  light 
transmitted  through  the  tumor  to  calculate  concen¬ 
tration  changes  in  Hb02,  Hb,  and  Hb,  (expressed  as 
AHb02,  AHb,  AHty,  respectively)  of  the  tumor  that 
are  caused  by  respiratory  intervention.  These 
changes  can  be  derived25  and  expressed  as  (see  Ap¬ 
pendix  A  for  derivations  and  justifications) 

AHb  =  Hb(transient)  -  Hb(baseline) 

~6{k-iog(^) 

£(€Hb  €Hb02  -  €Hb  eHb02) 
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AHb02  =  Hb02(transient)  -  Hb02 (baseline) 

^(€Hb  €Hb02 ~  eHb  €Hb02) 

where  e^h  and  €nbo2  are  extinction  coefficients35  of 
deoxygenated  and  oxygenated  hemoglobin,  respec¬ 
tively,  at  wavelength  X;  the  variable  Ab  is  a  constant 
amplitude  of  baseline;  At  is  the  transient  amplitude 
under  measurement;  and  L  is  the  optical  path  length 
between  the  source  and  the  detector. 

Using  the  approach  suggested  by  Cope  and 
Delpy,10  we  can  express  L  as  L  =  DPF  X  d ,  where  d 
is  the  direct  source- detector  separation  in  centime¬ 
ters  and  DPF  is  the  ratio  between  the  optical  path 
length  and  the  physical  separation  and  is  tissue  de¬ 
pendent.  The  DPF  for  tumors  has  not  been  well 
studied;  for  simplicity,  we  assume  the  DPF  to  be  1  in 
our  calculations.  The  justification  for  this  simplifi¬ 
cation  is  given  in  Section  5.  After  substituting  the 
extinction  coefficients35  at  758  and  782  nm  in  Eqs.  (5) 
and  (6)  with  values  of  e^fb  =  0.359,  €^0;j  -  0.1496,  ejfb2 
=  0.265  and  e^bcx  =  0.178,  respectively,  in  units  of 
inverse  millimoles  times  inverse  centimeters,  we  ar¬ 
rive  at 

a  tti_  [7.34  log(At/A,)758  -  6. 17  log(At/A()782] 


Fig.  2.  Results  of  a  drift  test  of  the  NIR  instrument  by  use  of  a 
tissue  phantom.  The  thicker  solid  curve  represents  relative 
changes  in  the  oxygenated  hemoglobin  concentration,  i.e.,  AHb02, 
and  the  thinner  solid  curve  represents  relative  changes  in  the  total 
hemoglobin  concentration,  i.e.,  AHbr  AHb02  and  AHb,  were  cal¬ 
culated  by  use  of  Eqs.  (8)  and  (9),  respectively. 


centrations,  as  calculated  from  Eqs.  (8)  and  (9).  In 
this  example,  the  standard  deviations  over  the  entire 
period  of  100  min  were  less  than  0.007  and  0.004  mM, 
respectively,  for  AHb02  and  AHb,.  F urthermore,  we 
calculated  uncertainties  for  both  of  these  quantities 
on  the  basis  of  the  propagation  of  errors,  and  the 
results  are  consistent  with  those  shown  in  Fig.  2. 


AHb02  = 

[-10.92  log(A„/At)758  +  14.80  log(A,/Af)782] 

L 

AHb,  =  A(Hb02  +  Hb) 

[-3.58  logCAfc/A,)758  +  8.63  log(At/A()782] 


where  the  units  are  in  millimoles.  Equations  (7) 
and  (8)  permit  the  calculation  of  changes  in  Hb  and 
Hb02  that  are  due  to  respiratory  challenge,  respec¬ 
tively,  whereas  Eq.  (9)  quantifies  a  relative  increase 
in  the  total  hemoglobin  concentration  that  is  caused 
by  the  intervention.  The  last  quantity  also  reflects  a 
change  in  blood  volume  because  it  is  proportional  to 
the  total  Hb  concentration. 

3.  Results 

A.  Instrument  Drift  Tests 

The  stability  of  the  NIR  instrument  was  tested  in 
terms  of  baseline  drift  after  a  warm-up  period  of 
30  min  by  use  of  a  tissue  phantom25’36  with  stable 
optical  properties.  Figure  2  shows  an  example  of 
a  phantom  measurement  that  displays  the  variation 
of  relative  changes  in  apparent  Hb02  and  Hb,  con- 


B.  Breast  Tumors 

Figure  3(a)  shows  the  results  taken  from  a  breast 
tumor  (4.5  cm3)  with  a  source- detector  separation  of 

1.8  cm.  The  data  were  smoothed,  and  the  measure¬ 
ment  uncertainties  are  shown  at  only  discrete  loca¬ 
tions.  The  figure  shows  the  relative  changes  in  total 
hemoglobin  concentration  AHb,  and  oxygenated  he¬ 
moglobin  concentration  AHb02.  The  arterial  Hb 
saturation  was  also  obtained  to  show  a  relatively 
rapid  change  in  arterial  signals  when  the  inhaled  gas 
was  switched  from  33%  02  to  carbogen.  Respiratory 
challenge  caused  a  sharp  rise  in  AHb02  (p  <  0.01 
after  1  min,  p  <  0.0001  by  1.5  min)  that  was  followed 
by  a  further  slow,  gradual,  but  significant,  increase 
over  the  next  25  min  (p  <  0.001).  AHb,  also  changed 
significantly  (p  <  0.001)  within  the  first  minute,  but 
the  total  change  was  only  approximately  10%  of  that 
of  AHb02.  Given  the  exponential  appearance  of  the 
rising  part  of  AHb02,  we  used  single-exponential  and 
double-exponential  expressions  to  fit  the  data  in  the 
rising  portion  to  better  understand  and  quantify  the 
dynamic  features  of  AHb02.  The  unsmoothed  data 
and  the  fitted  curves  are  shown  in  Figure  3(b).  The 
double  exponential  appears  to  give  a  much  better  fit, 
as  is  confirmed  by  the  respective  R  values  (0.98  ver¬ 
sus  0.81).  Time  constants  of  0.18  ±  0.02  min  and 

27.8  ±  3.9  min  were  obtained  for  fast  and  slow  dy¬ 
namic  changes,  respectively,  in  the  tumor  Hb02  con¬ 
centration. 

Figure  4(a)  was  obtained  from  a  second  breast  tu¬ 
mor  (5.9  cm3)  with  a  source- detector  separation  of 
1.6  cm.  Here  AHb02  increased  rapidly  after  the  ini- 
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Fig.  3.  (a)  Results  obtained  with  the  NIR  instrument  from  a 

4.5-cm3  rat  breast  tumor  while  the  breathing  gas  was  switched 
from  33%  02  to  carbogen.  The  thicker  solid  curve  represents 
AHb02,  the  thinner  solid  curve  represents  AHb*,  and  the  dashed 
curve  with  the  filled  circles  represents  arterial  saturation,  (b) 
The  unsmoothed  data  and  the  fitted  curves:  The  solid  curves 
represent  the  best  fits  to  the  AHb02  data  at  the  rising  portion. 
The  best  fit  to  the  two  exponential  terms  is  0.143(1  —  exp[-(£  - 
12.5J/0.18]}  +  0.36(1  -  exp[-(Z  -  12.5)/27.8]},  with  R  =  0.98, 
whereas  the  best  fit  with  one  exponential  term  is  expressed  as 
0.322(1  -  exp [-(f  -  12.5)/5.1]},  with  R  =  0.81. 


tial  gas  switch  but  did  not  exhibit  the  continued  slow 
rise  afterward.  AHb,  was  found  to  increase  with  car¬ 
bogen  inhalation,  although  the  magnitude  was 
smaller  than  that  of  AHb02  during  the  period  of  the 
intervention.  Again,  changes  in  AHb02  were  mod¬ 
eled  by  a  single-exponential  term  that  yielded  a  time 
constant  of  2.00  ±  0.04  min  (. R  =  0.97)  and  by  a 
double-exponential  formula  with  two  time  constants 
of  0.8  ±  0.2  min  and  3.0  ±  0.3  min  (R  -  0.98).  In 
this  case  both  expressions  fit  the  data  well,  as  shown 
in  Fig.  4(b). 

To  demonstrate  the  reproducibility  of  the  dynamic 
changes  in  response  to  respiratory  challenge,  we  sub¬ 
jected  one  animal  to  repeat  carbogen  inhalation. 
Figure  5(a)  shows  measurements  taken  from  a  breast 
tumor  (6.7  cm3)  with  a  source- detector  separation  of 
2  cm.  In  this  case  air  with  1.2%  isoflurane  (anes¬ 
thetic)  was  used  as  the  baseline  instead  of  33%  02. 
This  figure  shows  a  very  consistent  pattern  in  two 
repeated  time  responses  with  a  fast  and  a  slow  in¬ 
crease  in  AHb02.  Here  AHb,  shows  a  similar  dy- 


Time  (min) 

Fig.  4.  (a)  Results  obtained  with  the  NIR  instrument  from  a 

5.9-cm3  rat  breast  tumor  while  the  breathing  gas  was  switched 
from  33%  02  to  carbogen.  (b)  The  solid  curves  show  the  best  fits 
to  the  Hb02  data  at  the  rising  portion.  In  this  case  a  double 
exponential  with  values  of  0.09(1  -  exp[-(£  -  9.8)/0.8]}  + 
0.16(1  -  exp[-(*  -  9.8)/3.0]},  with  R  =  0.98,  and  a  single  expo¬ 
nential  with  a  value  of  0.250(1  -  exp[-(£  -  9.8)/2.00]},  with  R  = 
0.97,  provided  similarly  good  fits. 


namic  pattern,  i.e.,  a  rapid  rise  followed  by  a  slow 
continuation.  Figures  5(b)  and  5(c)  show  the  un¬ 
smoothed  data  together  with  the  fitted  curves  for  the 
rising  portions  of  the  two  repeated  increases  in 
AHb02.  Again,  the  double-exponential  expression 
with  two  time  constants  produced  much  better  fits 
than  did  the  single-exponential  term  in  both  pro¬ 
cesses  with  two  averaged  time  constants  of  tx  (mean) 
=  0.26  ±  0.11  min  and  t2  (mean)  =  8.2  ±  1.8  min. 
Individual,  respective  time  constants  and  coefficients 
are  summarized  in  Table  1.  Furthermore,  single¬ 
exponential  and  double-exponential  expressions  were 
fitted  to  obtain  time  constants  for  the  decay  processes 
after  the  inhaled  gas  was  switched  repeatedly  back  to 
the  baseline  conditions.  Similarly,  the  double¬ 
exponential  expression  fits  the  data  better  with  two 
mean  time  constants  of  Tldecay  (mean)  =  0.17  ±  0.07 
min  and  T2decay  (mean)  =  12.2  ±  0.7  min  for  the  two 
decay  processes. 

To  further  validate  our  experimental  observations, 
we  subjected  some  rats  to  cardiac  arrest  (with  KC1)  to 
observe  the  changes  in  Hb02  and  Hb,  on  death.  Fig¬ 
ure  6  shows  an  example  of  cardiac  arrest  on  a  rat  with 
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(b) 


Time  (min) 
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Fig.  5.  (a)  Relative  changes  in  the  Hb02  detected  with  the  NIR 

instrument  from  a  rat  breast  tumor  (6.7  cm3)  while  the  breathing 
gas  was  alternated  between  air  (21%  02)  and  carbogen.  The  best 
fits  to  the  Hb02  data  by  use  of  both  the  double-exponential  and  the 
single-exponential  expressions  for  (b)  the  first  and  (c)  the  second 
respiratory  challenges  are  shown.  The  fitted  equations  that  were 
obtained  from  (b)  are0.232{l  -  exp[-(£  -  19.5)/0.18]}  +  0.368{1  - 
exp[-(£  -  19.5)/6.93]},  with  R  =  0.98,  and  0.550(1  -  exp[-(£  - 
19.5)/2.80]},  with  R  =  0.89,  respectively.  The  fitted  equations 
that  were  obtained  from  (c)  were  0.321(1  —  exp[— (t  —  67)/0.332]} 
+  0.233(1  -  exp[-(*  -  67)/9.47]},  with  R  =  0.99,  and  0.485(1  - 
exp[-(i  -  67)/1.38]},  with  R  =  0.81,  respectively. 
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Fig.  6.  Influence  of  KCl-induced  cardiac  arrest  on  the  values  of 
Hb02  and  Hbt  of  a  breast  tumor  (5.3  cm3),  while  the  rat  was 
breathing  air. 


a  breast  tumor  (5.3  cm3).  Both  AHb,  and  AHb02 
dropped  significantly,  immediately  after  KCL  was 
admitted  intravenously.  Within  1  min  AHfy 
reached  a  plateau,  whereas  AHb02  decreased  rapidly 
within  the  first  30  s  and  then  was  followed  by  a  slow 
prolongation. 

C.  Prostate  Tumors 

Figure  7(a)  was  obtained  from  a  large  prostate  tumor 
(8.2  cm3).  In  common  with  the  breast  tumors, 
AHb02  showed  a  rapid  initial  increase  that  was  fol¬ 
lowed  by  a  slower  continuation.  AHb,  increased  rap¬ 
idly  and  then  reached  a  plateau.  Figure  7  (b)  shows 
that  the  double-exponential  equation  fits  the  un¬ 
smoothed  data  better  ( R  =  0.96)  than  does  the  single¬ 
exponential  term  (R  =  0.82).  Here  the  fast  and  the 
slow  time  constants  are  0.265  ±  0.007  min  and  6.02  ± 
0.15  min,  respectively. 

Figure  8(a)  was  obtained  from  another  large  pros¬ 
tate  tumor  (10.8  cm3  with  a  source- detector  separa¬ 
tion  of  2.5  cm).  Here  AHb02  displayed  a  gradual 
increase  throughout  the  entire  period  of  carbogen  in¬ 
halation,  whereas  the  increase  in  AHb,  was  consid¬ 
erably  delayed.  Variations  in  arterial  hemoglobin 
saturation  Sa02  are  also  shown  and  were  very  rapid 
in  comparison  with  AHb02,  in  common  with  Fig.  3. 
AHb02  dropped  rapidly  when  the  inhaled  gas  was 
switched  back  from  carbogen  to  33%  02.  Both  the 
single-exponential  and  the  double-exponential  ex¬ 
pressions  were  used  to  obtain  time  constants  for  the 
rising  portion  of  AHb02  that  was  due  to  carbogen 
intervention.  In  this  case  both  expressions  gave 
equally  good  fits,  as  shown  in  Fig.  8(b)  and  Table  1. 
For  the  decay  process,  we  obtained  Tldecay  =  0.6  ±  0.2 
min  and  T2decay  =  6.6  ±  1.7  min  with R  =  0.94  for  the 
double-exponential  fitting,  whereas  the  single¬ 
exponential  fitting  resulted  in  t  =  2.8  ±  0.4  min  with 
R  =  0.88.  For  comparison  the  rat  was  also  chal¬ 
lenged  with  100%  02. 

In  summary,  we  observed  dynamic  changes  in 
Hb02  that  were  due  to  carbogen  intervention  for  both 
breast  and  prostate  tumors.  In  most  cases  these 
changes  were  modeled  better  by  a  double-exponential 


(b) 

Fig.  7.  (a)  Influence  of  respiratory  challenges  (switching  from  air 

to  carbogen)  on  the  values  of  Hb02  and  Hb,  of  a  large  rat  prostate 
tumor  (8.2  cm3),  (b)  The  best-fitted  equations  are  0.090{1  - 
exp[-(f  -  12)/0.265]}  +  0.064{1  -  exp| [-(f  -  12)/6.02]},  with  R  = 
0.96,  and  0.140{1  -  exp[-(£  -  12)/1.13]},  with  R  =  0.82,  for  the 
double-exponential  and  the  single-exponential  expressions,  respec¬ 
tively. 


expression  with  a  fast  and  a  slow  time  constant  than 
they  were  by  a  single-exponential  fitting.  Dynamic 
changes  in  arterial  saturation  preceded  those  in 
Hb02.  The  detailed  parameters  regarding  tumor 
size,  fitted  time  constants,  corresponding  magni¬ 
tudes,  and  R2  are  listed  in  Table  1. 

4.  Model  for  the  Blood  Oxygenation  Dynamics  of 
Tumors 

As  was  shown  in  Section  3,  the  temporal  changes  in 
Hb02  caused  by  respiratory  challenge  can  be  fitted 
with  an  exponential  equation  that  has  either  one  or 
two  time  constants  (fast  and  slow).  In  this  section, 
we  further  derive  and  simplify  a  hemoperfusion 
model  to  interpret  these  time  constants  and  to  corre¬ 
late  the  experimental  findings  with  the  physiology  of 
the  tumors. 

To  develop  the  model,  we  follow  an  approach  used 
to  measure  regional  cerebral  blood  flow  (rCBF)  with 
diffusible  radiotracers,  as  originally  developed  by 
Kety37  in  the  1950’s.  The  basic  model  was  modified 
in  a  variety  of  ways  to  adapt  it  to  positron  emission 
tomography  studies.38  39  By  analogy,  we  can  evalu¬ 
ate  tumor  hemodynamics  such  as  tumor  blood  flow 
(perfusion)  by  using  the  respiratory-intervention  gas 
as  a  tracer. 
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(a) 


Fig.  8.  (a)  Variations  in  Sa02,  AHb02,  and  Hb,  of  which  the  latter 
two  were  detected  with  the  NIR  instrument,  from  a  large  rat 
prostate  tumor  (10.8  cm3)  during  respiratory  challenge,  (b)  The 
solid  curves  represent  the  best  fits  to  the  AHb02  data  at  the  rising 
portion  during  carbogen  inhalation.  The  fitted  equations  are 
0.004{1  -  exp[-(f  -  21.5)/0.3]}  +  0.38{1  -  exp [-(*  -  21.5)/15.6]} 
and0.37{l  -  exp[-(£  -  21.5)/14.8]}  for  the  double-exponential  and 
the  single-exponential  expressions,  respectively,  with  R  =  0.97  for 
both. 


In  general,  the  Fick  principle  can  be  stated  as  fol¬ 
lows38  :  The  rate  of  change  of  tracer  concentration  in 
a  regional  area  of  an  organ  equals  the  rate  at  which 
the  tracer  is  transported  to  the  organ  in  the  arterial 
circulation  minus  the  rate  at  which  it  is  carried  away 
into  the  venous  drainage,  i.e., 


dCt 

at 


=  f(Ca-C0), 


(10a) 


In  Eq.  (10b),  f  and  X  are  constants,  whereas  Ct  is  a 
time-dependent  variable  that  is  written  as  Ct(t).  In 
principle,  the  arterial-tracer  concentration  Ca  is  a 
time-varying  quantity.  If  a  certain  concentration  of 
the  arterial  tracer  is  administrated  continuously 
starting  at  time  0,  Ca  can  be  expressed  mathemati¬ 
cally  as  a  constant  value  of  Ca(0)  after  time  0.  Then 
Eq.  (10b)  can  be  solved  as 

Ct(t)  =  XCfl(0)[l  -  exp (-/*/X)].  (ID 

Equation  (11)  indicates  that,  at  time  t  after  the  onset 
of  tracer  administration,  the  local  tissue  (traditional¬ 
ly  brain)  Ct(t)  concentration  depends  on  the  blood 
flow  f,  the  arterial  time-activity  curve  Cc(0),  and  the 
partition  coefficient  X. 

In  response  to  respiratory  intervention,  a  sudden 
small  change  is  introduced  into  the  arterial  02  satu¬ 
ration  Sa02,  and  the  resulting  increase  in  arterial 
Hb02  concentration  (AHb02artery)  can  be  considered 
as  an  intravascular  tracer.40  Following  Kety’s 
method  and  assuming  that  changes  in  dissolved  02 
are  negligible,40  we  have 


—  (AHb02vasculature)  = 

d£ 


(12) 

where  f  still  represents  blood  flow  (or  perfusion  rate) 
and  y  is  defined  as  a  vasculature  coefficient  of  the 
tumor.  The  coefficient  y  is  the  ratio  of  the  Hb02 
concentration  change  in  the  vascular  bed  to  that  in 
veins  and  equals  (AHb02vascu,ature)/(AHb02vein). 
This  definition  implies  that  a  change  in  the  venous 
blood  oxygenation  AHb02vein  is  proportional  to  a 
change  in  the  Hb  oxygenation  in  the  vascular  bed, 
AHb02vasculature. 

In  Eq.  (12),  f  and  y  are  constants,  whereas 
AHb02vascu  ature  is  a  time-dependent  variable.  By 
analogy  to  Eq.  (11),  AHb02vasculature  can  be  solved 
rigorously  given  a  constant  input  H0  for  AHb02artery 
after  time  0.  Our  data  (Figs.  3  and  8)  demonstrate 
that  changes  in  the  arterial  Hb02  (Sa02)  are  much 
faster  than  in  the  vascular  bed.  Then  solving  Eq. 
(12)  leads  to 

AHb02vasc,,lotureW  =  yH0[l  -  exp  (-/i/7)].  (13) 


where  fis  the  blood  flow  (or  perfusion  rate),  Ct  is  the 
tracer  concentration  in  tissue,  and  Ca  and  Cv  are  the 
time-varying  tracer  concentrations  in  the  arterial  in¬ 
put  and  the  venous  drainage,  respectively.  Ca  can 
be  measured  from  a  peripheral  artery,  but  Cv  is  rel¬ 
atively  difficult  to  obtain  regionally.  Therefore  a 
brain- blood  partition  coefficient,  X  =  Ct/Cv ,  was  de¬ 
veloped  by  Kety,  and  Eq.  (10a)  becomes37 


Equation  (13)  indicates  that,  at  time  t  after  the  onset 
of  respiratory  intervention,  the  change  in  oxygenated 
hemoglobin  concentration  in  the  tumor  vasculature 
AHb02vasculature  (£)  depends  on  the  blood  perfusion 
rate  f *  the  arterial  oxygenation  input  HQy  and  the 
vasculature  coefficient  of  the  tumor  y. 

As  indicated  by  Eq.  (8),  our  NIR  instrument  is  able 
to  measure  an  increase  in  the  vascular  Hb02  concen¬ 
tration  AHb02vasculature.  Equation  (13)  gives  an  ex¬ 
ponential  of  the  same  form  as  that  used  to  fit  our 
experimental  data,  indicating  that  the  measured 
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Fig.  9.  Schematic  diagram  showing  (1)  a  tumor  model  with  two 
vascular  perfusion  regions,  (2)  the  source  and  the  detector  fibers 
and  their  geometry  with  respect  to  the  tumor  model,  and  (3)  the 
light  patterns  that  propagate  in  the  tumor  tissue.  A  represents  a 
portion  of  the  detected  signal,  which  interrogates  the  well-perfused 
region,  and  B  represents  another  portion  of  the  detected  signal, 
which  passes  mainly  through  the  poorly  pei-fused  region.  We 
have  assumed  that  the  total  detected  signal  is  the  sum  of  A  and  B. 

time  constant  is  associated  with  the  blood  perfusion 
rate  f  and  the  vasculature  coefficient  7  of  the  tumor  in 
the  measured  area.  If  the  measured  volume  in¬ 
volves  two  distinct  regions,  we  then  involve  with  two 
different  blood-perfusion  rates  fx  and  f2>  two  different 
vasculature  coefficients  y1  and  72,  or  all  four.  Here 
it  is  reasonable  to  assume  that  the  measured  signal 
results  from  both  regions,  as  illustrated  in  Fig.  9. 
Consequently,  Eq.  (13)  can  be  modified  with  a  double¬ 
exponential  expression  and  two  time  constants  as 

AHb02vasculature(f)  =  7,ffo[  1  -  exp( -A  */«,,)] 

+  y2H0[l  -  exp (-f2t/y2)] 

=  A1[l-exp(-A«/'Vi)] 

+  A2[l  -  exp(- ftt /y2)],  (14) 

where  fx  and  yx  are  the  blood-perfusion  rate  and  the 
vasculature  coefficient,  respectively,  in  region  1,  f2 
and  72  are  the  same  for  region  2,  Ax  =  yx H0,  and  A2  = 
y2 H0.  The  two  time  constants  are  equal  to  tx  =  yjfx 
and  tx  =  y2/f2.  Then,  if  Al9  A2,  and  the  two  time 
constants  are  determined  from  our  measurements, 
we  arrive  at  the  ratios  for  the  two  vasculature  coef¬ 
ficients  and  the  two  blood-perfusion  rates: 

h  _^-i  /i  _  Ax/ A2  ,ieA 

A  *  /•  /  ’  (lb) 

72  A2  /2  Ti/Tg 

With  these  two  ratios,  we  can  obtain  insight  into 
the  tumor  vasculature  and  blood  perfusion.  For  ex¬ 
ample,  a  ratio  of  7i/72  near  1  from  a  measurement 
implies  that  the  vascular  structure  of  the  measured 
tumor  volume  is  rather  uniform.  Then  the  coexist¬ 
ence  of  two  time  constants  reveals  two  mechanisms  of 
regional  blood  perfusion  in  the  tumor.  A  large  time 


constant  implies  slow  perfusion  through  a  poorly  per¬ 
fused  area,  whereas  a  small  time  constant  indicates 
fast  perfusion  through  a  well-perfused  area.  In  the 
meantime,  the  ratio  of  the  perfusion  rates  in  these 
two  areas  can  also  be  obtained  quantitatively.  Fur¬ 
thermore,  a  ratio  of  7i/72  >  1  (i.e.,  AJA2  >  1)  means 
that  the  measured  signal  results  more  from  region  1 
than  from  region  2  within  the  measured  tumor  vol¬ 
ume.  Therefore,  by  studying  tumor  blood  oxygen¬ 
ation  dynamics  and  obtaining  time  constants 
together  with  their  amplitudes,  we  can  gain  impor¬ 
tant  information  on  regional  blood  perfusion  and  vas¬ 
cular  structures  of  the  tumor  within  the  measured 
volume. 

Our  experimental  data  (Table  1)  reveal  that  all  the 
measurements  can  be  fitted  with  the  double¬ 
exponential  model  equivalently  to  or  better  than  the 
single-exponential  fitting.  Ratios  of  7i/72> 

and  fi/f2  are  also  shown  in  Table  1  for  respective 
cases. 

5.  Discussion 

Using  NIRS,  we  have  measured  relative  changes  in 
Hb,  and  Hb02  in  breast  and  prostate  rat  tumors  in 
response  to  respiratory  intervention.  We  have  ob¬ 
served  that  respiratory  challenge  caused  the  Hb02 
concentration  to  rise  promptly  and  significantly  in 
both  breast  and  prostate  tumors  but  that  the  total 
concentration  of  hemoglobin  sometimes  increased 
and  sometimes  remained  unchanged.  The  dynamic 
changes  of  tumor  oxygenation  can  be  modeled  by  ei¬ 
ther  one  exponential  term  with  a  slow  time  constant 
or  two  exponential  terms  with  fast  and  slow  time 
constants.  This  relation  suggests  that  there  may  be 
two  vascular  mechanisms  in  the  tumor  that  are  de¬ 
tected  by  the  NIRS  measurement.  As  indicated  by 
Eqs.  (13)  and  (14),  these  time  constants  are  inversely 
proportional  to  the  blood-perfusion  rates  of  the  mea¬ 
sured  volumes  of  the  tumors.  Based  on  the  double¬ 
exponential  model,  determination  of  the  two  time 
constants  and  their  corresponding  amplitudes  allows 
us  to  determine  the  relations  between  the  two  perfu¬ 
sion  rates  and  between  the  vascular  structures,  as 
expressed  in  Eq.  (15).  Further  investigation  with 
more  measured  quantities  may  lead  to  quantification 
of  each  parameter  individually  by  use  of  the  NIR 
technique. 

To  develop  a  model  for  interpreting  the  NIR  data 
taken  during  carbogen  inhalation,  we  have  defined  a 
vasculature  coefficient  y.  It  is  a  proportionality  fac- 
tor  between  AHb02vein  and  AHb02vasculature,  i.e., 
AHb02vein  =  AHb02vasculature/'y.  We  expect  that  y 
depends  on  (1)  the  oxygen  consumption  and  (2)  the 
capillary  density  of  the  tumor.  If  the  oxygen  con¬ 
sumption,  the  capillary  density,  or  both  of  the  tumor 
are  large,  changes  in  the  venous  Hb02  concentration 
will  be  small;  if  the  oxygen  consumption,  the  capillary 
density,  or  both  of  the  tumor  are  small,  changes  in  the 
venous  Hb02  concentration  will  be  large.  Further 
studies  are  necessary  to  learn  more  about  this  coef¬ 
ficient  and  to  confirm  our  speculation. 

Our  current  NIR  system  allows  us  to  quantify  the 
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ratio  of  y/f  by  using  the  single-exponential  model  or 
to  quantify  the  ratios  of  yl/y2  and  f1/f2  by  using  the 
double-exponential  model.  We  can  obtain  impor¬ 
tant  information  on  the  blood  perfusion  of  the  tumor: 
a  large  time  constant  usually  represents  slow  blood 
perfusion,  whereas  a  small  time  constant  indicates 
fast  blood  perfusion.  The  coexistence  of  two  time 
constants  implies  a  combination  of  well-perfused  and 
poorly  perfused  mechanisms  of  blood  perfusion.  In¬ 
deed,  some  tumor  lines  have  only  20%  to  85%  of  ves¬ 
sels  perfused.41  Furthermore,  tumor  structures  and 
oxygen  distribution6’42  are  highly  heterogeneous. 
Therefore  it  is  likely  that  our  measurement  detects  a 
well-perfused  region,  or  a  poorly  perfused  region,  or  a 
mixture  of  both  in  the  tumor,  depending  on  the  posi¬ 
tion  or  the  location  of  the  source  and  the  detector  of 
the  NIR  instrument. 

Hull  et  al22  reported  carbogen-induced  changes  in 
rat  mammary  tumor  oxygenation  by  using  spatially 
resolved  NIRS.  To  compare  our  results  to  theirs,  we 
applied  our  curve-fitting  procedure  to  their  published 
hemoglobin  saturation  curve  and  obtained  a  time 
constant  of  0.27  min  (or  16  s)  for  the  rising  edge, 
which  is  consistent  with  our  fast  component.  Their 
data  do  not  show  a  slow  component,  suggesting  that 
their  measurement  was  dominated  by  active  tumor 
vasculature.  This  difference  may  be  explained  as 
follows: 

(1)  The  tumor  volume  mentioned  in  the  paper  by 
Hull  et  al 22  was  approximately  1.5  cm3,  much 
smaller  than  the  volumes  of  the  tumors  that  we  mea¬ 
sured  (in  Figs.  3  to  8,  the  largest  tumor  was  10.8  cm3, 
whereas  the  smallest  tumor  was  4.5  cm3). 

(2)  Their  measurement  was  in  reflectance  geome¬ 
try  with  multiple  detectors  located  at  distances  1  to 
20  mm  away  from  the  source,  and  in  the  calculation 
the  tumor  was  assumed  to  be  homogenous  in  order  to 
use  diffusion  theory.  Thus  their  measurement  was 
more  sensitive  to  the  superficial  area  of  the  tumor, 
emphasizing  the  tumor  periphery,  which  is  often  bet¬ 
ter  vascularized  than  the  central  part  of  the  tu¬ 
mor.42*43 

The  fast  component  observed  in  our  tumors  is  con¬ 
sistent  with  the  rapid  changes  detected  in  Sa02  in  the 
leg  as  measured  by  use  of  the  pulse  oximeter,  provid¬ 
ing  further  evidence  that  this  relates  to  the  well- 
vascularized,  highly  perfused  region  of  the  tumor. 

The  data  shown  in  Fig.  8  in  this  paper  are  repre¬ 
sentative  of  the  measurement  of  a  poorly  perfused 
region  in  which  the  measured  tumor  was  large  and 
the  portion  for  the  fast  oxygenation  response  was 
small.  The  data  given  in  Figs.  3-5  and  7  resulted 
from  a  mixture  of  well-perfused  and  poorly  perfused 
areas  in  the  tumors  and  exhibited  a  mixture  of  fast 
and  slow  oxygenation  responses  to  hyperoxygen  con¬ 
ditions.  It  is  reasonable  to  expect  that  larger  tumors 
have  more  poorly  perfused  regions  than  do  smaller 
tumors.  The  time  constant  of  the  slow  component 
observed  here  approaches  that  observed  previously 
for  changes  in  tissue  p02  in  an  ATI  prostate  tumor 


measured  by  use  of  19F  nuclear  magnetic  resonance 
spectroscopy  to  interrogate  interstitial  oxygenation.44 
In  general,  well -oxygenated  tumor  regions  had  a 
large  and  rapid  response  to  respiratory  challenge, 
whereas  poorly  oxygenated  regions  were  much  more 
sluggish.45  One  would  indeed  expect  changes  in 
vascular  oxygenation  to  precede  changes  in  the  tis¬ 
sue,  and  combined  investigations  by  NIRS  and  nu¬ 
clear  magnetic  resonance  spectroscopy  in  the  future 
will  provide  further  insight  into  the  delivery  of  oxy¬ 
gen  to  tumors. 

Dynamic  changes  in  vascular  oxygenation  have 
been  assessed  previously  by  several  other  techniques. 
Following  the  infusion  of  Green  2W  dye  intrave¬ 
nously  into  EMT-6  tumor-bearing  mice,  Vinogradov 
et  al.4  were  able  to  image  changes  in  the  surface 
vascular  p02.  On  switching  from  air  to  carbogen 
inhalation,  they  observed  a  very  rapid  increase  in  p02 
with  a  rate  similar  to  the  fast  component,  which  we 
have  seen  here.  Although  the  phosphorescence 
method  provides  vascular  p02,  NIR  methods  gener¬ 
ally  provide  Hb02  or  S02  because  there  is  some  un¬ 
certainty  in  the  local  affinity  of  hemoglobin  for  tumor 
oxygen:  the  p02-S02  dissociation  curve  is  subject  to 
pH,  temperature,  and  other  allosteric  effectors,  such 
as  2,3-diphosphoglycerate  in  the  immediate  milieu. 
A  promising  new  approach  is  the  blood-oxygen-level- 
dependent  (BOLD)  contrast  *H  MRI,  which  is  sensi¬ 
tive  to  vascular  perturbations.  Robinson  et  al46 
explored  the  response  to  respiratory  challenge  in  var¬ 
ious  tumors  and  showed  reversible  regional  changes 
on  switching  from  air  to  carbogen  inhalation.  In 
common  with  our  NIR  data,  their  changes  were  often 
biphasic  with  a  large  change  occurring  within  the 
first  2  min  and  followed  by  slower  increases.46  How¬ 
ever,  interpreting  the  BOLD  MRI  results  is  compli¬ 
cated  by  variations  in  vascular  volume  and  flow,  and 
there  is  no  direct  measure  of  Hb02  in  tumors. 

The  time  constants  are  not  source- detector  sepa¬ 
ration  sensitive.  Equations  (8)  and  (9)  have  demon¬ 
strated  that  AHb02  and  AHb,  are  proportional  to  1/d, 
where  d  is  the  source- detector  separation.  This  re¬ 
lation  indicates  that  a  different  d  value  will  only 
stretch  or  compress  the  entire  temporal  profile  of 
AHb02,  but  it  does  not  change  the  transient  behavior 
of  the  time  response.  The  same  argument  can  apply 
to  the  DPF.  In  this  study,  we  have  assumed  a  DPF 
ratio  of  1  for  simplicity.  If  the  DPF  value  is  larger 
than  1,  the  values  of  AHb02  and  AHb,  will  decrease 
by  a  factor  of  DPF.  But  this  modification  does  not 
affect  the  time  constants  tx  and  t2,  which  constitute 
the  dynamic  responses  of  AHb02  of  the  tumors  to 
respiratory  intervention. 

Given  the  evidence  for  intratumoral  heterogeneity 
from  MRI6*46  and  histology,47  we  believe  it  will  be 
important  to  advance  our  NIR  system  to  have  multi¬ 
ple  sources,  multiple  detectors,  or  both  to  study  not 
only  dynamic  but  also  spatial  aspects  of  blood  oxy¬ 
genation  in  tumor  vasculature.  Nonetheless,  we  be¬ 
lieve  the  preliminary  results  described  here  are  a 
proof  of  principle  for  the  technique,  laying  a  founda¬ 
tion  for  more  extensive  tests  to  correlate  tumor  size 
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where  d  is  the  source- detector  separation  and  fxa  and 
are  the  absorption  and  the  reduced  scattering  co¬ 
efficients,  respectively.  Because  in  the  NIR  region 
the  ix's  of  tissue  is  not  sensitive  to  either  wavelength 
or  perturbation,  we  assume  that  a  change  in  L  results 
from  only  a  change  in  |xa,  which  is  both  wavelength 
and  perturbation  dependent.  With  this  assumption, 
Eq.  (A7)  leads  to 


A L  =  lA,ia 
L  2\ya  ' 


(A8) 


Equation  (A8)  allows  us  to  determine  the  relative 
errors  of  A L/L  that  are  caused  by  (a)  the  wavelength 
dependence  of  |xa  and  (b)  the  perturbation  depen¬ 
dence  of  |xa  in  the  tumor.  For  case  (a),  we  calculated 
this  error  by  using 


758 

I h* 


782 


2^ 


758 

a 


under  the  baseline  and  the  perturbed  conditions;  for 
case  (b),  we  employed 


I ±a  (transient)  -  |x£  (baseline) 

2 (baseline) 

at  both  \  —  758  nm  and  \  =  782  nm  for  the  error 
calculation.  The  |xa  values  used  here  were  taken 
from  Hull  et  al.22  Although  the  rat  tumor  used  in 
their  study  was  different  from  ours,  the  absorption 
coefficients  of  the  tumors  should  be  in  a  similar  order 
and  follow  a  similar  dynamic  trend.  The  calculation 
shows  that,  with  758  and  782  nm  under  carbogen 
perturbation,  the  maximum  value  of  AL/L  is  12%. 
This  result  implies  that  the  assumption  of  a  constant 
path  length  that  was  used  for  Eqs.  (5)  and  (6)  gives 
rise  to  a  maximal  relative  error  of  12%  in  L. 

On  the  basis  of  Eqs.  (5)  and  (6)  [or  Eqs.  (7)  and  (8)], 
we  arrive  at  AX/X  =  —AL/L,  where X  can  be  AHb02, 
AHb,  or  AHbt.  Thus  the  assumption  of  a  constant 
path  length  leads  to  a  maximal  relative  error  of  12% 
for  the  magnitude  of  the  changes  th  at  we  detected 
with  regard  to  respiratory  challenge.  Although  12% 
is  not  completely  negligible,  the  measurement  and 
the  calculation  with  the  assumption  of  a  constant 
path  length  are  still  worthwhile.  Such  an  approach 
makes  it  possible,  as  a  first-order  approximation,  to 
quantify  the  AHb,  and  the  AHb02  of  tumors  under 
respiratory  intervention,  providing  deep  insight 
into  tumor  vascular  phenomena  and  mechanisms 
of  modulating  tumor  physiology  for  therapeutic 
enhancement. 
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ABSTRCT 

Oxygen  plays  a  key  role  in  tumor  therapy  and  may  be  related  to  tumor  development:  e.g., 
angiogenesis  and  metastasis.  Using  noninvasive  techniques  to  accurately  measure  tumor  oxygenation 
could  assist  in  developing  novel  therapies.  Here*  we  have  used  the  FREDOM  (Fluorocarbon 
Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  approach  based  on 
hexafluorobenzene  (HFB)  to  monitor  tissue  oxygen  tension  (p02 )  of  rat  breast  and  prostate  tumors 
and  compared  the  results  with  changes  in  tumor  vascular  hemoglobin  saturation  (s02)  and 
concentration  observed  using  a  new  dual  wavelength  homodyne  near-infrared  (NlR)  system.  The 
dynamic  changes  in  p02  and  SO2  were  assessed  while  rats  were  breathing  various  gases.  NIR  showed 
significant  changes  in  vascular  oxygenation  accompanying  respiratory  interventions.  WF  MR-EPI 
also  showed  significant  changes  in  tissue  p02  and  revealed  considerable  regional  heterogeneity  in 
both  absolute  values  and  rate  of  change  accompanying  interventions.  Generally,  changes  in  vascular 
s02  preceded  tissue p02,  particularly  for  smaller  tumors. 

Keywords:  Oxygen  tension.  Echo  planar  imaging,  MRI,  Tumor,  NIR  spectroscopy 


INTRODUCTION 

The  growth  and  development  of  tumors  are 
greatly  influenced  by  oxygen  tension  (p02),  e.g.,  tumor 
hypoxia  leads  to  increased  expression  of  vascular 
endothelial  growth  factor  (VEGF),  and  thus, 
angiogenesis  [1].  Hypoxia  reduces  radiosensitivity  [2], 
but  chemotherapeutic  approaches  have  been  proposed  to 
exploit  the  tumor  hypoxia  based  on  selective  cytotoxicity 
of  bioreductive  drugs  [3,  4],  In  addition,  increasing 
evidence  from  clinical  trials  has  revealed  that  poorly 
oxygenated  tumors  have  poor  prognosis  for  patients  [5, 
6],  Therefore,  accurate  measurements  of  oxygenation 
could  enhance  cancer  treatment  planning.  Here,  we 
present  two  methods  of  measuring  tumor  oxygenation: 
the  FREDOM  approach  to  measure  tumor  tissue  oxygen 
tension  (p02 )  and  NIR  spectroscopy  to  measure  changes 
in  tumor  vascular  hemoglobin  saturation  (s02)  and 
concentration  [Hb].  By  comparing  these  two  techniques, 
we  can  examine  die  relationship  between  tumor  tissue 
pO 2  and  vascular  s02- 
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The  FREDOM  approach  is  based  on  19F  PBSR- 
EPI  of  hexafluorobenzene  (HFB).  It  has  been  shown  that 
the  spin-lattice  relaxation  rate,  tfl  (1AF1),  is  linearly 
proportional  to  dissolved  oxygen  concentration  [7].  HFB 
offers  exceptional  sensitivity  to  changes  \np02  while 
having  little  response  to  temperature  [8].  Because  of 
structural  symmetry,  HFB  has  a  single  resonance  and 
thus,  is  free  from  chemical  shift  artifact,  providing  an 
optimal  signal-to-noise  ratio  (SNR).  Maps  of  the  tumor 
tissue  p02  were  obtained  in  8  minutes  with  millimeter 
resolution,  allowing  the  fate  of  individual  voxels  to  be 
traced.  NIR  spectroscopy  can  be  used  to  measure  tumor 
vascular  s02  because  the  absorption  coefficients  of 
deoxy-hemoglobin  differ  from  those  of  oxy-hemoglobin 
at  the  wavelengths  selected  (758  nm  and  782  nm)  [9]. 
The  system  has  many  attractive  features:  completely 
non-invasive,  inexpensive,  portable,  and  real-time. 


METHODS 

Tumor  Model 

NF  13762  breast  and  Dunning  prostate  R3327- 
AT1  adenocarcinomas  were  implanted  in  skin  pedicles 
on  the  forebacks  of  adult  female  Fischer  and  male 
Copenhagen  rats  (-250  g),  respectively,  as  described 
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previously  [10].  Once  the  tumors  reached  ~1  cm 
diameter,  the  rats  were  anesthetized  with  200  pi 
ketamine  hydrochloride  (100  mg/ml)  and  maintained 
under  general  gaseous  anesthesia  (33%  02, 66%  N20  and 
0.5%)  through  a  mask  placed  over  the  mouth  and  nose. 
Body  temperature  was  maintained  at  37°C  by  a  thermal 
blanket.  A  fiber  optic  pulse  oximeter  was  placed  on  the 
hind  foot  to  monitor  arterial  hemoglobin  saturation  and 
heart  rate,  and  a  fiber  optic  probe  was  inserted  rectally  to 
monitor  core  temperature.  NIR  spectroscopy  and  19F 
PBSR-EPI  measurements  were  then  performed 
sequentially,  while  inhaled  gas  was  alternated  between 
33%  02,  carbogen  (95%  02  +  5%  C02),  and  100%  02. 


NIR  Spectroscopy 

The  tumor  vascular  s02  was  assessed  by  NIR 
spectroscopy  using  a  new  dual  wavelength,  homodyne 
system  (wavelengths  758  nm  and  782  nm).  These 
wavelengths  were  selected  because  they  not  only  allow 
the  calculation  of  s02,  but  also  fall  into  the  range  of 
wavelengths  compatible  with  the  low  cost  photo 
multiplier  tube  (PMT).  The  system  uses  only  one  RF 
source  to  determine  amplitude  and  phase  changes  of 
light.  Figure  1  shows  a  schematic  diagram  of  the  system. 
An  RF  source  modulates  the  light  from  two  laser  diodes 
at  140  MHz.  The  light  passes  through  fiber  optic  cables, 
is  transmitted  through  the  tumor  tissue,  and  is  collected 
by  a  second  fiber  bundle.  The  light  is  then  detected, 
amplified,  filtered,  and  demodulated  into  I  and  Q 
components.  Amplitude  and  phase  changes  caused  by  the 
tumor  are  related  to  changes  in  hemoglobin 
concentration  [Hb]  and  hemoglobin  saturation  [Hb02], 
i.e.,  s02. 


Figure  1.  A  schematic  diagram  of  the  NIR  IQ  system 


To  obtain  the  absorption  coefficients  of  deoxy¬ 
hemoglobin  and  oxy-hemoglobin,  we  have  assumed 
background  absorbance  to  be  negligible  and  estimated 


the  absorption  coefficients  by  multiplying  the 
extinction  coefficients  for  deoxy-hemoglobin  and  oxy¬ 
hemoglobin  with  their  respective  concentrations. 

Pa75®  =  EH0758[Hb]  +  EHb02758[HbO2]  (1) 

Pa782  =  EHb782[Hb]  +  SHbC2782[Hb02]  (2) 

where  pa758  and  pa782  are  the  absorption  coefficients, 
sHb758  and  enb782  the  extinction  coefficients  for  deoxy- 
hemoglobin,  Enboz758  and  SHb02782  the  extinction 
coefficients  for  oxy-hemoglobin  at  the  wavelengths  758 
nm  and  782  nm,  respectively,  and  [Hb]  and  [HbOJ  are 
the  deoxy-hemoglobin  and  oxy-hemoglobin  concen¬ 
trations,  respectively. 

Although  the  IQ  system  does  give  both  phase 
and  amplitude  values,  given  the  tumors'  small  sizes  and 
our  fiber  configuration,  conventional  diffusion  theory 
doesn't  hold.  To  overcome  this  difficulty,  we  modified 
Beer-Lambert’s  law  and  used  the  amplitude  values  to 
find  trends  in  the  changing  absorption  coefficients. 

PaC  -  Pal  =  (l/L)*log  (Aj/Ac)  (3) 

where  Ai  is  the  initial  amplitude  (amplitude  of  baseline), 
Ac  the  current  amplitude,  and  L  the  optical  pathlength 
between  source  and  detector. 

By  manipulating  equations  1-3,  we  can  compute 
changes  in  blood  volume  and  saturation  between  the 
initial  state  and  the  intervention  state  from  the 
transmitted  amplitude  of  the  light  through  the  tumor. 

A[Hb],otai  =  -  (3.63/L)*  log  (A,/Ac)758 

+  (8.68/L)*  log  (At/Ac)782  (4) 


A[Hb02]  -  A[Hb]  =  -  (18.49/L)*  log  (A,/Ac) 

+  (2 1 ,20/L)  *  log  (Ai/Ac) 782  (5) 

where  A[j  represents  change  in  concentration.  The 
constants  were  computed  with  extinction  coefficients  for 
oxy-  and  deoxy-hemoglobin  at  the  two  wavelengths 
used. 

Once  stable  baseline  measurements  were 
achieved,  the  inhaled  gas  was  altered  to  pure  oxygen  or 
carbogen  and  dynamic  changes  were  observed  over  a 
period  of  two  hours.  Both  the  magnitude  and  rate  of 
change  of  s02  were  examined.  Following  the  NIR 
experiments,  the  MRI  experiments  were  performed. 

19F  MR-EPI 

All  MRI  experiments  were  performed  on  an 
Omega  CSI  4.7  T  40  cm  system  with  actively  shielded 
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gradients.  A  tunable  2  cm  single  turn  solenoid  coil 
was  placed  around  the  tumor  and  40  jal  HFB  were 
injected  directly  into  both  central  and  peripheral  regions 
of  the  tumor  using  a  32  G  needle.  Shimming  was 
performed  on  the  signal  of  the  tissue  water  to  a 
typical  linewidth  of  50  Hz*  3D  spin-echo  (SE)  *H  images 
were  acquired  for  anatomical  reference  and 
corresponding  l9F  images  were  then  obtained  to  show  the 
distribution  of  HFB  in  the  tumor.  Regional  tumor  p02 
maps  were  generated  using  19F  PBSR-EPI  based  on  the 
relationship:  /?02(torr)  =  [i?l~  0.074]/  0.0016,  where 
is  the  spin  lattice  relaxation  rate  of  HFB,  as  described  in 
detail  previously  [11].  Twenty-three  p02  maps  were 
produced  in  3  hours  with  respect  to  respiratory  challenge. 


RESULTS 


NIR  Results 

Figure  2  shows  the  time  course  of  changes  in 
tumor  vascular  hemoglobin  saturation  and  concentration 
accompanying  alterations  in  inhaled  gases  for  a  breast 
tumor  and  Figure  3  shows  the  result  for  a  prostate  tumor. 
Hemoglobin  saturation  and  concentration  are  presented 
as  relative  millimolar  changes.  Both  tumors  show 
significant  changes  in  vascular  oxygenation 
accompanying  respiratory  interventions.  Hemoglobin 
saturation  increased  almost  immediately  after  a  gas 
switch  from  baseline  (33%  02)  to  either  carbogen  or 
100%  02  and  increased  steadily  for  several  minutes,  and 
then  gradually  returned  to  baseline  after  the  gas  was 
switched  back  to  baseline.  In  contrast,  total  hemoglobin 
change  is  insignificant,  indicating  relatively  constant 
blood  volume  in  the  tumor* 


33%  Ch  carbogen  33%Q2  100%Q.  33%  Oz 


Figure  2.  Hemoglobin  saturation  and  concentration 
change  in  a  4.0  cm3  breast  tumor. 


33%  cartx^n  33%  100%  33% 


Figure  3*  Hemoglobin  saturation  and  concentration 
change  in  a  4.0  cm3  prostate  tumor. 


Figure  4.  *H  and  19F  coronal  images  of  a  breast  tumor.  FOV  «  48  x  48  mm,  matrix  size  » 
128  x  64,  and  slice  thickness  -  4  mm. 
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MRI  Results 

Figure  4  shows  conventional  spin  echo  (SE)  lH 
images  (upper)  and  corresponding  %  images  of  a 
representative  breast  tumor.  Figure  5a  shows  a  19F 
PBSR-EPI  projection  image  obtained  from  the  tumor 
shown  in  Figure  4  in  a  single  acquisition  (t  *  90  s)  and 
Figure  5b  shows  corresponding  pOt  map  (expanded). 


a) 


Figure  5,  (a)  A 19F  PBSR-EPI  projection  image  obtained 
from  the  tumor  shown  in  Figure  4  in  a  single  acquisition 
=  90  s).  Fourteen  images  were  acquired  with  variable 
relaxation  delays  (t)  ranging  from  200  ms  to  90  sec  and 
1.5  mm  in  plane  resolution,  (b)  Corresponding  pOi  map 
(expanded). 


The  ,9F  MR-EPI  oximetry  of  tumor  has  the 
distinct  advantage  over  other  techniques  that  subsequent 
measurements  are  completely  non-invasive.  The  greatest 
strength  of  this  method  is  the  ability  to  trace  the  fate  of 
individual  voxels  (regions)  with  respect  to  therapeutic 
interventions.  Figure  6  shows  dynamic  changes  in  p02  of 
six  specific  voxels  of  a  breast  tumor  with  respect  to 
different  inhaled  gases.  It  is  noteworthy  that  voxels  with 
high  baseline  p02  had  significantly  different  response 
characteristics  from  those  with  initially  low  p02,  which 
showed  small  changes.  Figure  7  shows  p02  histograms 
obtained  by  FREDOM  of  HFB  for  a  representative  breast 
tumor.  The  histograms  show  the  heterogeneity  of  p02 
values  within  the  tumor  as  well  as  the  mean  p02  values. 


Figure  6,  Dynamic  changes  in  p02of  six  specific  voxels 
of  a  breast  tumor. 


pO 2  (torr) 


Figure  7.  Histograms  of  oxygen  tension  for  an  NF13762 
breast  tumor  determined  by  FREDOM  of  HFB.  (A)  p02 
distribution  while  the  rat  breathed  33%  O2.  Mean  pCh 
-44  ±3  torr.  (b)  p02  distribution  while  the  rat  breathed 
catbogen  (95%  02  +  5%  C02).  Mean  p(h  =  99  ±  4  torr. 
(p<0.0001).  (c)  pOt  distribution  while  the  rat  breathed 
100%  02.  Mean p02  *  14$  ±  4  torr. 
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Comparison 

While  absolute  p02  values  are  important  for 
investigating  tumor  hypoxia,  dynamic  changes  may  be 
more  valuable  for  investigating  tumor  response  to 
therapeutic  interventions.  Figure  8  shows  comparisons 
between  the  dynamic  changes  in  p02  and  s02.  Since  the 
IQ  system  provides  a  global  s02  value,  each  p02  data 
point  is  represented  as  an  average  over  all  voxels  of  each 
p02  map. 


33%  carbogen  33% 


33%  carbogen  33% 


Figure  8.  Comparison  of  pOi  and  hemoglobin  saturation  (s02) 
of  a  breast  tumor,  (a)  s02  determined  by  NIR  spectroscopy. (b)  pCh 
determined  by  FREDOM.  Both  pOi  and  s02  increase  with  a  tran¬ 
sition  from  33%  to  carbogen  and  begin  to  decrease  when  switched 
back  to  33%  02.  Changes  in  vascular  s02  precede  tissue  pOi.  The 
tumor  size  ~  1 .9  x  2.2  x  2.0  cm. 

To  further  investigate  the  response,  we  modeled 
the  temporal  response  in  p02  and  s02  using  exponential 
equations: 

y  =  a  +  b  •  (1  -  e~t,T ) ,  for  increasing  trend  (6) 
y-a  +  b'e~tlx  ,  for  decreasing  trend  (7) 

where  y  is  p02  or  s02,  a  and  b  are  two  constants,  t  is 
time,  and  t  is  the  time  constant. 

Changes  in  s02  were  faster  than  p02  (Figure  9), 
especially  in  the  case  of  larger  tumor,  which  was  found 


to  be  less  well  oxygenated  and  presumably  less  well 
perfused. 


Gtos  ctengp 


Figure  9,  Comparison  of  p02  and  .v02  time  constants 
of  two  breast  tumors  (small  tumor  -  1 .9  x  2.2  x  2,0 
cm,  large  tumor  -  2.4  x  2.2  x  2. 1  cm).  The  temporal 
responses  in  p02  and  s02  were  modeled  using  expo¬ 
nential  equations. 


DISCUSSION 

Since  poorly  oxygenated  tumors  tend  to  resist 
conventional  therapy,  there  have  been  many  efforts  to  re¬ 
oxygenate  tumors  prior  to  therapy.  A  simple  intervention 
is  respiratory  challenge,  i.e.,  attempting  the  elevated 
oxygen  concentration  of  the  inhaled  gas.  Past  clinical 
trials  were  often  disappointing,  but  it  is  now  thought  that 
results  were  significantly  influenced  by  the  inability  to 
identify  hypoxic  tumors  (i.e.,  those  which  would  benefit 
from  manipulation),  a  priori.  As  techniques  become 
available  to  measure  tumor  oxygenation,  it  is  appropriate 
to  reevaluate  approaches  to  manipulating  tumor 
oxygenation.  By  increasing  the  oxygen  tension  of  the 
inspired  gas,  the  arterial  s02  should  increase,  leading  to 
increased  hemoglobin  saturation  of  the  tumor 
vasculature,  and  hence,  increased  tumor  tissue  p02.  Our 
data  indicate  that  breathing  elevated  02  did  indeed  have  a 
significant  effect  on  both  tumor  vascular  s02  and  tissue 
p02.  For  a  typical  large  breast  tumor,  significant  changes 
were  found  in  vascular  s02  (p<0.05)  and  in  tissue  p02 
(p<0.0001)  by  ANOVA  in  the  case  of  carbogen 
inhalation,  and  (p<0.01)  and  (pO.OOOl)  in  the  case  of 
oxygen  inhalation.  There  has  been  a  debate  as  to  whether 
carbogen  is  more  effective  at  modulating  tumor 
oxygenation  than  100%  oxygen  since  C02  is  a 
vasodilator.  Indeed,  recent  work  in  a  human  glioma 
xenografts  suggested  that  oxygen  alone  had  no  influence 
on  tumor  vascular  oxygenation,  whereas  carbogen 
produced  a  pronounced  effect  (12J.  Our  preliminary 
results  indicate  that  both  gases  produced  changes  in  s02 
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and  p02 ,  with  vascular  sOz  increasing  more  with  100% 
02  inhalation  than  with  carbogen. 

We  recently  showed  that  the  FREDOM 
approach  indicates  similar  p02  values  in  tumors  to 
electrodes  [13],  especially  on  interstitial  tissue  p02.  By 
contrast,  NIR  investigates  vascular  oxygenation.  As 
expected,  changes  in  vascular  oxygenation  (s02)  were 
found  to  be  more  rapid  than  tumor  tissue  p02i  with 
greater  differences  in  large  tumors.  This  probably 
reflects  the  extensive  perfusion  of  the  small  tumors  with 
lesser  perfusion  of  large  tumors  as  reflected  by  lower 
mean p02  and  larger  hypoxic  fraction  [13]. 

Regional  tumor  tissue  p02  and  blood  s02  are 
important  physiological  parameters.  The  capability  to 
measure  them  will  provide  insight  into  progressive 
physiological  changes  in  a  tumor  accompanying 
interventions.  NIR  has  the  advantages  of  being  entirely 
nonin vasive,  inexpensive,  portable,  and  real-time.  But 
the  MRI  approach  clearly  reveals  detailed  oxygenation 
heterogeneity.  The  correlation  of  NIR  and  MR 
technologies  will  provide  insight  into  issues  of  oxygen 
delivery  and  consumption.  We  believe  that  application  of 
multiple  approaches  to  tumor  oxygenation  can  lead  to 
better  understanding  of  tumor  physiology  and  probably 
optimized  tumor  therapy. 
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ABSTRACT 

Studies  have  shown  that  hypoxic  tumor  cells  are  relatively  more  resistant  to  radiotherapy,  chemotherapy,  and 
photodynamic  therapy.  Tumor  oximetry,  e.g.,  measurement  of  oxygen  tension  (pC>2)  of  tissue  and/or  blood 
oxygenation  (SO2)  of  the  vascular  bed,  could  be  valuable  for  optimizing  treatment  plans. 

In  this  study,  we  employed  a  recently  developed  homodyne  system  to  measure  changes  in  hemoglobin 
saturation  (SO2)  and  concentration  in  the  vascular  bed  of  rat  prostate  and  breast  tumors.  For  comparison,  tissue  p02 
values  were  measured  using  19F  MR  EPI  of  hexafluorobenzene,  providing  a  map  of  regional  tumor  oxygenation 
tension.  Both  S02  and  p02  measurements  were  taken  while  the  inhaled  gas  was  alternated  between  33%  oxygen, 
100%  oxygen  and  carbogen  (95%  oxygen,  5%  CO2). 

The  results  obtained  for  both  techniques  showed  significant  changes  in  tumor  oxygenation  accompanying 
respiratory  challenge,  with  changes  in  vascular  SO2  preceding  tissue  p02  change.  The  combined  use  of  these  two 
techniques  provides  new  insight  into  the  dynamics  of  tumor  oxygenation  by  making  available  a  method  of  obtaining 
regional  information  of  the  state  of  the  tissue,  as  well  as  a  non-invasive,  real-time  method  for  determining  changes 
in  the  vascular  bed. 


Keywords:  Frequency-Domain  Spectroscopy,  NIR  spectroscopy,  19F  MRI,  Hexafluorobenzene,  Oximetry 


1.  INTRODUCTION 

Frequently,  blood  vessel  formation  is  unable  to  keep  up  with  the  rapid  growth  of  a  tumor.  If  this  occurs,  the  cells  in 
the  tumor  furthest  from  a  fresh  blood  supply  will  suffer  a  lack  of  oxygen  and  hypoxic  areas  will  form  (chronic  hypoxia). 
These  regions  can  be  as  much  as  3  times  more  resistant  to  radiotherapy.  In  addition  to  studies  in  vitro  and  in  animal  tumors, 
there  is  increasing  evidence  from  clinical  trials  that  poorly  oxygenated  tumors  indicate  poor  prognosis  for  patients.2,3 
Methods  of  determining  the  oxygen  content  of  a  tumor  could,  therefore,  be  helpful  in  the  development  of  an  optimal 
treatment  plan.  This  paper  will  present  the  experimental  results  of  two  such  methods:  NIR  spectroscopy  to  determine  blood 
oxygenation  (S02)  of  the  tumor’s  vascular  bed  and  19F  MRI  of  hexafluorobenzene  (HFB)  to  determine  tissue  p02. 

NIR  spectroscopy,  through  use  of  a  recently  developed  frequency-domain  system,  based  on  an  in-phase  and 
quadrature  (IQ)  demodulator  chip4,  is  attractive  as  a  non-invasive,  inexpensive,  portable,  real-time  system  that  can  provide 
global  S02  values.  We  show  that  this  IQ  system  can  be  used  to  determine  the  S02  in  a  tumor’s  vascular  bed  and  measure  the 


response  of  blood  volume  and  oxygen  saturation  to  inhaled  gas.  The  technique  of  using  19F  MRI  relaxometry  to  map  tissue 
p02  is  also  relatively  new.5  The  spin-lattice  relaxation  rate  of  hexafluorobenzene  is  particularly  sensitive  to  oxygen  while 
being  insensitive  to  temperature.6  Following  direct  injection  of  HFB  into  a  tumor,  F  MRI  maps  tissue  pOz  at  millimeter 
resolution.  This  method  facilitates  measurements  of  dynamic  changes  in  p02  accompanying  therapeutic  interventions  and 
allows  the  fate  of  individual  voxels  to  be  traced. 

Through  comparison  of  these  two  techniques,  it  is  possible  to  examine  the  relationship  between  S02  of  the  vascular 
bed  and  p02  of  the  tissue.  Blood  oxygenation,  blood  volume,  arterial  S02  and  temperature  may  also  be  compared. 


2.  METHODS  AND  INSTRUMENTATION 


2.1  Tumor  Model 

Dunning  prostate  adenocarcinoma  R3327-AT1  was  implanted  in  adult  male  Copenhagen  rats  and  NF13762  breast 
tumor  in  female  Fisher  rats.  The  tumors  were  grown  in  pedicles  on  the  forebacks  of  the  rats  until  they  were  approximately  2 
cm  in  diameter.  Rats  were  anesthetized  with  200  |il  ketamine  hydrochloride  (100  mg/ml)  and  maintained  under  general 
gaseous  anesthesia  with  33  %  inhaled  02  [0.3  dm5/min  02,  0.6  drnVrnin  N20,  and  0.5%  methoxyflurane]  through  a  mask 
placed  over  the  mouth  and  nose.  Body  temperature  was  maintained  by  a  warm  water  blanket.  A  fiber  optic  pulse  oximeter 
was  placed  on  the  hind  foot  to  monitor  arterial  oxygenation  (AS02)  and  a  fiber  optic  probe  was  inserted  rectally  to  measure 
temperature.  Inhaled  gas  was  alternated  between  33%  oxygen,  100%  oxygen  and  carbogen  (95%  oxygen,  5%  carbon 
dioxide).  NIR  and  EPI  measurements  were  performed  sequentially  for  comparison. 

2.2  NIR  Spectroscopy 

As  shown  in  figure  1,  we  used  a  new  homodyne  system  able  to  determine  amplitude  and  phase  changes  of  light.4  In 
this  setup,  an  RF  source  modulates  the  light  from  two  laser  diodes  (wavelengths  758  nm  and  782  nm)  at  140  MHz.  The  light 
passes  through  fiber  optic  cables,  is  transmitted  through  the  tissue,  and  is  collected  by  a  second  fiber  bundle.  The  light  is 
then  detected  by  a  PMT  and  is  demodulated  with  a  commercially  available  in-phase  and  quadrature  (IQ)  demodulator  chip 
into  I  and  Q  components.  Once  these  components  are  put  through  a  low  pass  filter,  they  can  be  used  to  calculate  amplitude 
and  phase.  These  steps  can  be  seen  mathematically  in  equations  1-4. 

(1)  I(t)  =  2A  sin  (cot+0)  sin(cot)  =  A  cos  (0)  -  A  cos(cot+0)  — >  low  pass  filter  - >  Idc  =  A  cos(0) 

(2)  Q(t)  =  2A  sin  (cot+0)  cos(cot)  =  A  sin  (0)  +  A  sin(cot+0)  — >  low  pass  filter  — >  Qdc  =  A  sin(0) 

(3)  0  =  tankQdc/Idc) 

(4)  A=  (Ide'  +  Qdc2)1'2 

A  =  amplitude  of  detected  light;  0  =  phase  of  detected  light;  co=  modulation  frequency  (140  MHz) 


Figure  1:  Setup  for  NIR  experiment. 


NIR  spectroscopy  can  be  used  to  determine  hemoglobin  saturation  because  the  extinction  coefficient  values  of 
deoxygenated  hemoglobin  differ  from  those  of  oxygenated  hemoglobin  at  the  wavelengths  selected  (758  nm  and  782  nm). 
At  this  point  in  our  algorithm  calculations,  we  have  assumed  background  absorbance  to  be  negligible  and  estimated  that  the 
absorption  coefficients  were  composed  of  the  extinction  coefficients  for  deoxy-hemoglobin  and  oxy-hemoglobin  multiplied 
by  their  respective  concentrations  (equations  5&6). 

(5)  fia758  =  SHb758[Hb]  +  £Hb02758[Hb02] 

(6)  Pa78  =  £Hb782[Hb]  +  eHb02782[Hb02] 

The  IQ  system  does  give  both  phase  and  amplitude  values,  but  given  the  tumor’s  small  size  and  our  fiber 
configuration,  we  haven’t  yet  derived  a  suitable  algorithm  to  compute  Pa  and  p,,.  The  data  presented  in  this  paper  was 
analyzed  using  Beer-Lamberf  s  law  and  the  amplitude  values  to  find  trends  in  the  changing  absorption  coefficients  (equation 
7).  By  manipulating  equations  5-7,  algorithms  were  formed  which  allow  the  calculation  of  blood  volume  and  saturation 
trends  from  the  transmitted  amplitude  of  the  light  through  the  tumor  (equations  8&9). 

(7)  Pac-Itai=l/L*log(AI/Ac) 

(8)  A[Hb]totaI  =  -3.63*  log  (Aj/Ac)758  +  8.68*  log  (Aj/Ac)782 

(9)  A[Hb02]  -  A[Hb]  =  -18.49*  log  (Aj/Ac) 758  +  21.20  *  log  (Ai/Ac) 782 

Aj=  initial  amplitude  (amplitude  of  first  sample  taken);  Ac=  current  amplitude;  L  =  source/detector  separation;  The  constants 
were  computed  with  extinction  coefficients  for  oxy-  and  deoxy-  hemoglobin  at  the  two  wavelengths  used. 


2.3  MRI  Instrumentation  and  Procedure 

MRI  experiments  were  performed  on  an  Omega  CSI  4.7  T  40  cm  system  with  actively  shielded  gradients.  A 
homebuilt  tunable  2  cm  *H/19F  single  turn  solenoid  coil  was  placed  around  the  tumor.  The  first  step  of  the  MRI  procedure  is 
injection  of  hexafluorobenzene.  HFB  (40  pi)  was  administered  directly  into  the  tumor  using  a  fine  sharp  (32  G)  needle  with 
deliberate  dispersion  along  several  tracts  to  interrogate  both  central  and  peripheral  tumor  regions.  HFB  is  ideal  for  the 
imaging  of  p02  because  it  has  a  single  resonance  and  its  relaxation  rate  varies  linearly  with  oxygen  concentration.  *H  images 
were  acquired  for  anatomical  reference  using  a  traditional  3D  spin-echo  pulse  sequence  as  seen  in  figure  2a.  Conventional 
19F  MR  images  (figure  2b)  were  then  taken  to  show  the  3D  distribution  of  the  HFB  in  the  tumor.  Figure  2b  may  be  directly 
overlaid  over  figure  2a  to  show  the  position  of  the  HFB  in  that  slice. 


Figure  2a:  Proton  (*H)  coronal  image  of  a  representative 
slice  through  a  breast  tumor  (NF13762):  TR  =  250  ms, 
TE  =  8  ms,  NA  =  2,  FOV  =  48  x  48  mm,  slice 
thickness  =  4  mm,  and  matrix  size  =  128  x  64  x  8. 


Figure  2b:  Corresponding  19F  MR  Image  showing 
distribution  of  HFB  within  the  tumor: 

FOV  =  48  x  48  mm,  matrix  size  =  128  x  64 


Tumor  oxygenation  was  assessed  using  19F  PBSR-EPI  of  HFB.8  The  PBSR  preparation  pulse  sequence  consists  of  a 
series  of  20  non-spatially  selective  saturating  90  *  pulses  with  20  ms  spacing  to  saturate  the  ,9F  nuclei.  Following  a  variable 
delay  time  t  ,  a  single  spin  echo  EPI  sequence  with  “blipped”  phase  encoding  was  applied.9  A  PBSR-EPI  image 
corresponding  to  the  images  shown  in  figures  2a  and  2b  is  shown  in  figure  2c.  Fourteen  32  x  32  PBSR-EPI  images,  with  t 
ranging  from  200  ms  to  90  sec  and  an  FOV  of  40x40  mm,  were  acquired  in  eight  minutes.  An  R1  map  was  obtained  by 
fitting  signal  intensity  of  each  voxel  of  the  fourteen  images  to  a  three  parameter  relaxation  model  by  Levenberg-Marquardt 
least  squares  algorithm  (equation  10): 

(10) 

yn  (*-  J)  =  Ad,  j)  ■  [1  -  (1  +  W)  •  exp(-/?l(i,  j)  •  r„ )] 

(n  =  1,2, . .14) 

(ij  =  1,2, . ,32) 


where  yn  ( i ,  j)  is  the  measured  signal  intensity  corresponding  to  delay  time  tn  (the  nth  images)  for  voxel  A(i,j)  is  the 
fully  relaxed  signal  intensity  amplitude  of  voxel  (i,  j),  IF  is  a  dimensionless  scaling  factor  allowing  for  imperfect  signal 
conversion,  and  Rl(i,  j)  is  the  relaxation  rate  of  voxel  (i,j)  in  unit  of  sec"1 


Figure  2c:  l9F  PBSR-EPI  projection  image  obtained  from  the  tumor  in  Figure  I  in  a 
single  acquisition  (x  =  90  s).  Fourteen  images,  were  acquired  with  variable  relaxation 
delays  (x)  ranging  from  200  ms  to  90  sec.  Using  a  32  x  32  matrix,  FOV  of  ^x4§%nm, 
p02  maps  were  generated  with  1.25  x  1.25  mm  resolution. 


p02  maps  were  then  generated  by  applying  the  calibration  curve:  pC^Oorr)  =  [Rl(s  *)  -  0.074]/ 0.0016  to  the 
/?1  maps.10  The  map  shown  in  figure  2d  focuses  in  on  a  region  of  the  same  slice  that  was  presented  in  figures  2a-2c. 
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Figure  2d:  Typical  p02  map,  composed  of  fourteen  PBSR  images  from  the  tumor  presented  in  figure  2a-c.  Using  a  32  x  32  matrix,  FOV 
of  ^x4gmm,  p02  maps  were  generated  with  1.25  x  1.25  mm  resolution. 


3.  RESULTS 


3.1  NIR  Results 


The  effects  of  the  inhaled  gas  on  hemoglobin  saturation  and  concentration  as  recorded  by  the  IQ  system  are  shown 
below  in  figures  3a  &  3b,  .  The  X-axis  shows  time  in  minutes  from  the  start  of  the  experiment  and  the  dotted  vertical  lines 
mark  the  point  when  the  gas  was  changed.  Hemoglobin  saturation  and  concentration  are  presented  as  unit- ess  re  a  ive 
trends.  It  can  easily  be  seen  that  hemoglobin  saturation  begins  to  increase  almost  immediately  after  a  gas  sw 
baseline  (33%  oxygen)  to  either  carbogen  or  100%  oxygen  and  increases  steadily  for  several  minutes.  Tota  emogo  m 
change  is  quite  small  in  comparison,  indicating  relatively  constant  blood  volume  in  the  tumor.  These  tren  s  y 

consistent  for  both  breast  and  prostate  tumors.  Typical  responses  of  a  breast  tumor  and  prostate  tumor  are  presented  below. 


Figure  3b:  Hemoglobin  saturation  and  concentration  change  in  a  prostate  tumor. 


It  is  also  worthwhile  to  compare  the  saturation  changes  measured  in  the  tumors  vascular  bed  by  the  IQ  system  with 
the  arterial  saturation  changes  measured  by  a  pulse  oximeter  from  the  rat’s  hind  foot.  One  such  comparison  taken  from  a 
prostate  tumor  is  presented  below  as  figure  4.  Hemoglobin  saturation  in  the  vascular  bed  is  again  represented  as  a  unit-less 
trend,  arterial  saturation  values  are  presented  to  the  right.  Again,  the  X-axis  gives  the  time  from  the  beginning  of  the 
experiment  in  minutes  and  the  dotted  lines  mark  the  time  of  gas  change.  In  the  case  presented,  the  arterial  saturation  follows 
the  same  trend  as  the  tumor  vascular  bed’s  hemoglobin  saturation,  but  shows  a  faster  change.  It  should  be  noted  that  not  all 
the  rats  we’ve  studied  have  behaved  this  way.  We  saw  a  great  amount  of  variability. 


Figure  4:  In  the  prostate  tumor  presented  in  Figure  4,  both  arterial  S02  and  the  S02  in  the  tumor  increased  for  inhaled  gases  carbogen  and 
100%  02  and  decreased  for  33%  02.  Arterial  S02  in  hind  foot  measured  by  commercial  pulse  oximeter  and  S02  in  tumor  using  IQ 
system. 


3.2  MR1  Results 

MRI  provides  the  advantage  of  being  able  to  look  at  regional  changes  in  p02  values.  Histograms,  such  as  those 
presented  in  figure  5,  are  able  to  show  the  heterogeneity  of  p02  values  within  the  tumor  as  well  as  the  average  p02  values. 
The  data  presented  here  were  taken  from  a  breast  tumor  and  show  the  average  values  from  several  p02  maps.  In  figure  5a, 
we  see  that  when  the  rat  was  breathing  33%  oxygen,  the  average  p02  value  was  about  40  torr.  When  the  rat  was  breathing 
carbogen  (figure  5b),  there  was  a  large  shift  towards  higher  p02  values  leading  to  a  mean  value  of  about  99  torr.  These 
values  increased  further  while  the  rat  was  breathing  100%  oxygen  (figure  5c)  such  that  the  average  voxel  now  had  a  p02 
value  of  about  145  torr.  The  time  course  of  these  changes  will  be  presented  later  in  the  paper  in  figure  7b. 
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Figure  5a:  percent  voxels  in  a  certain 
p02  range  while  the  rat  was  breathing 
33%  oxygen. 

Average  value  =  39.96  ±  3.09  torr 


Figure  5b:  percent  voxels  in  a  certain 
p02  range  while  the  rat  was  breathing 
carbogen. 

Average  value  =  99.30  ±  3.70  torr 


Figure  5c:  percent  voxels  in  a  certain 
p02  range  while  the  rat  was  breathing 
100%  oxygen. 

Average  value  =  145.17  ±  3.52  torr 


p 02  (torr) 


3.3  Comparison 

While  absolute  p02  values  are  important  for  investigating  hypoxia,  dynamic  changes  may  be  more  interesting  for 
investigation  of  response  to  intervention.  In  Figure  6  &  7,  the  dynamic  changes  in  p02  and  S02  are  compared.  These  plots 
show  that  p02  reacts  in  a  very  similar  fashion  to  blood  saturation,  but  that  the  effect  is  slower.  Since  the  IQ  system  provies  a 
global  measure  across  the  whole  tumor,  p02  measurements  are  presented  as  mean  values  attained  from  the  pixels  of  each 
p02  map.  The  NIR  and  MRI  data  shown  are  for  the  same  rat  undergoing  the  same  procedure,  but  on  consecutive  days.  Data 
from  two  breast  tumors  of  various  size  are  presented.  The  larger  tumor  was  about  2.1  x  2.4  x  2.2  cm,  the  smaller  tumor  was 
1.9  x  2.2  x  2  cm.  Carbogen  and  100%  oxygen  each  consistently  showed  increased  S02  and  p02  over  33%  oxygen.  The  p02 
and  S02  measurements  were  plotted  against  the  same  X-axis  to  allow  for  rate  of  change  comparisons. 

Figure  6a:  Saturation  from  IQ 


Figure  6b:  p02  from  MRI 


Figure  10:  These  two  sets  of  data  were  taken  from  the  small  rat  tumor  on  consecutive  days.  Both  the  p02  and  blood  saturation  increased 
with  a  transition  from  33%  to  carbogen,  and  decreased  when  switched  back  to  33%. 


Figure  7a:  Saturation  from  IQ 


Inhaled 

gas  -->  33%  carbogen  100%  33% 


Figure  7b:  p02fromMRI 


Figure  7:  These  two  sets  of  data  were  taken  from  the  large  rat  tumor  on  consecutive  days.  Both  the  p02  and  blood  saturation  increased 
with  a  transition  from  33%  to  carbogen,  and  then  further  increased  slightly  when  switched  from  carbogen  to  100%  oxygen. 
Values  began  to  decrease  when  switched  back  to  33%. 


To  further  study  the  temporal  response,  the  changes  in  p02  and  S02  were  modeled  using  the  exponential  equations  11  and  12 
to  provide  time  constant  T. 

(11)  For  increasing  values:  Saturation  =  a  +  b*(l-e_t/T) 

(12)  For  decreasing  values:  Saturation  =  a  +  b*(e't/x) 

Generally,  blood  saturation  effects  had  a  much  shorter  time  constant  than  oxygen  tension  in  the  tissue  (Figure  8).  For  the 
smaller  tumor,  the  rate  of  increase  and  decrease  were  much  faster  for  S02  than  p02.  Less  difference  was  seen  in  the  larger 
tumor. 


33%02~  1 00%02  --  33%02  -  Carbogen -- 

>Carbogen  >33%02  >Carbogen  >33%02 

gas  change 


Figure  8:  Time  constants  for  the  data  presented  in  Figures  6&7. 


4.  CONCLUSION 

The  data  suggest  interesting  correlations  between  several  vascular  parameters.  Both  tumor  vascular  Hb  saturation 
and  mean  p02  increased  in  response  to  inhaling  an  elevated  percent  02,  either  through  carbogen  or  100%  02.  Arterial  S02 
and  tumor  S02  respond  similarly  to  changes  in  inhaled  gas,  with  arterial  changes  preceding  changes  in  the  tumor  vascular 
bed.  Changes  in  S02  are  generally  faster  than  p02,  though  absolute  values  are  highly  variable  and  suggest  heterogeneity 
amongst  the  tumor  population.  Given  the  distinct  heterogeneity  among  tumors  even  of  a  given  type  and  size,11  further 
investigations  are  required  to  form  a  sound  picture  of  the  interplay  of  multiple  vascular  parameters. 
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ABSTRACT 


We  recently  described  a  novel  approach  to  measuring  regional  tumor  oxygen  tension 
using  pulse  burst  saturation  recovery  NMR  echo  planar  imaging  relaxometry  of 
hexafluorobenzene.  We  have  now  compared  oxygen  tension  measurements  in  a  group  of 
size  matched  Dunning  prostate  rat  tumors  R3 327- ATI  made  using  this  FREDOM 
(Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping), 
approach  with  a  traditional  polarographic  method:  the  Eppendorf  Histograph.  We  also 
compare  MR  and  electrode  approaches  to  monitoring  dynamic  changes  with  respect  to 
interventions  and  demonstrate  extension  of  the  MR  technique  to  rat  breast  tumors. 


Key  words:  echo  planar  imaging,  electrode,  MRI,  oxygen,  prostate,  tumor 

Abbreviations  ARDVARC  (Alternated  Relaxation  Delays  with  Variable  Acquisitions  to 
Reduce  Clearance  effects);  EPI  (echo  planar  imaging);  FREDOM  (Fluorocarbon 
Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping);  HFB 
(hexafluorobenzene);  i.t  (intra  tumoral) 

Acknowledgments 

This  work  was  supported  in  part  by  The  American  Cancer  Society  (RPG-97-116- 
010CCE;  RPM),  DOD  Breast  Cancer  initiative  (YS),  Verein  zur  Forderung  der 
Krebserkennung  and  Krebshandlung  e.V.  Heidelberg  (PP)  and  the  NIH  BRTP  Facility 
#5  -P4 1 -RR025 84 . 


2 


/ 


Introduction 

It  is  widely  appreciated  that  tumor  oxygenation  may  significantly  influence 
therapeutic  success.  In  particular,  the  efficacy  of  radiotherapy  [1],  photodynamic  therapy 
[2]  and  hypoxia  selective  chemotherapeutic  agents  [3]  depend  on  p02-  It  had  been 
suggested  that  the  ability  to  measure  tumor  oxygenation  in  patients  could  allow  therapy  to 
be  individualized  and  optimized  [4],  and  indeed,  several  recent  studies  have  found 
significant  prognostic  value  based  on  the  Eppendorf  Histograph  in  assessing  clinical 
tumors  [5-7].  While  electrodes  may  be  considered  a  "gold  standard",  they  have  certain 
shortcomings  and  there  is  clearly  a  need  for  alternative  methods  [8].  We  have  been 
developing  a  new  approach  based  on  19F  NMR  of  perfluorocarbons  [9-12]  and  believe 
the  method  can  now  provide  useful  measurements  of  tumor  oxygen  dynamics  in  vivo. 

The  FREDOM  approach  exploits  the  exceptional  response  of  the  ^9F  NMR  spin 
lattice  relaxation  rate,  Rl,  of  fluorocarbons  to  changes  in  oxygen  tension.  Fluorocarbons 
act  as  ideal  liquids  and  solvation  of  gases  is  directly  proportional  to  the  partial  pressure  of 
the  gas  (Henry'  law).  Since  oxygen  (02)  is  paramagnetic  it  induces  relaxation  in  solution 
directly  proportional  to  the  concentration  of  oxygen,  and  hence,  p02  [13].  The  highly 
hydrophobic  nature  of  fluorocarbons  ensures  both  a  high  solubility  of  gases,  providing 
molecular  amplification,  and  minimal  solvation  of  other  materials  (e.g.,  metal  ions) 
minimizing  interference  from  other  environmental  factors.  We,  and  others,  have  explored 
the  use  of  numerous  PFC  reporter  molecules  and  various  routes  of  administration  [14]. 
We  believe  that  direct  intra  tumoral  (i.t.)  injection  of  HFB  provides  an  optimal  approach 
to  tumor  oximetry,  and  should  provide  measurements  comparable  to  those  obtained  using 
electrodes.  In  addition,  this  minimally  invasive  approach  facilitates  mapping  of  dynamic 
changes  in  p02  with  respect  to  interventions. 

Methods 

Dunning  prostate  adenocarcinomas  (R3327-AT1)  were  implanted  in  male 
Copenhagen  rats  (-250  g),  as  described  in  detail  previously  [15].  Tumors  were  divided 
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into  two  groups  and  allowed  to  grow  to  about  ~2  cm3  or  >  3.5  cm3  volume.  For  MR 
investigations  rats  were  placed  under  general  gaseous  anesthesia  with  33%  inhaled  O2 
(0.3  dm3/min  O2, 0.6  dm3/min  N2O,  and  0.5%  methoxyflurane.  Hexafluorobenzene  (25  - 
40  pi)  was  injected  directly  into  the  tumors  in  both  central  and  peripheral  regions  using  a 
Hamilton  syringe  with  a  custom  made  fine  sharp  needle  (32  G).  A  fiber  optic  probe  was 
placed  rectally  to  monitor  core  temperature.  NMR  experiments  were  performed  using  an 
Omega  CSI  4.7  Tesla  horizontal  bore  magnet  system  with  actively  shielded  gradients 
with  a  tunable  (^^F)  single  turn  size-matched  solenoid  coil  placed  around  the  tumor. 
Following  traditional  imaging  to  establish  the  distribution  of  HFB,  tumor  oxygenation 
was  estimated  on  the  basis  of  19F  PBSR  EPI  relaxometry  of  the  HFB  [10]  with  a  typical 
1.25  mm  in  plane  resolution.  For  initial  work  three  consecutive  R1  measurements  were 
made  over  a  period  of  1  hour  to  investigate  reproducibility,  and  stability  of  the  system 
when  the  rats  breathed  33%  O2  (baseline).  Since  R1  is  a  linear  function  of  pC>2  at 
constant  temperature,  pC>2  was  estimated  on  a  voxel  by  voxel  basis  using  the  relationship 
pC>2  (torr)  =  (R1  -  0.074)/0.0016  [10].  The  inhaled  gas  was  then  altered  to  100%  oxygen, 
and  relaxation  measurements  (three)  were  immediately  repeated  over  a  period  of  1  hour. 
Finally,  the  gas  was  switched  back  to  the  baseline  state  and  three  further  pC>2 
determinations  were  immediately  performed  over  1  hour.  Our  initial  studies  required  20 
mins  to  produce  a  p02  map,  but  more  recent  introduction  of  the  ARDVARC  acquisition 
protocol  [12]  provides  enhanced  maps  in  8  mins.  Breast  13762  NF  adenocarcinomas  were 
examined  similarly. 

Histography  was  applied  to  groups  of  size  matched  tumors,  which  did  not  receive 
HFB.  Halothane  was  used  in  place  of  methoxyflurane.  Using  the  Eppendorf  Histograph 
100  to  200  individual  pC>2  determinations  were  made  in  each  tumor  ,as  recommended  by 
the  manufacturer.  For  dynamic  measurements  a  Diamond  General  micro-electrode  (700 
pm)  was  inserted  to  a  specific  location.  Baseline  p02  was  measured  and  the  inhaled  gas 
altered  to  100%  O2  or  carbogen  (95%02/5%C02)  for  30  mins.  At  this  stage  p02  was 
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again  measured.  Following  a  series  of  measurements  with  different  gases  at  one  location, 
the  needle  was  moved  and  the  gases  cycled  again. 

Statistical  significance  of  changes  in  oxygenation  was  assessed  using  analysis  of 
variance  (ANOVA)  on  the  basis  of  Fisher  PLSD.  Experiments  were  approved  by  the 
Institutional  Animal  Care  and  Advisory  Committee  conducted  in  accordance  with 
National  Laws. 

Results 

Both  FREDOM  and  electrode  methods  indicated  similar  oxygen  tension  distributions 
for  the  ATI  tumors  (Fig.  1).  Moreover,  both  techniques  showed  that  tumors  with  volume 
>3.5  cm3  were  significantly  (p  <  0.0001)  less  well  oxygenated  than  smaller  tumors 
(volume  <  2  cm3).  For  the  large  tumors  FREDOM  indicated  median  p02  =  2  torr  and 
fraction  <  10  torr  (HFio)  =  82  %,  while  the  Eppendorf  electrode  indicated  median  p02  = 
3  torr  and  HFio  =  84%.  For  the  small  tumors  the  match  was  less  good  with  median  =  15 
v  8  torr  and  HFio  =  44  versus  66%  for  FREDOM  and  electrode,  respectively. 
Examination  of  the  MR  images  showed  that  for  1  small  tumor  most  of  the  HFB  resided 
very  close  to  the  tumor  edge  and  may  have  biased  the  apparent  p02-  Indeed,  if  this  tumor 
was  excluded  there  was  no  significant  difference  between  the  respective  p02 

distributions. 

Using  the  FREDOM  approach  we  also  examined  response  to  respiratory  challenge. 
Increasing  the  concentration  of  inspired  oxygen  from  33%  to  100%  02  produced  a 
significant  increase  (p  <  0.0001)  in  tumor  oxygenation  for  the  group  of  small  tumors.  In 
contrast  no  change  was  observed  in  the  mean  p02  for  the  group  of  large  tumors.  A 
strength  of  the  FREDOM  approach  is  the  ability  to  follow  individual  tumor  regions,  with 
respect  to  intervention,  in  this  case  respiratory  challenge.  Six  representative  regions  were 
selected  from  a  single  tumor  (Fig.  2a).  Three  regions,  which  were  initially  well 
oxygenated  (p02  >  10  torr)  showed  rapid  and  significant  increases  within  8  minutes  of 
switching  from  33%  02  to  100%  02-  Changes  in  relatively  poorly  oxygenated  regions 
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were  much  slower,  although  2  of  3  regions  did  show  a  significant  change  in  p02  after  24 
mins. 

Electrode  investigation  of  dynamic  changes  in  pC>2  also  showed  3  of  six  regions  with 
significant  changes  in  switching  from  33%  02  to  100%  02,  but  only  1  region  was  also 

significantly  different  with  carbogen  (Fig.  2b). 

In  a  representative  large  breast  tumor  (~  4  cm3)  we  found  significant  changes  in  p02 
(p  <  0.0001)  with  respect  to  respiratory  challenge  with  baseline  mean  p02  =  40+3  rising 
to  mean  p02  =  99±4,  when  rat  inhaled  carbogen  and  mean  p02  =  145  ±4  for  oxygen 

inhalation. 

Discussion 

These  results  demonstrate  the  similarly  of  measurements  obtained  using  traditional 
electrodes  or  the  new  FREDOM  approach  to  tumor  oximetry.  In  each  case  there  was  a 
significant  difference  in  pC>2  observed  in  small  versus  large  ATI  tumors.  For  larger 
tumors  the  hypoxic  fraction,  mean  and  median  were  very  similar,  together  with  the  range 
of  typical  pC>2  values.  In  smaller  tumors  MR  suggested  a  larger  range  with  a  number  of 
measurements  in  excess  of  100  torr.  This  may  have  arisen  from  measurements  close  to 
the  tumor  periphery,  which  are  less  common  using  electrodes. 

A  significant  strength  of  the  FREDOM  approach  is  the  ability  to  monitor  dynamic 
changes  in  regional  p02  in  response  to  acute  interventions.  Others  have  used  the 
Eppendorf  system  to  examine  acute  chnages  [16],  but  this  required  reintroduction  of  the 
needle  electrode  and  generation  of  new  tracks.  Not  only  was  this  invasive,  but  it  also  led 
to  sampling  of  parallel  tissue  regions  rather  than  the  fate  of  specific  regions.  Given  the 
extensive  heterogeneity  encountered  in  tumors  and  steep  local  gradients  in  p02  we 
believe  it  will  be  valuable  to  follow  individual  tumor  regions.  Historically,  regional 
response  to  intervention  was  assessed  by  placing  an  electrode  at  a  specific  location  and 
monitoring  changes  in  p02  [17].  e  have  now  performed  such  experiments  with  a  micro 
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electrode  and  found  a  range  of  baseline  p02  values  and  response  to  respiratory  challenge 
similar  to  those  seen  using  MRI. 

We  have  now  shown  both  that  there  is  distinct  intra  tumoral  heterogeneity  in  baseline 
oxygenation  in  the  Dunning  prostate  ATI  tumor  and  also  in  the  response  to  intervention. 
In  common  with  our  previous  observations  a  three  fold  change  in  FO2  seems  to  lead  to  a 
threefold  response  in  tumor  p02-  However,  the  rate  of  change  in  highly  variable. 
Preliminary  data  with  8  min  time  resolution  suggest  that  well  oxygenated  regions  respond 
rapidly,  whereas  those  poorly  oxygenated  require  much  longer.  Such  observations  could 
have  significant  implications  for  patient  inhalation  times  prior  to  therapy:  while  previous 
work  had  shown  that  Pre  Irradiation  Breathing  Times  (PIBT)  could  substantially 
influence  the  effect  of  oxygen  or  carbogen  breathing  [18],  the  differential  response  of 
individual  tumor  regions  may  not  have  been  fully  appreciated. 

In  developing  a  new  technique  it  is  important  demonstrate  its  reliability,  robustness 
and  general  application.  We  and  several  other  groups  have  now  applied  the  FREDOM 
approach  to  tumor  oximetry.  Initially  investigators  favored  intra  venous  or  intra 
peritoneal  administration  of  emulsions  of  fluorocarbons.  While  material  became  trapped 
in  tumors  and  could  be  used  to  report  pC>2  [19-22],  it  became  increasingly  apparent  that 
material  delivered  via  the  vasculature  tended  to  bias  measurements  towards  well  perfused 
tumor  regions  [22].  Indeed,  recent  measurements  by  Griffiths  et  al.  have  confirmed  such 
a  bias  [23].  Furthermore,  the  use  of  emulsions  to  carry  the  PFCs  tend  to  lead  to  extensive 
uptake  by  the  reticuloendothelial  system  with  hepatomegaly.  Intra  tumoral  administration 
is  minimally  invasive  provided  that  a  fine  sharp  needle  is  applied,  as  we  have  used  here. 
We  have  now  extended  our  work  from  the  Dunning  prostate  R3327-AT1  tumor,  which  is 
poorly  differentiated,  has  only  microscopic  necrotic  foci  and  is  firm,  to  the  13762  breast 
tumor,  which  has  less  structure  and  considerable  cystic  fluid.  Here,  we  have  simply 
reported  the  ability  to  measure  dynamic  changes  in  the  breast  tumor  oxygenation,  but  in 
the  accompanying  work  (Song  et  al.,  this  volume),  we  show  more  extensive  results. 
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Since  the  MR  and  electrode  approaches  appear  to  give  similar  results  one  may  debate 
the  relative  their  merits.  Clearly,  MR  is  very  expensive,  with  a  typical  imaging  system 
costing  upwards  of  $1  M,  compared  with  $60  000  for  the  Eppendorf  and  <  $5  000  for  a 
laboratory  micro  electrode  system.  However,  MR  facilitates  the  simultaneous 
measurements  of  dynamic  changes  in  response  to  intervention  at  multiple  points  within  a 
tumor.  While  we  were  able  to  follow  changes  in  pC>2  at  specific  regions  using  a  needle 
electrode  with  placement  at  sequential  locations  accompanied  by  cycling  of  the 
intervention,  such  an  approach  would  be  less  satisfactory  for  other  interventions,  and 
even  here,  may  have  led  to  some  conditioning  or  hysteresis.  The  FREDOM  approach  may 
be  readily  combined  with  other  measurements  such  as  blood  flow/perfusion  [24],  pH  [25] 
or  metal  ions  by  infusion  of  appropriate  reporter  molecules  [26]. 

As  a  reporter  molecule  HFB  has  many  advantageous  properties.  It  is  cheap,  readily 
available,  and  exhibits  minimal  acute  toxicity  (LD50  >  25  g/kg)  [27].  No  signs  of  renal  or 
hepatic  toxicity  have  been  found  [28]  and  others  have  tested  doses  as  high  as  50  g/kg 
(twice  weekly)  orally  in  rats  over  35  weeks  [29].  We  typically  find  substantial  clearance 
from  tumors  within  24  h,  though  this  does  limit  our  measurements  to  acute  response  to 
interventions  [12].  High  symmetry  within  the  molecule  leads  to  a  single  19F  MR 
resonance  providing  optimal  SNR.  The  R1  (=1/T1)  is  highly  sensitive  to  p02  while 
showing  little  response  to  temperature  [9].  Long  Tls  up  to  14  s  appear  to  make  HFB  less 
efficient  for  spin  lattice  relaxometry,  but  use  of  the  pulse  burst  saturation  recovery 
approach  minimizes  the  length  of  th  experiment  [10]  and  a  large  range  of  T1  values  is  a 
requisite  for  sensitivity  to  changes  in  pC>2-  The  long  transverse  relaxation  time  (T2)  is 

ideally  suited  to  echo  planar  imaging. 

The  ultimate  value  of  a  novel  technique  will  depend  on  its  adoption  by  multiple 
laboratories,  and  the  significance  of  the  results  that  can  be  generated.  We  believe  that  the 
FREDOM  approach  is  versatile  and  we  are  demonstrating  increasing  applications,  and 
thus,  we  foresee  expanded  future  application  of  the  technique. 
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Figure  legends 
Figure  1 

Comparison  of  oxygenation  in  size-matched  groups  of  ATI  tumors  based  on  MR  EPI 
relaxometry  (left)  and  electrode  polarography  (right),  when  rats  inhaled  33%  02-  Small 
tumors  are  shown  at  top  (volume  <  2  cm^)  and  large  tumors  below  (volume  >  3.5  cm^). 
Each  method  shows  a  significant  difference  in  tumor  oxygenation  for  small  versus  large 
tumors  (p<  0.0001). 

Figure  2 

a)  Dynamic  changes  in  pC>2  of  six  specific  regions  of  an  ATI  tumor.  The  three  high  p02 
regions  had  significantly  different  pC>2  (*  p<  0.05) )  from  those  with  low  pC>2  at  each 
time  point.  Within  8  mins  of  elevating  inspired  02  the  three  high  pC>2  voxels  had 
significantly  increased  p02  (p<  0.05)  while  the  low  p02  voxels  required  >  24  mins  to 
show  significant  changes.  All  six  regions  were  observed  simultaneously  using  the 
FREDOM  approach. 

b)  Dynamic  changes  in  pC>2  of  six  specific  regions  of  an  ATI  tumor.  The  electrode  was 
placed  in  one  location  at  a  time  and  inhaled  gases  cycled  for  subsequent  locations. 
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Regional  Tumor  Oxygen  Tension  and  Blood  Flow: 
Correlation  Studies  Using  WF  PBSR-EPI  of  Hexafluorobenzene 


Y.  Song,  R.  P.  Mason,  S.  Hunjan,  A.  Constantinescu,  E.  W.  Hahn,  and  P.  P.  Antich, 
Department  of  Radiology,  UT  Southwestern  Medical  Center,  Dallas,  Texas,  USA 

Introduction:  It  is  recognized  that  therapeutic  efficacy  may  be  influenced  by  tumor 
oxygenation.  In  particular,  hypoxic  tumors  resist  radiotherapy.  We  have  recently  shown  the 
feasibility  of  monitoring  tumor  oxygen  tension  based  on  19F  PBSR-EPI  of  hexafluorobenzene 
(HFB)  [1].  We  also  found  that  HFB  clears  from  tumors  over  a  period  of  hours  [2].  Since  HFB  is 
a  non-ionic  freely  diffusable  tracer,  it  appeared  that  clearance  rate  would  provide  an  indication  of 
relative  tumor  blood  flow  (TBF).  We  have  now  investigated  the  feasibility  of  mapping  the 
clearance  rate  of  HFB  and  correlating  this  putative  blood  flow  marker  with  corresponding  p02. 

Methods:  Dunning  prostate  R3327-AT1  or  breast  13762  NF  adenocarcinoma  was  implanted  in 
a  skin  pedicle  on  the  foreback  of  a  rat.  When  the  tumor  reached  1-2  cm  diameter,  40  |il  HFB 
were  injected  directly  into  the  tumor  (IT),  both  centrally  and  peripherally.  The  rat  was 
maintained  under  general  gaseous  anesthesia  (33%  O2,  66%  N2O  and  0.5%  methoxyflurane).  A 
homebuilt  tunable  2  cm  *H/19F  single  turn  solenoid  coil  was  placed  around  the  tumor  and  MR 
experiments  were  performed  using  a  4.7  T  magnet  equipped  with  actively  shielded  gradients.  3D 
*H  images  were  acquired  for  anatomical  reference  and  corresponding  19F  images  were  obtained 
to  show  the  distribution  of  HFB.  Tumor  oxygenation  was  assessed  using  19F  PBSR-EPI  of  HFB. 
By  applying  the  acquisition  protocol  ARDVARC  (Alternated  Relaxation  Delays  with  Variable 
Acquisitions  to  Reduce  Clearance  effects)  [2],  we  achieved  /?1  maps  in  8  min.  A  series  of  maps 
were  acquired  over  a  period  of  2  hours  with  respect  to  respiratory  challenges.  pOj  maps  were 
then  generated  by  applying  the  relationship:  /702(torr)  =[/?l(.y_1)-0.074]/0.0016  to  the  R 1 
maps.  The  data  also  allowed  us  to  produce  a  clearance  map  based  on  EPI  images  with  the  longest 
delay  (90  s). 

Results:  pOi  maps  were  generated  with  a  typical  precision  of  2  -  5  torr  and  30  -  100  individual 
voxels  within  a  tumor.  For  many  regions,  the  HFB  signal  intensity  was  found  to  decline 
exponentially  with  a  typical  clearance  half-life  ranging  from  T\a  —  700  to  1200  min,  though 
many  voxels  indicated  no  apparent  changes. 

Discussion:  Regional  tumor  oxygen  tension  and  blood  flow  are  important  physiological 
parameters  and  the  opportunity  to  measure  both  simultaneously  would  be  of  value  in 
physiological  research.  Based  on  the  preliminary  data  presented  here,  we  believe  that  clearance 
of  HFB  provides  an  indication  of  relative  TBF  by  analogy  with  studies  of  cerebral  blood  flow 
using  freon-23  [3].  In  future  studies,  such  measurements  will  be  rigorously  evaluated. 
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TUMOR  OXIMETRY:  COMPARISON  OF  19F  MR  EPI  AND  ELECTRODES 


Mason,  Ralph;  Hunjan,  Sandeep;  Constantinescu,  Anca;  Song,  Yulin;  Zhao,  Dawen;  Hahn,  Eric;  Antich, 

Peter;  and  Peschke,  Peter+.  U.T.  Southwestern  Medical  Center,  Dallas,  TX  and  +DKFZ,  Heidelberg, 
Germany 

Therapeutic  efficacy  may  be  influenced  by  tumor  oxygenation.  In  particular,  hypoxic  tumors  resist 
radiotherapy  and  may  be  good  candidates  for  hypoxia  selective  cytotoxic  agents.  We  recently  described  a 

novel  approach  to  measuring  regional  tumor  oxygen  tension  using  pulse  burst  saturation  recovery 
(PBSR)  nuclear  magnetic  resonance  (NMR)  echo  planar  imaging  (EPI)  relaxometiy  of  hexafluorobenzene 
(HFB)  (1).  We  have  now  compared  oxygen  tension  measurements  in  a  group  of  size  matched  Dunning 
prostate  rat  tumors  R3327-AT1  made  using  this  new  method  with  a  traditional  polarographic  method:  the 
Eppendorf  Histograph.  We  also  demonstrate  extension  of  the  MR  technique  to  rat  breast  tumors. 

Methods:  Dunning  prostate  R3327-AT1  or  breast  13762  NF  adenocarcinoma  was  implanted  in  a  skin 

pedicle  on  the  foreback  of  a  rat.  When  the  tumor  reached  a  volume  ~2  cm^  or  >  3.5  cm^,  40  pi  HFB  were 
injected  directly  into  both  central  and  peripheral  regions  of  the  tumor.  The  rat  was  maintained  under 

general  anesthesia  (33%  C>2,  66%  N2O  and  0.5%  methoxyflurane).  A  tunable  2  cm  ^H/l^F  single  turn 
solenoid  coil  was  placed  around  the  tumor  and  MR  experiments  were  performed  using  a  4.7  T  magnet 
equipped  with  actively  shielded  gradients.  3D  *H  images  were  acquired  for  anatomical  reference  and 
corresponding  images  were  obtained  to  show  the  distribution  of  HFB.  Tumor  oxygenation  was 
assessed  using  l^F  PBSR-EPI  of  HFB.  A  series  of  p02  maps  was  acquired  over  a  period  of  2  hours  with 
respect  to  respiratory  challenge  using  the  relationship:  pC>2  (torr)  =  (Rl-0.074)/0.0016.  Histography  was 
applied  to  groups  of  size  matched  tumors,  which  did  not  receive  HFB. 

Results:  Similar  oxygen  tension  distributions  were  found  using  each  method  and  both  techniques 

showed  that  tumors  with  volume  >  3.5  cm^  were  significantly  (p  <  0.0001)  less  well  oxygenated  than 

smaller  tumors  (volume  <  2  cm^).  Using  the  l^F  EPI  approach  we  also  examined  response  to  respiratory 
challenge.  Increasing  the  concentration  of  inspired  oxygen  from  33%  to  100%  O2  produced  a  significant 
increase  (p  <  0.0001)  in  tumor  oxygenation  for  a  group  of  small  tumors.  In  contrast  no  change  was 
observed  in  the  mean  pC>2  for  a  group  of  large  tumors.  Consideration  of  individual  tumor  regions, 
irrespective  of  tumor  size  showed  a  strong  correlation  between  the  maximum  pC>2  observed  when 
breathing  100%  O2,  as  compared  with  mean  baseline  p02- 

Discussion:  These  results  demonstrate  the  similarity  of  results  obtained  using  electrode  or  MR 

approaches  to  tumor  oximetry.  They  also  indicate  the  feasibility  of  l^F  MR  to  monitor  dynamic  changes 
in  regional  p02  in  response  to  acute  interventions.  The  ability  to  measure  pC>2  could  be  valuable  in  pre- 
clinical  evaluation  of  novel  therapies  and  could  allow  therapy  to  be  individualized  and  optimized  for 
patients. 
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TUMOR  OXYGEN  DYNAMICS:  COMPARISON  BETWEEN  19F  MR-EPI  OF 
HEXAFLUOROBENZENE  AND  FREQUENCY  DOMAIN  NIR  SPECTROSCOPY 


Song1,2,  Y.;  Worden1,  K.  L.;  Jiang1,  X.;  Zhao2,  D.;  Constantinescu2,  A.;  Liu1,  H.;  and  Mason2,  R.  P. 
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Introduction:  Oxygen  plays  a  key  role  in  tumor  therapy  and  may  be  related  to  tumor  development: 
e.g.,  angiogenesis  and  metastasis.  Using  noninvasive  techniques  to  accurately  measure  oxygenation 
could  assist  in  developing  novel  therapies.  Here,  we  have  used  19F  MR-EPI  relaxometry  of 
hexafluorobenzene  (HFB)[1]  to  monitor  tissue  oxygen  tension  (p02)  of  rat  breast  tumors  and  compared 
the  results  with  changes  in  hemoglobin  saturation  (SO2)  and  concentration  in  the  vasculature  of  the 
tumors  observed  using  a  new  dual  wavelength  homodyne  near-infrared  (NIR)  system. 


Methods:  Breast  13762  NF  adenocarcinomas  were  implanted  in  skin  pedicles  on  the  forebacks  of 

adult  female  Fischer  rats.  Once  the  tumors  reached  ~lcm  diameter,  the  tumor  blood  SO2  was  assessed 
by  NIR  spectroscopy  using  a  dual  wavelength  NIR  system  (758  nm  and  782  nm)  in  transmission 
geometry  [2].  The  tumor  blood  volume  and  SO2  were  calculated  from  the  light  amplitude.  The  rats 
were  maintained  under  general  gaseous  anesthesia  (33%  O2,  66%  N2O  and  0.5%  methoxyflurane). 
Once  stable  baseline  measurements  were  achieved,  the  inhaled  gas  was  altered  to  pure  oxygen  or 
carbogen  and  dynamic  changes  were  observed  over  a  period  of  two  hours.  Both  the  magnitude  and  rate 
of  change  of  SO2  were  examined.  Following  the  NIR  experiments,  40  pi  HFB  were  injected  directly 
into  both  central  and  peripheral  regions  of  the  tumors.  A  tunable  2  cm  *H/19F  single  turn  solenoid  coil 
was  placed  around  the  tumor  and  MR  experiments  were  performed  using  a  4.7  T  magnet.  Regional 
tumor  pOz  was  estimated  using  the  relationship:  /?02(torr)  =  [R1  —  0.074]/ 0.00 16 ,  where  Rl  is  the 

spin  lattice  relaxation  rate  of  HFB.  Twenty-three  p(h  maps  were  produced  in  3  hours  with  respect  to 
respiratory  challenge. 


Results:  NIR  showed  significant  changes  in  vascular  oxygenation  accompanying  respiratory 

interventions.  19F  MR-EPI  also  showed  significant  changes  in  tissue  pOi,  with  considerable  regional 
heterogeneity  in  both  absolute  values  and  rate  of  change  accompanying  interventions.  Generally, 
changes  in  vascular  SO2  preceded  tissue  pOi,  particularly  for  smaller  tumors. 


Discussion:  Regional  tumor  /?02and  blood  SO2  are  important  physiological  parameters.  The 

capability  to  measure  them  will  provide  insight  into  progressive  physiological  changes  in  a  tumor 
accompanying  interventions.  NIR  has  the  advantage  of  being  entirely  noninvasive,  but  the  MRI 
approach  clearly  reveals  detailed  oxygenation  heterogeneity.  We  believe  that  better  understanding  and 
monitoring  of  tumor  oxygenation  can  lead  to  improved  tumor  therapy. 
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Molecular  Determinants  of  Sensitivity  to  Anti  tumor  Agents 

Tumor  oxygenation  and  measurement  of  regional  dynamic  changes 
Ralph  P.  Mason.  Sandecp  Hunjan,  Anca  Constantinescu,  Yulin  Song,  Eric  W.  Hahn,  and  Peter  P. 
Andch,  Advanced  Radiological  Sciences,  U.T.  Southwestern  Medical  Center,  Dallas.  Texas  and 
Christian  Blum  and  Peter  Peschke,  Deutsches  Krebsforschungszentrum,  Heidelberg,  Germany 

Therapeutic  efficacy  may  be  influenced  by  tumor  oxygenation.  In  particular,  hypoxic 
tumors  resist  radiotherapy  and  may  be  good  candidates  for  hypoxia  selective  cytotoxic  agents. 
We  recently  described  a  novel  approach  to  measuring  regional  tumor  oxygen  tension  using 
pulse  burst  saturation  recovery  (PBSR)  nuclear  magnetic  resonance  (NMR)  echo  planar  imaging 
(EPI)  relaxometry  of  hexafluorobenzene  (HFB)  (1).  We  have  now  compared  oxygen  tension 
measurements  in  a  group  of  size  matched  Dunning  prostate  rat  tumors  R3327-AT1  made  using 
this  new  method  with  a  traditional  polarographic  method:  the  Eppendorf  Histograph.  Wc  also 
demonstrate  extension  of  the  MR  techniques  to  rat  breast  tumors. 

M&lhp.ds:  Rat  Dunning  prostate  R3327-AT1  or  breast  13762  NF  adenocarcinomas  were 
examined  at  a  volume  <2  cm3  or  >  3,5  cm3:  for  MRI  40  jil  HFB  were  injected  directly  into  the 
tumor,  both  centrally  and  peripherally.  The  rat  was  maintained  under  general  gaseous  anesthesia 
(33%  O2, 66%  N2O  and  0.5%  methoxyflurane).  A  tunable  2  cm  ^H/^F  single  turn  solenoid  coil 
was  placed  around  the  tumor  and  MR  experiments  were  performed  using  a  4.7  T  magnet 
equipped  with  actively  shielded  gradients.  Tumor  oxygenation  was  assessed  using  PBSR- 
EPI  of  HFB.  By  applying  the  acquisition  protocol  ARDVARC  (Alternated  Relaxation  Delays 
with  Variable  Acquisitions  to  Reduce  Clearance  effects),  we  achieved  R1  maps  in  8  min.  A 
series  of  maps  was  acquired  over  a  period  of  2  hours  with  respect  to  respirutory  challenges.  pC>2 
maps  were  then  generated  by  applying  the  relationship:  pC>2  (tore)  =  (R1  -0.074)/0.0016.  In 
parallel  experiments  p02  was  determined  polarographically 

Results:  Similar  oxygen  tension  distributions  were  found  using  MRI  or  polarography  and 
both  techniques  showed  that  tumors  with  volume  >  3.5  cm3  wcrc  significantly  (p  <  0.0001)  less 
well  oxygenated  than  smaller  tumors  (volume  <  2  cm3),  Using  the  ppj  approach  we  also 
examined  response  to  respiratory  challenge.  Increasing  the  concentration  of  inspired  oxygen 
from  33%  to  100%  O2  produced  a  significant  increase  (p  <  0.0001)  in  tumor  oxygenation  for  a 
group  of  small  tumors.  In  contrast  no  change  was  observed  in  the  mean  p02  for  a  group  of  large 
tumors.  Consideration  of  individual  tumor  regions,  irrespective  of  tumor  size  showed  a  strong 
correlation  between  the  maximum  pC>2  observed  when  breathing  100%  O2.  as  compared  with 
mean  baseline  pC>2. 

Conclusions:  These  results  further  demonstrate  the  usefulness  of  ^F  EPI  to  assess  changes  in 
regional  tumor  oxygenation.  The  ability  to  measure  p02  could  be  valuable  in  pre-clinical 

evaluation  of  novel  therapies  and  ultimately  allow  therapy  to  be  individualized  and  optimized  for 
patients. 

This  work  was  supported  in  part  by  The  American  Cancer  Sociery  (RPM),  The  Whitaker 
Foundation  (RPM),  DOD  Breast  Cancer  Initiative  (YS),  Verein  zur  Ffirderung  der 
Krebserkennung  und  Krebshandlung  e.V.  Heidelberg  (PP)  and  the  NIH  BRTP  Facility  #5-P41- 
RR02584.  7 
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INTRODUCTION  Tumor  oxygenation  has  a  critical  influence  on  success  in 
radiotherapy  and  is  thought  to  play  a  key  role  in  angiogenesis  and  metastasis.  However,  no 
non-invasive  procedure  for  tumor  oximetry  has  been  established  in  clinical  practice,  so  far. 
Here,  we  present  and  compare  two  oximetry  techniques:  the  FREDOM  (Fluorocarbon 
Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  approach  based  on 
hcxafluorobenzene  (HFB)  to  measure  tumor  tissue  oxygen  tension  (p02 )  and  near  infrared 
(NIR)  spectroscopy  based  on  a  new  I/Q  system  to  measure  changes  in  tumor  vascular 
hemoglobin  saturation  (5O2)  and  concentration  [Hb], 

METHODS  NF  13762  breast  adenocarcinomas  were  implanted  in  skin  pedicles  on  the 
forebacks  of  adult  female  Fischer  rats.  Once  the  tumors  reached  -1  cm  diameter,  the 
dynamic  characteristics  of  tumor  vascular  5O2  in  response  to  respiratory  challenge  were 
monitored  by  NIR  spectroscopy  Then,  40  p.1  HFB  were  injected  directly  into  the  tumor.  3-D 
MR  spin-echo  H  images  were  acquired  for  anatomical  reference  and  corresponding  19F 
images  were  then  obtained  to  reveal  the  distribution  of  HFB  in  the  tumor.  The  dynamic 
characteristics  of  regional  tumor  p02  were  assessed  using  FREDOM. 

RESULTS  Tumor  vascular  y02  increased  rapidly  after  a  gas  switch  from  baseline  (33% 
O2)  to  either  carbogen  or  100%  O2  and  then  slowly  returned  to  baseline  after  the  gas  was 
switched  back  to  baseline.  Total  hemoglobin  change  was  insignificant.  Tumor  tissue  p02 
also  showed  significant  changes,  but  with  considerable  regional  heterogeneity  in  both 
absolute  values  and  rate  of  change.  Changes  in  ^02  preceded  those  in  p02  Tumor  voxels 
with  high  baseline  p02  had  significantly  different  response  characteristics  from  those  with 
initially  low  p02 ,  with  voxels  of  high  baseline  p02  showing  significant  changes  in  p02 , 
while  voxels  of  low  baseline  p02  showing  small  changes.  Strong  correlation  existed 
between  the  maximum  p02  value  attained  during  the  course  of  an  experiment  and  mean 
baseline  p02 

CONCLUSIONS  Synergistic  application  of  the  FREDOM  and  NIR  techniques  will 
provide  new  insight  into  issues  of  tumor  angiogenesis  and  perfusion,  and  could  lead  to 
improved  tumor  therapy. 
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DIVERSE  APPROACHES  TO  MONITORING  OXYGEN  DYNAMICS  IN  RAT  BREAST  AND 
PROSTATE  TUMORS 

Dawen  Zhao,  Yulin  Song,  Hanli  Liu+,  Anca  Constantinescu,  Eric  W.  Hahn,  and  Ralph  P. 
~  son’  Department  of  Radiology,  UT  Southwestern,  Dallas,  TX  75390  and  *BME  UT- 
Arlington,  TX76109 

Summary:  The  increasing  evidence  that  direct  measurement  of  p02  in  patient's  tumors  has 
prognostic  value  provides  strong  impetus  to  develop  robust  methods  for  measuring  tumor 
oxygen  dynamics.  Using  the  FREDOM  NMR  approach  (Fluorocarbon  Relaxometry  using  Echo 
planar  imaging  for  Dynamic  Oxygen  Mapping)  we  are  able  to  examine  regional  variations  in 
tumor  p02  in  response  to  interventions.  To  obtain  a  more  complete  view  of  vascular  dynamics 
we  now  incorporate  Near  Infrared  (NIR)  measurements  of  hemoglobin  concentration  [Hb]  and 
saturation  (sHb02)  in  tumors. 

Methods:  Dunning  prostate  R3327  (ATI  and  HI)  and  13762NF  breast  tumors  were  implanted  in 
pedicles  (1)  on  the  foreback  of  syngeneic  Copenhagen  and  Fisher  rats,  respectively.  Once 
tumors  reached  ~  1  cm  diameter,  rats  were  anesthetized  (generally,  1.2%  isoflurane  in  air)  and 
examined  using  a  homodyne  dual  wavelength  NIR  device  in  transmission  mode,  to  monitor 
relative  A[Hb]  and  AsHb02  (2).  A  pulse  oximeter  simultaneously  provided  sa02  in  the  leg. 
Variations  in  these  values  were  observed  with  respect  to  respiratory  challenge  and  vasoactive 
.  ru^s  (oxy9®n.  carbogen  and  hydralazine).  The  following  day  the  same  tumors  were 
interrogated  by  F  NMR  EPI  with  respect  to  the  same  interventions  (3).  Hexafluorobenzene  (45 
jil)  was  injected  directly  into  both  central  and  peripheral  regions  of  the  tumors  using  a  fine  sharp 
needle  (32G)  and  p02  maps  determined  with  8  minute  time  resolution  and  1.25  mm  in  plane 

SIn?f'U!IOn  USi°9  pulse  burst  saturation  NMR  echo  planar  imaging  relaxometry 
{FREDOM).  In  some  instances,  additional  studies  were  performed  using  oxygen 
microelectrodes  and  the  OxyLite™  optical  fiber  sensor. 

Results:  As  expected,  the  fastest  changes  in  oxygenation  in  response  to  respiratory  challenge 
occurred  in  the  arterial  saturation  (sa02)  in  the  leg  with  a  typical  time  constant  x<20  s  (Fig.  1). 

Small  tumors  often  showed  biphasic  vascular  response  with  a  rapid  AHb02  component  (T<30  s) 
approaching  that  of  sa02  and  a  more  sluggish  component  increasing  over  20  mins.  The  slow 
component  alone  was  typical  of  larger  tumors.  Changes  in  A[Hb]  indicated  that  vascular  volume 
was  also  modulated.  Changes  in  tumor  tissue  p02  were  highly  variable.  Well  oxygenated 

regions  responded  significantly  with  a  time  constant  in  the  range  x=  10-30  mins.  The 
response  of  initially  poorly  oxygenated  regions  was  highly  tumor  dependent.  ATI  tumor  regions 
initially  <10  torn  showed  little  response,  whereas  most  regions  in  the  slower  growing  HI 
responded,  significantly  reducing  the  hypoxic  fraction  with  oxygen  or  carbogen  (see 
accompanying  abstract  by  Zhao  et  a!.).  Although  the  breast  tumors  grow  rapidly,  hypoxic 
regions  were  also  found  to  respond  in  these  tumors  (Fig.  2).  Both  oxygen  microelectrodes  and 
tne  OxyUte  showed  changes  in  p02,  which  were  comparable  to  the  FREDOM  results  in  terms 
of  both  rates  and  magnitude  of  change. 

Conclusions:  The  FREDOM  approach  allows  dynamic  changes  in  regional  tumor  p02  to  be 
m  'versf  ^umors-  These  data  indicate  that  regional  response  to  respiratory  challenge 
depends  not  only  on  the  baseline  tumor  oxygenation,  but  also  the  degree  of  tumor 
differentiation.  By  adding  NIR,  we  are  additionally  able  to  probe  changes  in  vascular  volume  and 
hemoglobin  oxygen  saturation.  These  preliminary  data  were  obtained  sequentially,  but 
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simultaneous  NIR  and  FREDOM  investigations  are  imminently  feasible  to  provide  deeper  insight 
into  tumor  vascular  phenomena.  We  believe  the  complimentary  applications  of  the  two 
techniques  will  provide  deeper  insight  into  tumor  physiology  and  mechanisms  of  modulati  g 
tumor  physiology  for  therapeutic  enhancement. 


Fig.  1  Dynamic  response  in  sa02, 
AHbOZl  and  AHbTotai  in  4.5-cm3  rat 
breast  tumor  with  respect  to 
respiratory  challenge  observed 
using  NIR  instrument.  A  bi¬ 
exponential  curve  fit  for  changes 
in  AHb02  gave  time  constants 

Tt=1  1  s  and  i2  =  27.8  mins, 
respectively. 


Fig.  2  Dynamic  changes  in  p02 
measured  using  FREDOM  of 
representative  regions  in  a  second 
breast  tumor  (3.5  cm3).  The  local 
response  depended  strongly  on  initial 
baseline  p02,  both  in  terms  of  rate 
and  magnitude  (measurements  at  8 
min  intervals). 
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FREDOM  (Fluorine  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen 
Mapping):  a  contextual  review. 

Ralph  P.  Mason,  Dawen  Zhao,  Anca  Constantinescu,  Yulin  Song,  Lan  Jiang  &  Eric  W. 
Hahn,  Department  of  Radiology,  UT  Southwestern  Medical  Center,  Dallas,  Texas  75390 

Tumor  oxygenation  has  been  widely  recognized  as  a  potent  factor,  which  influences  tumor 
response  to  therapy.  Recent  development  of  methods  to  measure  p02  in  the  clinic  has  generated 
data  confirming  a  prognostic  relationship  between  tumor  pC>2  and  clinical  outcome.  Such  results 
provide  a  powerful  stimulus  for  development  of  enhanced  methods  for  monitoring  tumor 
oxygenation.  Many  diverse  techniques  exist  with  characteristic  virtues,  providing  unique  insight, 
as  represented  by  this  session  on  Tumor  Oximetry.  We  have  developed  an  approach  using  19F 
NMR  echo  planar  relaxometry  of  hexafluorobenzene  (HFB)  following  direct  intra  tumoral 
administration.  Currently,  this  provides  maps  of  p02  in  eight  minutes  with  a  typical  precision  of 
1-3  torr  at  30  -  100  locations  within  a  plane  of  an  individual  tumor.  We  will  describe 
applications  to  interrogate  p02  dynamics  in  response  to  a  variety  of  interventions,  such  as 
respiratory  challenge  and  vasoactive  drugs  in  rat  prostate  and  breast  tumors  and  human 
lymphoma. 

FREDOM  exploits  the  well  known  sensitivity  of  the  spin  lattice  relaxation  rate,  Rl,  of 
perfluorocarbons  to  oxygen.  The  choice  of  HFB  as  the  interrogating  molecule  reflects  its 
remarkable  characteristics:  single  resonance,  which  is  exceedingly  sensitive  to  changes  in  p02 
[p02  (torr)=  (Rl-0.0835)/0.001876  at  37  °C],  while  being  insensitive  to  changes  in  temperature 
[at  5  torr  and  37  °C,  a  2  °C  error  in  temperature  estimate  would  yield  An  error  in  p02  =??]. 
Hexafluorobenzene  is  readily  available,  cheap  and  has  well  documented  lack  of  toxicity. 

Correlative  studies  show  that  results  achieved  using  FREDOM  or  the  Eppendorf  electrode 
system  provide  a  similar  representation  of  a  tumor  (1).  Investigation  of  dynamic  changes  at 
specific  locations  show  similarity  with  polarographic  electrodes  or  fiber  optic  probes.  Most 
tumors  exhibit  great  heterogeneity  in  baseline  oxygenation  with  regions  ranging  from  relatively 
hypoxic  (<  5  torr)  to  well  oxygenated  (>  20  torr).  Most  tumors  show  reduced  baseline  pC>2  at 
larger  size,  e.g.,  in  the  anaplastic  Dunning  prostate  R3327-AT1  subline  (Volume  doubling  time 
(VDT)  ~  5  days)  tumors  with  a  vol.  >  3.5  cm3  were  significantly  (p  <  0.0001)  less  well 
oxygenated  than  tumors  (<  2  cm3).  Similar  results  were  found  in  the  slower  growing  moderately 
well  differentiated  R3327-HI  subline  (VDT  ~  10  days).  Changes  in  pC>2  in  response  to 
respiratory  challenge  with  oxygen  or  carbogen  were  found  to  be  dependent  on  baseline  pC>2  in 
the  ATI  tumor.  However,  by  contrast,  even  the  poorly  oxygenated  areas  of  the  HI  tumors  were 
very  sensitive,  and  responded  with  highly  significant  transient  increases  in  pC>2.  These  findings, 
if  confirmed,  in  human  tumors  could  be  valuable  in  treatment  planning. 

We  believe  the  ability  to  probe  oxygen  tension  dynamics  quantitatively  at  multiple 
selected  locations  of  interest  within  a  tumor  will  be  valuable  in  probing  tumor  development  and 
the  efficacy  of  adjuvant  interventions  to  modulate  tumor  physiology. 
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INTRODUCTION 

Tumor  oxygenation  has  a  profound  effect  on  growth 
and  development,  and  hypoxia  reduces  radiosen-sitivity. 
Moreover,  increasing  evidence  from  clinical  trials  has  indicated 
that  poorly  oxygenated  tumors  have  poor  prognosis.  Therefore, 
development  of  non-in vasive  techniques  to  accurately  measure 
oxygenation  is  essential  in  cancer  treatment  planning  and 
prognosis  prediction.  Here,  we  present  and  compare  two  such 
techniques:  the  FREDOM  (Fluorocarbon  Relaxometry  using 
Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  approach 
based  on  hexafluorobenzene  (HFB)[1]  to  measure  tumor  tissue 
oxygen  tension  (p02)  and  NIR  spectroscopy  to  measure  changes 
in  tumor  vascular  hemoglobin  saturation  (s02)  and 
concentration  [Hb],  The  synergistic  application  of  MR  and  NIR 
techniques  could  provide  new  insight  into  issues  of  tumor 
angiogenesis  and  perfusion. 

METHODS 

NF  13762  breast  and  Dunning  prostate  R3327-AT1 
adenocarcinomas  were  implanted  in  skin  pedicles  on  the 
forebacks  of  adult  female  Fischer  and  male  Copenhagen  rats 
(-250  g),  respectively.  Once  the  tumors  reached  -1  cm  diameter, 
the  rats  were  anesthetized  with  200  pi  ketamine  hydrochloride 
(100  mg/ml)  and  maintained  under  general  gaseous  anesthesia 
(33%  02,  66%  N20  and  0.5%  methoxyflurane).  The  tumor 
vascular  s02  was  assessed  by  NIR  spectroscopy  using  a  new 
dual  wavelength,  homodyne  system  (wavelengths  758  nm  and 
782  nm)[2],  while  inhaled  gas  was  alternated  between  33%  02, 
carbogen  (95%  02  +  5%  CQ2),  and  100%  02.  The  light 
amplitude  and  phase  changes  caused  by  the  tumor  are  related  to 
changes  in  hemoglobin  concentration  [Hb]  and  hemoglobin 
saturation  [Hb02],  i.e.,  sfy: 

A[Hb]toiii  =  -  [3.56*  log  (A,/Ac)758 

+  8.59*  log  (A/Ac)  j/d  (1) 

A[Hb02]  -  A[Hb]  =  -  [18.27*  log  (A,/Ac) 7S8 

+  20.92  *  log  (A]/ Ac)  782]/d  (2) 

where  Ai  is  the  initial  amplitude  (amplitude  of  baseline),  Ac  the 
current  amplitude,  d  the  direct  source-detector  separation  in  cm, 
and  AO  the  change  in  concentration  in  mM. 

Following  the  NIR  experiments,  a  tunable  2  cm  1H/19F 
single  turn  solenoid  coil  was  placed  around  the  tumor  and  40  pi 
HFB  were  injected  directly  into  both  central  and  peripheral 
regions  of  the  tumor  using  a  32  G  needle.  3D  spin-echo  (SE)  *H 
images  were  acquired  for  anatomical  reference  and 
corresponding  l9F  images  were  then  obtained  to  show  the 
distribution  of  HFB  in  the  tumor.  Regional  tumor  p02  maps 
were  generated  using  l9F  PBSR-EPI  based  on  the 
relationship:  pO  2  (torr)  =  [R1  -  0.0836]  /  0.00 188 ,  where  R\ 
(1/T)  is  the  spin  lattice  relaxation  rate  of  HFB  in  1/sec. 


RESULTS 

Tumor  vascular  s02  increased  almost  immediately  after 
a  gas  switch  from  baseline  (33%  02)  to  either  carbogen  or  100% 
02  and  increased  steadily  for  several  minutes,  and  then  gradually 
returned  to  baseline  after  the  gas  was  switched  back  to  baseline. 
In  contrast,  total  hemoglobin  change  was  insignificant, 
indicating  relatively  constant  blood  volume  in  the  tumor.  Tumor 
tissue  p02  also  showed  significant  changes,  but  with 
considerable  regional  heterogeneity  in  both  absolute  values  and 
rate  of  change.  Both  carbogen  and  100%  02  inhalation  produced 
significant  changes  in  ,s02  and  pC h,  especially  with  100%  02 
inhalation.  Temporal  dynamic  response  in  both  s02  and  p02 
were  modeled  and  exponential  time  constants  were  determined. 
The  FREDOM  technique  also  allowed  us  to  compute  the  time 
constant  on  a  voxel-by-voxel  basis.  It  was  found  that  s02  had  a 
faster  time  constant  than  p02,  especially  in  the  cases  of  larger 
tumors,  which  were  found  to  be  less  well  oxygenated  and 
presumably  less  well  perfused.  It  was  also  found  that  some 
tumors  showed  a  bimodal  response  (slow  plus  fast)  in  s02,  as 
compared  to  a  unimodal  response  (slow)  in  p02. 

DISCUSSION 

An  increase  in  inspired  gas  F02  should  lead  to 
increased  tumor  vascular  s02,  and  hence,  increased  tumor  tissue 
p02.  Our  data  indicate  that  breathing  elevated  02  did  indeed 
have  a  significant  effect  on  both  tumor  vascular  s02  and  tissue 
p02.  Vascular  s02  values  were  found  to  have  a  faster  time 
constant  than  tumor  tissue  p02,  with  greater  differences  in  large 
tumors.  This  probably  reflects  the  extensive  perfusion  of  the 
small  tumors  with  lesser  perfusion  of  large  tumors  since  time 
constant  should  be  inversely  proportional  to  tumor  blood  flow 
(TBF),  as  reflected  by  lower  mean  p02  and  larger  hypoxic 
fraction.  We  believe  that  application  of  multiple  approaches  to 
tumor  oxygenation  can  lead  to  better  understanding  of  tumor 
physiology  and  probably  optimized  tumor  therapy. 
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